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The  spectrum  of  the  stimulated  far  infrared  emission  from  p-Germanium  in 
crossed  electric  and  magnetic  fields  is  studied  by  means  of  a  tunable  narrow- 
band  GaAs-detector .  A  multimode  spectrum  is  observed  from  polished  high 
parallel  samples.  The  mode  spectrum  is  quantitatively  explained  in  terms  of 
internal  waveguide-like  modes. 

With  the  help  of  an  external  resonator  the  number  of  modes  was  drastically 
reduced.  The  external  resonator  has  a  confocal  configuration.  For  the  first 
time  a  single  mode  operation  is  demonstrated. 

The  problem  of  the  electron  distribution  in  crossed  electric  and  quantizing 
magnetic  dields  is  treated  for  free  electrons  in  a  parabolic  band.  The  master 
equation  for  the  diagonal  elements  of  the  density  matrix  is  solved  including 
the  kz~dependence  of  the  electron  distribution.  It  is  shown  that  population 
inversion  between  Landau  levels  can  be  achieved  if  scattering  by  optical 
phonons  and  a  quasi-elastic  scattering  process  are  taken  into  account.  -V4— — 
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The  energy  relaxation  of  2D  electrons  in  GaAs/GaAlAs  structures  has  been 'in¬ 
vestigated  by  analysing  the  electric  field  dependence  of  Shubnikov-de  Haas 
oscillations,  the  far  infrared  emission  and  photoluminescence  spectra.  A  quite 
general  behavior  of  the  electron  heating  AT  =  Te  -  TL  as  a  function  of  the 
input  power  per  electron  P0  is  found:  T«/Pe.  The  corresponding  energy  re¬ 
laxation  times  in  the  range  of  nsec  are  independent  of  the  electron  temperature 
up  to  30  K  and  inversly  proportional  to  the  electron  density.  At  higher  elec¬ 
tron  temperatures  the  energy  relaxation  is  governed  by  optical  phonon  emission. 
However,  the  onset  depends  on  electron  concentration  and  is  different  for 
heterostructures  and  quantum  wells.  From  intensity  dependent  cyclotron  resonanci 
transmission  experiments  Landau  level  lifetimes  between  0.2  ns  and  1  ns  depen¬ 
ding  on  the  electron  density  are  found  in  agreement  with  data  from  time- 
resolved  photoluminescence. 

Smith  Purcell  radiation  is  coupled  out  from  a  drifting  2D-carrier  gas  with 
the  help  of  sinusoidal  gratings.  From  the  emitted  radiation  the  form  of  the 
distribution  function  can  be  determined.  As  a  result  a  drifted  Fermi-distribu- 
tion  with  rather  low  temperatures  is  found  at  the  saturated  drift  velocity 
of  2  x  107cm/s. 
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a)  A  detailed  investigation  of  various  emission  processes  in  the  far  infrared 
with  the  technical  goal  to  develope  a  narrowband  and  powerful  solid  state 
source . 

b)  The  study  of  the  hot  electron  emission  spectra  in  MOS-type  devices  to 
obtain  information  on  the  electron  distribution  function  in  the  2D  channel 

Several  methods  and  techniques  have  been  suggested  in  the  proposal.  Not  all  of 
them  have  been  persued  in  detail.  However,  we  can  state  that  both  claims  have 
been  achieved  within  the  project. 

The  final  results  are  the  following: 

We  were  able  to  repeat  as  the  first  group  in  the  western  world  the  findings  by 
a  Russian  group  (Andronov  et  al.):  the  generation  of  stimulated  emission  from 
streaming  carriers  in  p-Ge.  (Xir  main  contribution  to  this  field  was  the 
explanation  of  the  multimode  structure  of  the  lasing  process  and  the  first 
demonstration  of  a  single  mode  operation  with  external  mirrors.  In  addition  a 
fundamental  theoretical  study  was  perf(  rned  to  explain  the  electron 
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distribution  under  streaming  conditions  for  a  single  band  case. 

The  also  suggested  method  of  generating  far  infrared  emission  from 
2-dimensional  plasma  oscillations  mas  not  continued  further  due  to 
requirements  which  we  could  not  fulfill.  It  would  have  been  necessary  to 
produce  gratings  with  periods  below  1000  A  which  is  beyond  the  possibilities 
of  our  technology.  However,  we  have  done  an  Important  step  forward  in  this 
direction  by  demonstrating  for  the  first  time  Smith-Purcell  radiation  from  a 
drifting  carrier  gas  under  a  periodic  grating.  The  continuation  of  this  work 
has  the  potential  to  produce  Bloch-type  oscillations  in  the  future. 

The  analysis  of  the  far  Infrared  emission  from  MOS-type  devices  was  restricted 
to  GaAs/GaAlAs  heterostructures  since  this  system  is  considered  as  the 
candidate  for  the  highest  speed  devices.  Ve  have  studied  the  energy  loss  rate 
in  GaAs/GaAlAs  heterostructures  as  a  function  of  the  electrical  input  power  to 
get  information  on  the  dominant  scattering  mechanisms  at  low  temperatures .  As 
a  result  from  an  analysis  of  the  black  body  emission  we  found  a  sample 
independent  energy  loss  rate  over  the  input  power.  This  is  only  valid  for 
electron  temperatures  up  to  40  K  where  acoustic  deformation  potential 
scattering  is  dominant.  For  higher  electron  temperatures  optical  phonon 
emission  dominates  the  energy  loss.  However,  there  is  still  some  controversy 
left  about  the  importance  of  hot  phonon  effects  in  the  electron  temperature 
range  between  50  K  and  100  K.  The  results  from  the  black  body  emission 
indicate  the  influence  of  hot  phonon  effects  while  our  most  recent  work  using 
the  analysis  of  the  Smith-Purcell  radiation  indicates  considerably  lower 
electron  temperatures  and  no  evidence  for  hot  phonon  effects. 


On*  of  the  eain  achievements  of  the  project  is  the  observation  of  the 
Seltb-Purcell  radiation  from  drifting,  carriers  over  a  periodic  grating  as 
mentioned  before.  Ve  have  used  the  emitted  spectrum  to  determine  the  hot 
electron  distribution  function  and  the  drift  velocity  in  extremely  high 
mobility  GaAs/GaAlAs  heterostructures.  This  is  the  first  measurement  of  the 
drift  velocity  in  a  semiconductor  by  a  purely  optical -technique.  The  observed 
drift  velocities  are  about  10  X  lower  than  derived  from  current  measurements. 

The  electron  temperatures  at  drift  velocities  close  to  the  saturation  velocity 
are  about  50  K.  and  thus  considerably  lower  than  obtained  from  black  body 

f 

emission  experiments  (~  100  K).  As  a  result  we  find  a  drifted  Fermi -Dirac 

I 

distribution  function  with  a  large  drift  energy  and  low  electron  temperature.  J 

( 

The  final  report  will  consist  of  two  parts: 

I 

I 

In  Fart  I  the  results  on  stimulated  emission  from  p-Ge  will  be  presented  j 

together  with  a  quantum  mechanical  calculation  of  the  hot  carrier  distribution 
function  under  crossed  field  conditions. 

In  Fart  II  the  far  infrared  black  body  emission  from  GaAs/GaAlAs 
heterostructures  will  be  discussed  together  with  the  analysis  of  the 
Smith-Purcell  radiation  from  drifting  carriers  under  a  periodic  potential.  j 
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Introduction 


The  distribution  function  in  momentum  space  is  of  particular  interest  in 

the  investigation  of  hot  carriers  in  external  fields.  This  widely  used 

expression  ’hot  carrier*  implicitly  assumes  that  the  distribution  function  is 

of  quasi  -  Maxwellian  type,  characterized  by  an  electron  temperature  higher 

than  that  of  the  thermal  bath.  This  model  corresponds  to  a  quite  Isotropic 

distribution  function  and  is  justified  if  scattering  events  can  be  regarded  as 

quasi  elastic  processes.  On  the  other  hand,  a  quite  different  type  of  high 

field  phenomena  characterelzed  by  a  highly  anisotropic  distribution  can  be 

obtained,  when  optical  phonon  emission  predominates  over  all  other  scattering 

mechanisms.  In  fact,  in  pure  semiconductors  at  low  lattice  temperatures  it  is 

possible  to  accelerate  a  carrier  up  to  an  energy  equal  to  the  optical  phonon 

energy  Provided  that  the  optical  phonon  interaction  is  adequately 

strong,  the  carriers  having  reached  will  emit  almost  immediately  an 

optical  phonon  by  dissipating  all  of  its  kinetic  energy.  The  carrier,  thereby 

scattered  near  to  the  ground  state,  will  be  accelerated  again  to  t=ftu>  and 

op 

will  repeat  thereafter  the  same  process  until  this  is  Interrupted  by  some 
other  scattering  event.  Such  a  repeated  motion  is  called  ’streaming  motion’ 
and  it  results  in  a  non  -  Maxwellian  type  distribution  function  /l/. 

By  applying  a  magnetic  field  perpendicular  to  the  electric  field  a 


k 
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carrier  accumulation  in  a  limited  area  of  the  momentum  region  can  occur,  which 
leads  to  a  population  inversion  /2 /.  For  such  a  hot  carrier  system, 
fascinating  effects  are  theoretically  predicted,  such  as  negative  differential 
conductivity  in  d.c.  fields  and  amplification  in  the  range  from  im-wave  to  far 
infrared. 73/ 

Several  experimental  investigations  of  the  streaming  effect  have  been 
performed  in  different  materials  /4/.  Hall  and  current  measurements  showed  a 
clear  evidence  for  streaming  motion  and  for  carrier  accumulation  / 5 /.  In  the 
p-Ge  system  it  was  possible  to  observe  far  infrared  emission  /6/.  This  opened 
up  the  possibility  of  directly  probing  the  distribution  function  of  the 
population  inverted  hot  carriers,  but  also  encouraged  further  seeking  for  far 
infrared  amplification. 

The  first  report  on  'laser  like1  from  p-Ge  was  reported  by  a  Russian 
group  /7/.  Shortly  after  this  report  also  Komiyama  et  al.  /8/  observed 
stimulated  emission  from  p-Ge.  The  observed  spectra  showed  a  broad  band 
emission  without  structure  but  with  a  clear  threshold  behavior.  The  most 
recent  investigations  were  devoted  to  the  study  of  the  real  mode  structure  of 
this  new  laser  process  /9/./10 /.  At  present  several  groups  are  persuing  the 
extraction  of  one  single  narrow  line  from  the  broad  band  emission  spectrum. 
Since  only  a  single  line  laser  can  be  used  for  spectroscopy. 

In  the  first  two  sections  of  this  report  a  model  for  the  explanation  of 
the  population  inversion  in  p-Ge  within  the  streaming  motion  picture  will  be 
discussed.  The  experimental  methods  are  explained  in  the  following  sections. 
In  the  final  sections  the  results  are  discussed  and  our  recent  explanation  of 
the  mode  structure  is  presented. 


Streaming  motion  is  expected  under  the  condition  that  the  optical  phonon 
Interaction  predominates  over  all  other  scattering  processes.  We  will  show 
that  such  a  situation  is  realized  in  pure  crystals  at  low  temperatures.  In 
ionic  crystals  the  optical  phonon  scattering  is  determined  by  very  strong 
polar  optical  coupling.  In  covalent  materials  .such  as  Ge.  the  carriers 
interact  with  the  optical  deformation  potential,  which  is  expressed  in  the 
form 


Q 

where  the  coupling  strength  DtK  is  9*10  eV/cm.  the,  density  of  the  crystal  p 

3  M 

is  5.32  g/cm  .  the  effective  mass  m  =  0.043*mQ  for  the  ligth  holes  and 

m*  =  0.35*iDq  for  the  heavy  holes.  The  optical  phonon  energy  is  =  37  meV. 

Figure  1  shows  the  scattering  probabilities  for  a  heavy  hole  in  p-Ge  at  4.2  K. 

The  interaction  is  not  so  strong  as  in  the  polar  optical  case,  but  the 

probability  of  optical  phonon  emission  at  kinetic  energies  equal  to  the 

optical  phonon  energy  is  still  much  higher  than  that  due  to  ionized  impurity 

scattering  and  acoustical  phonon  scattering.  N  Is  supposed  to  be  very  low  at 

4.2  K  and  therefore  only  the  second  term  (emission  term)  in  eq.l  contributes 

to  the  rate.  Carrier-carrier  scattering  is  supposed  to  be  unimportant  at 

14  3 

carrier  densities  lower  than  5*10  cm  . 

Let  us  now  Imagine  that  a  free  carrier  with  an  effective  mass  os*  is  at 
rest  at  time  t  =  0.  With  an  electric  field  applied,  the  carrier  is  accelerated 


in  the  direction  of  the  electric  field  with  an  acceleration  rate  v  =  e* E/m*. 


The  velocity  of  the  carrier  reaches  the  critical  velocity  v 


J4n**v2  =  hu  (2) 

op  op  v  * 

which  corresponds  to  the  kinetic  energy  equal  to  the  optical  phonon  energy 

h<j  .  at  time  t  =  t 
op  op 

top  =  (2m»»«fiuop)J4/(e«E)  (3) 

Supposed  that  E  is  not  too  high  the  carrier  will  almost  immediately  emit 
an  optical  phonon  and  return  to  about  the  ground  state.  The  carrier  will 
perform  a  streaming  motion  with  emitting  optical  phonons  after  the  time 
interval  t  .  The  drift  velocity  of  an  ideally  streaming  carrier  is 

v .  =  J4v  (4) 

In  an  actual  experimental  situation,  streaming  motion  is  only  possible  if 
the  carrier  suffers  no  other  scatterings  before  reaching  t  .  Therefore  the 
condition 


imp  ac 

must  be  fullfilled.  where  Tlmp  and  are  scattering  rates  for  impurity  and 
acoustical  phonon  scattering,  averaged  for  energies  smaller  than  the  optical 
phonon  energy.  From  eg. 3  and  eq.4  the  condition  for  the  lower  limit  of  the 
electric  field  Enin  results.  On  the  other  hand,  to  justify  the  expression 
’immediatly  scattering*  the  electric  field  should  not  exceed  a  maximum  limit 
E _  ,  which  Komlyama  /ll/  tried  to  estimate  in  the  following  rough  measure 

fflBX 

t  =  t  (l.l«fi«  )  (6) 

op  op'  op'  '  ’ 

where  is  the  optical  phonon  scattering  probability  for  the  energy  e. 

Thus  streaming  motion  is  expected  in  the  range  E  .  <  E  <  E  .  Table  1  shows 

characteristic  quantities  for  streaming  motion  in  different  materials. 

When  a  transverse  magnetic  field  !?  =  (0.0.B)  is  applied  to  an  ideally 
streaming  carrier  in  an  electric  field  If  *  (E.0,0),  the  trajectory  of 


streaming  Is  curved  due  to  the  Lorentz  force.  The  free  motion  of  the  carriers 

is  described  by  the  classical  equation 

=«•(£+  vx3)  (7) 

The  motion  of  the  carriers  is  a  cyclotron  orbit  with  the  center  at 

C  =  (O.-E/B.O)  in  velocity  space  (see  Fig.  1).  The  motion  in  real  space  is  a 

‘cycloid  motion*  with  a  mean  drift  velocity  v^  =  -E/B  in  y-direction.  To  take 

into  account  for  the  interaction  with  the  optical  phonon  interaction  it  is 

useful  to  introduce  the  normalized  field  f 

£  =  vop.B/E  (8) 

which  specifies  the  position  of  the  center  of  the  cyclotron  orbit  relative  to 

the  circle  |v|  =  vQp.  At  this  circle  the  motion  of  the  carrier  is  Interrupted 

by  optical  phonon  emission  and  scattered  back  to  v  =  (0.0,0).  The  trajectory 

of  the  Ideally  streaming  carriers  is  called  the  'main  trajectory* 

If  £  <  1.  point  C  is 'situated  outside  the  circle  |v|  *  v  .  For  f  >  1. 

op 

the  point  C  lies  inside  the  circle.  For  £  >  2,  the  main  trajectory  is  closed 
within  the  circle  |v|  a  VQp  and  streaming  motion  becosies  prohibited.  In  the 
last  two  cases  a  round  area  is  formed  around  C  (Indicated  by  the  shaded  area 
in  Fig.  1),  in  which  the  trajectories  of  the  cyclotron  oscillation  do  not 
cross  the  circle  |v|  *  vQp:  A  carrier  in  this  area,  if  present,  performs  a 
sustained  cyclotron  oscillation  and  drifts  along  the  y-direction  in  real 
space.  In  the  case  of  ideally  streaming  for  1  <  £  <  2  this  'passive*  region  is 
expected  to  be  empty,  as  the  main  trajectory  is  outside  of  this  area.  In  a 
real  experiment,  however,  a  few  carriers  may  be  occasionally  scattered  into 
this  area  by  optical  phonon  emission  after  an  adequately  penetration  into  the 
high  energy  range  |v|  >  vQp.  Carriers,  once  scattered  into  this  region,  stay 
there  till  removed  by  other  scattering  processes;  as  a  consequence  a  carrier 
accumulation  is  formed  in  this  area.  In  an  early  theoretical  paper  Neada  and 


Kurosawa  /2 /  obtained  a  significant  accumulation  of  heavy  holes  in  this  region 
for  p-Ge . 


Population  Inversion  and  Far  Infrared  Emission 

In  p-type  Ge  the  ligth  and  heavy  lole  bands  are  degenerate  at  the  band 

edge.  In  the  case  of  a  carrier  accumulation  far  infrared  emission  can  be 

expected  in  a  wide  frequency  range (Fig.  2).  Under  the  special  condition 

f*1  <  1  <  f*  (h  for  heavy,  1  for  ligth  hole  band)  an  accumulation  takes  place 

only  in  ligth  hole  band  and  streaming  in  the  heavy  hole  band.  The  light  hole 

accumulation  is  addionally  Increased  by  scattering  Induced  redistribution  of 

carriers  from  the  heavy  hole  band.  This  accumulation  can  lead  to  a  population 

inversion  between  the  ligth  and  heavy  hole  bands  in  momentum  space.  In  thermal 

j  jj  1  h  3/2 

equilibrium  the  ratio  of  the  populations  is  n  /n  •  (a  /m  )  .  In  the  above 

situation  this  is  changed  to  n* /n1  *  (mVm*1)'^«(rVtkp),  where  r *  is  the  life 
time  of  the  holes  in  the  accumulation  area  and  tQjj  is  the  optical  phonon 
scattering  time  of  the  heavy  holes.  For  the  situation  <  1  <f*.  the  life 
time  of  the  ligth  holes  is  assumed  to  be  larger  them  the  heavy  hole  scattering 
time.  Therefore  a  population  inversion  is  built  up.  Far  infrared  emission  is 
expected  due  to  transitions  between  the  ligth  and  heavy  hole  bands.  Assuming 
parabolic  bands  the  maximum  energy  of  the  emission  is  given  when  the 

light  hole  energy  reaches  the  optical  phonon  energy. 

tMx  =  (1/2DB1  -  l/2mh)(k*p)aha  (9) 

<i/2*1><kiP>v'  -  o°) 


Fig.  1:  Classical  motion  of  carriers  in  crossed  fields  in  velocity  and 
real  space  for  different  values  of  f. 

Fig.  2-  Valence  bandstructure  of  Ge  and  drifted  distribution  functions 
of  heavy  and  light  holes. 


For  p-Ce  the  maximum  energy  fe  =32.5  meV.  This  type  of  population  inversion 

luBX 

is  of  particular  interest,  since  it  occurs  in  the  continuum  of  the  energy 
hands,  in  contrast  to  the  usual  population  inversion  realized  between  discrete 
energy  levels  in  many  existing  laser  systems. 

For  a  calculation  of  the  gain  for  far  infrared  amplification  it  is 
necessary  to  know  the  distrubtion  functions  of  both  carrier  types.  In  addition 
other  absorption  processes  have  to  be  considered,  e.g.  free  carrier 
absorption.  Pozhela  et  al./12/  used  Monte  Carlo  simulations  to  get  realistic 
distribution  functions  and  from  them,  he  calculated  the  emission  and 
reabsorption  to  get  a  gain  of  about  0.05  cm  *.  This  is  a  rather  low  gain.  This 
implies  a  high  quality  factor  of  the  resonator  to  achieve  stimulated  emission. 


The  crystals  of  p-Ge  used  in  this  work  are  Ga  doped  having  a 
concentration  p  *  |Na  -  Njj|  =  1.1*10**  cm*3.  The  resistivity  at  room 
temperature  is  25  0cm.  The  sample  parameters  have  to  be  choosen  to  achieve 
high  gain.  For  a  too  high  concentration  streaming  motion  becosws  inpossible  as 
impurity  scattering  will  become  dominant.  The  concentration  p  =  1.1*10**  cm-3 
seems  to  compromise  with  both  requirements.  The  samples  are  cut  from  the 
crystal  with  a  wire  saw.  cleaned  with  triclorethylene.  acetone  and  propanol 
and  finally  etched  with  CP4. 

The  sample  dimensions  are  about  5  x  1  x  0.5  mm3.  For  making  contacts  In 
is  evaporated  to  the  end  faces  and  thereafter  alloyed  at  400  °C  for  5  min. 

The  experiments  are  performed  in  a  helium  cryostat  with  two 
superconducting  magnets  with  the  sample  placed  in  the  middle  of  the  upper 
magnet.  For  the  current  and  spontaneous  emission  measurements  the  electric 


iifVXV  >TU  WJ  TU 


-14- 


field  is  appl ied  by  a  cable  discharge  pulser  with  50  Q  impedance  and  a  maximum 
output  voltage  of  1000  V.  The  pulselength  is  1  psec  to  avoid  sample  heating, 
the  repetition  rate  is  tunable  between  10  and  100  Hz.  The  current  through  the 
sample  is  measured  by  a  series  resistance,  which  is  also  immersed  in  liquid 
helium.  For  the  emission  experiments  the  far  infrared  is  guided  by  a  10  mm 
bore  brass  ligth  pipe  to  different  detectors,  which  are  placed  in  the  same 
cryostat  in  the  middle  of  the  lower  magnet.  A  Ga  doped  Ge-detector  is  used  as 
a  broad  band  detector  with  a  peak  sensitivity  around  100  pm.  The  spectra  of 
the  emission  are  investigated  with  a  GaAs-detector.  The  energy  of  the  Is  -*  2p« 
impurity  transition  is  linearly  tunable  in  a  wide  range  by  a  magnetic  field, 
the  linewith  is  0.25  cm'1  /13/.  An  InSb-cyclotron  resonance  detector  is 
addionally  used  for  spectroscopy.  All  detectors  are  driven  in  a  constant 
current  operation,  the  photoconductlve  signal  is  preamplified  and  thereafter 
Integrated  by  the  Boxcar  Integrator. 


For  stimulated  emission  the  sample  dimensions  have  to  be  larger  and  the 
side  faces  have  to  be  paralel 1  within  30*'.  This  processing  is  done  by  a 
comerclal  company  for  infrared  lenses.  Contacts  are  made  in  the  same  way  as 
described  above.  For  lowering  the  surface  resistance  and  to  secure  the 
homogenlty  of  the  electric  field  copper  plates  are  soldered  to  the  contact 
sides.  Copper  wires  are  connected  to  these  plates.  The  copper  plates  serve 
also  as  sample  holder  to  keep  the  sample  in  an  upright  position  and  in  the 
center  of  the  magnet.  The  angle  of  total  reflection  for  Ge  (n=4)  is  only  15°. 
As  a  result  a  resonator  is  formed  by  internal  total  reflections  at  the  sample 
surfaces  if  the  side  faces  are  polished.  The  laser  mode  will  be  similar  to 
wave  guide  like  mode.  The  sample  dimensions  have  to  be  quite  big 
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(50  x  5  x  5  on3)  for  this  type  of  laser.  As  a  consequence  the  sample  Impedance 

is  very  low  (about  211):  The  current  denislty  is  j  =  n^e'v^  a  200  A/cm2,  for 

n  ss  1.1»1014  cm"3  and  v^  =  1*107  cm/s.  For  a  contact  area  of  5  x  0.5  cm2  the 

current  through  the  sample  is  about  300  A.  To  realize  an  electric  field 

E  >  E  .  we  have  to  apply  voltages  of  about  1000  V.  This  means  for  a  sample 
min 

resistance  of  about  20a  required  electrical  power  of  200  kV.  We  use  a 
specially  designed  thyratron  tube  pulse  generator  which  supplies  5000  V  into 
25  (1.  An  impedance  transforming  network  is  used  in  connection  with  a  triax 
cable  operated  between  outer  and  inner  shielding  (impedance  ae  10  (2)  to  match 
the  low  Impedance  of  the  sample.  The  pulse  length  is  1  (is  and  the  repetlton 
rate  about  10  Hz.  for  heiger  rates  lasing  breaks  down  due  to  sample  heating. 

To  enhance  a  longitudinal  mode  two  approaches  have  been  used: 

a)  evaporation  of  metal  films  to  the  end  faces  of  the  sample 
so  that  they  can  serve  as  mirrors 

b)  attachment  of  external  mirrors  with  controlled  reflectivity 

In  the  first  case  we  evaporated  CaF2  as  an  Insulation  layer  and 
thereafter  1000  A  Au  for  the  10QJC  reflective  mirror  and  150  A  Ni-Cr  for  the 
semitransparent  one.  In  the  second  case  external  mirrors  are  used  in  a 
confocal  geometry:  One  mirror  is  plane  and  coated  with  1000  A  Au  on  one  side 
to  serve  as  100%  reflectivity  mirror.  The  second  mirror  is  piano  convex  and 
its  curved  surface  is  coated  with  different  thicknesses  of  Cr.  The  different 
mirror  reflectivities  were  measured  with  a  00a  pumped  FIR-gas  laser  and  a 
Golay  -  detector  to  be  between  80%  and  97%.  The  external  mirrors  are  attached 
to  the  sample  with  silicone  oil  and  fixed  with  nail  enamel. 


The  measurements  of  the  current  through  the  samples  show  the  expected 
behavior:  For  a  constant  applied  voltage  magnetic  field  the  current  increases 
sligthly  with  Increasing  magnetic  field  until  J;*1  =  1,  where  a  small  maximum 
occurs:  for  >  1  the  current  decreases  as  the  carriers  start  to  accumulate. 
This  behavolr  is  a  clear  proof  that  our  concentration  is  low  enough  to  enable 
streaming  motion. 

Fig.  3  shows  typical  results  of  the  spontaneous  far  infrared  emission, 
investigated  with  a  Ge(Ga) -detector.  The  magnetic  field  is  varied,  the 
different  traces  belong  to  different  applied  fields.  The  signal  increases,  as 
the  llgth  holes  start  to  accumulate  (f*  =  1).  The  maximum  is  between  f*1  s  1 
and  f*1  =  2.  where  the  heavy  holes  are  in  a  situation  between  streaming  motion 
and  accumulation.  For  higher  magnetic  fields  the  signal  decreases  as  an 
increasing  amount  of  heavy  holes  are  in  the  passive  region. 

Encouraged  from  the  streaming  behavior  evident  from  the  spontaneous 
emission  data,  a  Russian  group  /!/  and  Komlyama  et  al./8/  observed  stimulated 
emission  with  rather  long  samples  without  mirrors.  The  observed  spectra  showed 
a  clear  threshold  behavior  for  lasing.  However  the  spectra  were  quite  broad 
and  showed  many  unexplained  peaks.  The  mode  of  the  lasing  was  not  clear  at 
this  point.  At  this  stage  we  also  tried  to  obtain  lasing  in  samples  which 
showed  clear  evidence  for  streaming.  However  we  could  not  observe  lasing.  The 
reason  was  that  the  reflectivity  was  too  low. 

In  experiments  with  larger  samples  (50  x  5  x  5  mm)  we  were  able  to 
observe  stimulated  emission  even  without  mirrors.  Fig.  4  shows  the  magnetic 
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Fig.  3:  Magnetic  field  dependence  of  the  spontaneous  emission  signal  of 
a  p-Ce  sample  for  different  values  of  the  electric  field  recorded  with  a 
Ge'Ca.  -  detector. 

Fig.  4:  Magnetic  field  dependence  of  the  stimulated  emission  signal  from 
sample  1  for  different  electric  fields  detected  with  a  Ge:Ga  -  detector. 


field  dependence  of  the  detector  signal  for  different  electric  fields.  For 
certain  magnetic  fields  the  signal  Increases  for  two  orders  of  magnitude  (from 
0.5  mV  to  50  mV  for  Ge-detector  current  of  4*10‘5  A).  This  range  shifts  for 
higher  electric  fields  to  higher  magnetic  fields,  similar  to  the  maximum  of 
the  spontaneous  emission.  Emission  occurrs  for  electric  fields  E  between 
1000  V/cm  <  E  <  2000  V/cm  and  magnetic  fields  B  between  1  T  <  B  <  2  T. 
Expressed  by  the  the  normalized  field  f.  stimulated  emission  is  observed 

between  1.3  <  f*1  <  1.9.  These  data  agree  quite  well  with  the  results  of  other 

authors  /7/./8A  The  small  differences  in  the  onset  position  of  the  induced 
emission  can  be  explained  by  the  sllgthly  different  carrier  concentration  and 
an  other  crystal  orientation.  A  different  carrier  concentration  influences  the 
far  infrared  absorption  (free  carrier  absorption)  and  a  different  crystal 
orientation  parallel  to  the  magnetic  field  can  change  the  distribution 
function.  For  two  of  our  samples  the  orientation  of  the  long  axis  (this  is  the 
axis  parallel  to  the  magnetic  field)  is  (1,0.0).  For  a  third  sample  the 

orientation  of  the  long  axis  is  perpendicular  to  the  (1,0,0)  direction.  This 

sample  shows  no  lasing  activity.  For  maximum  output  power  in  the  literature 
the  orientation  parallel  to  the  magnetic  field  is  reported  to  be  the 
(1, l.l)-direction.  The  power  of  the  emitted  far  infrared  radiation  is 
estimated  by  a  comparision  with  other  FIR-sources  (002  pumped  gas  laser 
InSb-cyclotron  resonance  emitter)  to  be  about  300  mV.  This  is  lower  than 
Andronov  /10/  reported  for  his  samples  with  the  crystal  orientation  (1,1.1) 
parallel  to  the  magnetic  field. 

A  detailed  analysis  of  the  emission  spectra  was  performed  with  a 
GaAs-detector.  The  GaAs-detector  can  be  used  as  a  high  sensitivity 
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Fig.  5:  Spectrum  of  stimulated  emission  for  different  electric  and 
magnetic  fields. 

Sample  1:  E=  1.8  kV/cm.  B=  1.4  T 
Sample  3:  E=  1.9  kV/cm,  B=  1.7  T 

The  insert  schematically  shows  the  ray  path  in  the  Ce-crystal. 


pectrometer  when  tuned  with  the  magnetic  field.  Fig.  5  shows  emission  spectra 
(intensity  versus  frequency)  from  two  samples  for  different  values  of  the 
electric  and  magnetic  fields.  The  decreasing  background  level  is  due  to  the 
variation  of  the  detector  characteristic  with  magnetic  field.  The  spectra 
consist  of  10-20  lines  with  a  regular  mode  spacing.  Their  positions  remain 
fixed  for  changed  applied  fields.  The  measured  linewidth  is  about  0.4  cm"1.  We 
conclude  that  the  real  width  of  the  emission  is  of  the  same  order  as  the 
detector  linewidth  (0.25  cm"1).  By  varying  the  electric  and  magnetic  fields  it 
is  possible  to  obtain  lasing  between  75  pm  and  110  pm.  There  is  a  general 
trend  towards  emission  at  shorter  wavelength  for  both  higher  electric  and 
higher  magnetic  fields.  In  contrast  to  Andronov  et  al  /10/  and 
Komlyama  et  al./9/.  no  long  wavelength  lasing  (150-250  pm)  could  be  found  in 
our  samples,  which  is  probably  because  of  the  somewhat  higher  carrier 
concentration.  Andronov  et  al./7/  and  Komlyama  et  al./8/  suggested  that  the 
lasing  modes  are  due  to  total  reflections  at  the  sample  surfaces,  since 
without  mirrors  the  quality  factor  of  the  cavity  is  far  too  low  to  allow 
longitudinal  or  transverse  modes  (the  Ge  surface  has  36%  reflectivity  at 
normal  incidence).  According  to  Fig. 5  the  observed  mode  s pacings  are  1.74  cm"1 
and  1.20  cm"1  for  samples  1  and  3,  respectively,  and  1.40  cm"1  for  sample  2 
(not  shown).  Comparing  the  spaclngs  with  the  dimensions  of  the  samples,  it 
appears  to  be  very  likely  that  the  width,  a,  is  the  crucial  parameter 
determining  the  mode  spacing  (see  Table  I). 

The  mode  spaclngs  can  be  explained  by  considering  an  electromagnetic  wave 
in  the  germanium  crystal  (see  insert  in  Fig. 5)  propagating  with  components  of 
the  wavevector  in  the  direction  of  1  (z-direction)  and  a  (x-direction)/14/: 


transverse  resonance  condition  is  given  by  /15/  k^n  =  k  n  sin  0  =  Mir/a 
where  n  is  the  refractive  index,  M  is  an  integer,  and  0  is  the  angle  between  lc 
and  the  longitudinal  axis  of  the  crystal.  The  analogous  longitudinal  resonance 
condition  is  not  considered,  because  it  would  yield  a  very  narrow  mode 
structure  which  is  not  observed  in  the  experiment.  In  addition,  we  have  the 
geometrical  condition  tan  0  =  Na/1  (N  integer),  which  ensures  that  the  path 
of  the  ray  is  closed  within  the  sample.  This  leads  to  a  mode  spacing  (in 
wavenumbers)  Av  =  (2a  n  sin  0)”1.  0  must  satisfy  the  condition 

14.5°  <  0  <  75. 5°.  because  14.5°  is  the  angle  of  total  reflection  for 
germanium.  One  can  assume  that  modes  with  the  smallest  possible  0  have  the 
lowest  losses,  because  they  undergo  the  fewest  reflections.  Thus,  for  each 
sample  the  mode  spacing  is  determined  first  by  choosing  the  angle.  0,  with  the 
smallest  possible  N.  and  then  this  angle  is  used  in  the  expression  for  Av. 
This  procedure  gives  a  mode  spacing.  Av,  of  1.77  cm-1  (N=6).  1.41  cm-1  (rWS). 
and  1.14  cm"1  (N=4).  for  sample  1.2  and  3,  respectively,  in  good  agreement 
with  the  measured  values.  This  interpretation  is  confirmed  by  the  work  of 
Andronov  et  al./10/,  who  reported  emission  of  a  single  line  by  using  a  very 
thin  (a  =  0.7  nan)  sample,  where  the  mode  spacing  is  as  wide  as  the  gain 
region. 

After  attaching  external  mirrors  to  sample  3  (see  insert  in  Fig. 6),  the 
region  of  electric  and  magnetic  fields  where  lasing  occurs  stays  roughly  the 
same.  The  detector  signal  is  also  of  the  same  order  of  magnitude  with  and 
without  mirrors.  The  observed  spectra,  however,  exhibit  a  dramatic  change. 
Fig. 6  shows  spectra  of  sample  3  for  different  electric  and  magnetic  fields, 
recorded  with  the  GaAs-detector .  A  considerably  reduced  number  of  lines 
appears  with  wider  mode  spacing  and  higher  amplitude,  accompanied  by  some 
other  less  dominant  structure.  Varying  the  electric  and  magnetic  field  allows 
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B  =1.6  T 


Ei&i_6:  Emission  spectra  from  sample  3  after  attaching  external  mirrors 
for  different  electric  and  magnetic  fields.  The  Insert  shows  the  cavity 
configuration. 
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to  switch  on  and  off  lines  on  the  long-  or  short-  wavelength  side  of  the 
spec t rum  due  to  a  change  in  the  gain  spectrum.  In  this  way  it  was  possible  to 
obtain  operation  of  a  single  line  (lower  part  of  Fig. 6)  /16/.  This  was  the 
first  report  on  forcing  a  multimode  p-Ge  laser  to  oscillate  on  a  single  line 
by  attaching  external  mirrors.  The  emission  wavelength  is  87.3  pm  (114.5 
cm-1),  the  adjacent  line  at  83.3  pm  (120.0  cm-1)  is  just  at  threshold.  As  seen 
from  Fig. 6  the  single  line  contains  approximately  the  same  total  intensity  as 
the  three  lines  together.  However  one  cannot  exclude  the  coexistence  of 
several  cavity  inodes  within  this  line,  because  their  spacing  would  be  below 
the  resolution  of  the  detector.  The  spectra  observed  by  Andronov  et  al./10/ 
from  samples  with  external  mirrors  showed  a  very  irregular  and  complicated 
structure.  Komiyama  et  al./17/  reported  some  changes  in  the  spectrum  and  were 
able  to  reduce  the  number  of  oscillating  lines  by  sandwiching  a  thin  layer  of 
poly-4-methylpent-l-ene  (manufactured  under  the  trade  name  TPX  by  Imperial 
Chemical  Industries)  between  the  sample  and  one  mirror. 

The  mode  spacing  in  Fig. 6  is  equal  to  Au  =  5.5  cm-1,  which  is  too  large 
for  fundamental  longitudinal  or  transverse  cavity  modes.  A  possible 
explanation  is  the  formation  of  an  intracavity  interference  filter  by  the  two 
Ge-plates:  Since  they  are  not  subject  to  the  electric  field  like  the  sample, 
they  are  expected  to  have  a  slightly  different  index  of  refraction  and  also 
the  thin  layer  between  sample  and  mirror  Induces  a  change  of  the  refractive 
index,  which  provides  the  possibility  of  Fabry-Perot  like  reflections.  A 
thickness  of  dt=  0.9  nan  and  d2=  1.1  mm  (the  slight  curvature  of  the  convex 
plate  is  neglected)  would  give  a  spacing  of  only  1.4  cm-1  and  1.1  cm-1, 
respectively.  If.  however,  the  combined  resonance  condition  for  both  plates  is 
considered,  the  spacing  is  calculated  to  5.5  cm-1.  This  result  is  based  upon  a 
ratio  of  d2/dt  equal  to  5/4.  Considering  an  uncertainty  of  the  measured 


thicknesses  of  0.05  nan  and  a  low  finesse  of  the  Fabry-Perots.  one  gets 
spacings  between  4  cm'1  and  6  cm-1.  A  definitive  explanation  of  the  present 
spectra  would  require  knowledge  about  the  refractive  index  of  the  lasing 
sample  as  well  as  the  influence  of  the  thin  layer  between  sample  and  mirrors. 
The  Investigation  of  different  mirror  configurations,  especially  with 
different  thicknesses,  appears  necessary. 

We  also  checked  the  emission  when  the  plane  mirror  had  a  reflectivity  of 
only  90%.  In  this  case  no  lasing  was  achieved,  which  confirms  the 
interpretation  as  waveguide-like  modes  in  the  mirrorless  configuration:  These 
modes  are  destroyed  by  attaching  mirrors,  which  prevents  total  reflection  on 
both  end  faces.  Considering  only  the  mirror  output  and  no  other  losses  the 
gain  is  calculated  to  be  about  5*10~3  cm'1. 


Table  I :  Sample  dimensions  and  observed  mode  spacings:  1  is  the  sample  length 
((lll)-direction).  d  is  the  distance  between  the  electrical  contacts,  a  is  the 
width  and  At>  is  the  measured  mode  spacing. 


Sample 

1 

d 

a  (<m») 

Av  (i 

1 

46.9 

3.7 

2.4 

1.74 

2 

43.1 

4.2 

3.15 

1.40 

3 

48.4 

5.2 

3.75 

1.20 
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In  this  chapter  we  want  to  discuss  the  most  important  ideas  from  a 
theoretical  point  of  view  for  the  build-up  of  an  Inversion  in  the  crossed 
field  configuration  and  review  calculations  of  the  distribution  function.  The 
main  part  will  be  devoted  to  a  new  approach  for  the  calculation  of  the 
distribution  function  in  crossed  fields  which  explicitly  considers  the  LL 
structure  of  the  energy  spectrum.  Only  one  type  of  carriers  population 
Inversion  is  considered  within  this  carrier  system. 

The  existence  of  an  Inverted  distribution  in  crossed  fields  was  first 

predicted  by  Kaeda  and  Kurosawa  in  1972  /l/.  They  performed  Monte  Carlo 

calculations  and  also  tried  to  understand  the  situation  by  qualitative 

arguments:  They  found  out  the  importance  of  the  optical  phonon  scattering  for 

the  generation  of  a  non  thermal  distribution.  Their  model  was  refined  by  other 

authors  /2/,/3/  and  became  famous  under  the  name  "streaming  motion  and 

accumulation".  This  concept  is  based  on  the  classical  motion  of  carriers  in 

crossed  fields  considering  a  sufficiently  strong  interaction  with  optical 

phonons.  In  crossed  fields,  the  carriers  move  on  circles  in  velocity  space 

with  central  point  v^  =  -(E/B) .  The  energy  of  optical  phonons  hu^  corresponds 

to  another  circle  v  =  v  =  V2hu  /m*  centered  at  the  origin  v  =  0.  When  the 

op  op 

trajectory  of  a  carrier  crosses  the  circle  v  =  v  ,  the  carrier  can  emit  an 
optical  phonon  and  is  scattered  back  near  the  origin. 


In  the  Halt  of  large  electric  fields.  E/B  >>  v  .  all  trajectories  are  open. 

op 

l.e.  they  cross  the  circle  v  =  v  .In  this  case  the  motion  is  only  slightly 

op 

Influenced  by  the  magnetic  field  and  is  called  "streaming  motion"  because  of 

the  successive  acceleration  and  phonon  emission.  In  the  other  limiting  case 

E/B  <<  v  most  of  the  trajectories  are  closed  and  the  carriers  are  not 
op 

affected  by  optical  phonon  scattering.  In  the  intermediate  range  a  fraction  of 
the  carriers  performs  streaming  motion,  while  another  part  moves  on  closed 
trajectories.  The  streaming  carriers  have  a  finite  probability  of  being 
scattered  into  the  area  of  closed  trajectories,  the  so  called  passive  region, 
where  they  can  accumulate.  They  can  only  be  removed  from  there  by  acoustic 
phonon  or  impurity  scattering,  which  are  both  much  slower  than  the  optical 
mode  scattering.  As  a  consequence,  the  carriers  get  accumulated  (or  bunched) 
in  a  limited  area  of  momentum  space. 

Calculations  showed  that  for  f  <  2  a  torus-type  distribution  function  can 
arise  in  this  accumulation  area,  which  corresponds  to  a  population  inversion 
between  Landau  levels.  Analytical  approaches  based  on  the  Boltzmann  equation 
have  been  performed  for  this  situation  /4/  as  well  as  Monte  Carlo 
calculations  /5/.  mostly  applied  to  n-InSb,  n-GaAs  and  the  heavy  holes  of 
p-Ge.  Very  recently.  Abou  El  Ela  and  Ridley  /6/  made  a  Monte  Carlo  simulation 
for  n-GaAs.  including  electron-electron  scattering  and  Intervalley  transfer. 

Experimentally,  the  accumulation  of  carriers  was  proved  in  different  materials 
/7/.  but  it  was  not  clear,  if  a  population  inversion  between  Landau  levels 
could  be  achieved.  A  short  time  ago.  strong  far  infrared  emission  signals  have 
been  observed  in  p-Germanium  /&/ ,  which  were  ascribed  to  a  LL  inversion  of 
light  holes. 


Obviously,  the  most  important  application  of  an  inverted  distribution  function 
is  the  generation  and  amplification  of  coherent  radiation,  in  this  case  caused 
by  transitions  between  Landau  levels.  According  to  this,  the  main  drawback  of 
the  theoretical  work  mentioned  above  is  the  fact  that  it  does  not  account  for 
the  discrete  structure  of  the  energy  spectrum,  i.e.  the  Landau  levels. 

All  these  studies  were  performed  in  the  limit  of  nonquantizing  magnetic 
fields,  that  means  the  use  of  a  classical  distribution  function  smeared  out 
over  the  I .«nA»n  levels  and  the  use  of  the  zero-magnetic  field  scattering 
rates.  These  classical  approaches  are  not  suited  to  yield  information  about 
the  population  of  different  Landau  levels. 

In  this  work  we  treat  the  problem  of  carriers  in  crossed  fields  in  a  quantum 
mechanical  frame,  i.e.  the  Landau  levels  are  automatically  incorporated  in  the 
calculation. 

We  consider  the  Hamiltonian  for  a  carrier  with  charge  e  and  effective  mass  m* 
in  crossed  electric  and  magnetic  fields  Interacting  with  phonons  /9/./10/: 

h  -  -*?•?  *  %  ♦  Hel_ph  (i) 

Here  is  the  free  phonon  term: 

“ph'pq^V141  (2) 


where  the  b  and  b  are  the  phonon  creation  and  destruction  operators. 

q  q 


respectively,  and  Hel_ph  Is  the  electron-  (hole-)  phonon  coupling: 

H.i-ph  • 1  5  c<?>  <bq  'i5?  -  <  «'l5?)  <3> 

C(q)  depends  on  the  specific  nature  of  the  electron-phonon  coupling. 
Eigenvalues  and  eigenfunctions  of  the  electronic  part  of  the  Hamiltonian  (E  in 
z-direction.  £  in  x-direction)  are  given  by 

Eu  =  (n+H)  huc  +  f»2k2/2m*  -  eEX  +  K  m*(E/B)2  (4) 

where  the  index  v  summarizes  the  quantum  numbers  n.  k  .  and  k  ,  u  =  |e|B/ro*. 

y  x  c 

and  X  a  (hky/eB  +  m**E/eB2). 

*.(f)  *  e' V  A1  V5?)  <5> 

where  1  =*Ai/ |e |B  and  Is  the  n-th  harmonic  oscillator  wave  function. 

The  basic  quantum  transport  equation  is  the  equation  of  motion  for  the  density 
operator,  the  Liouvllle  -  von  Neumann  equation 

ih  =  [H.  p]  (6) 

Here,  p  is  the  density  operator  of  the  total  electron  -  phonon  system.  The 
electron  distribution  function  is  given  by  the  diagonal  elements  of  the 
reduced  density  operator,  f  ,  which  is  obtained  by  taking  the  trace  over  the 
phonon  states: 


Assuming  no  correlation  between  electron  and  phonons  in  the  diagonal  elements 
of  p  and  a  weak  electron  -  phonon  coupling  the  Pauli  master  equation  can  be 
derived  for  the  diagonal  elements  of  the  reduced  density  operator  (neglecting 
the  Pauli  exclusion  principle)/10/,/ll/,/12/ 


370  «  Z.  f  .  V  .  -f  W  . 
dt  v  v  v  v  v  vv 


The  W .  are  the  scattering  probabilities  from  state  i>  to  state  i>*  according 
to  Fermi’s  golden  rule  : 


1<" IC(^)  .lqr|»->|*  [N _WN +l)«(E1)-E1).-hU  )] 


q  '  q 


is  the  Planck  distribution  for  the  phonons,  the  matrix  elements  are  given 
by  (see  e.g.  Ref.  13) 


l<»|.‘qrl»,>l2  ■  l-WtvV1’  V.i 


VW 


with  |Jm.{f)|a  =  Siy  f"'"n  e'f  [L”''n(f)]a  Wn'), 


and  f  as  l2(q£  +  q^)/2.  The  L“  “  are  associated  Laguerre  polynomials  /14/.  It 
is  worth  noting  that  the  electric  field  only  appears  in  the  energy  conserving 
6-function  and  not  in  the  squared  matrix  elements. 


We  now  proceed  to  get  a  steady  state  solution  of  the  master  equation  (8)  which 
then  reads 


0  =  2.  f  .  W  .  -  f 

V  V  V  V  V 


vv 


(11) 


In  the  following  we  want  to  consider  the  situation  where  the  current  through 

the  sample  is  constant.  Therefore  we  assume  a  spatially  homogeneous 

distribution  function.  Thus,  the  f  have  to  be  independent  of  X.  and  therefore 

of  ky  /9/./19/.  This  enables  us  to  perform  the  summation  over  directly  and 

we  are  left  with  a  system  of  N  coupled  integral  equations  (for  N  LL)  in  the 

variable  k  : 

z 


o  =  j.  j-dk;  f„.(k;)  »n.n(k2.ig  -  fn(k2)  Vfk^)  (12) 


The  rates  W  .  (k  .k*)  are  now  also  k  -integrated.  The  distribution  functions 
nn  z  t.  y 

for  the  different  LL‘s  can  be  expressed  as 


W  =  pr^tr,  ■ 


(13) 


where  W  . (k  )  =  fdk'  W  ,(k  ,k'). 

nn  v  zy  z  nn  v  z  z' 


The  proper  normalization  condition  is  2j*dk  f  (k  )  =  1. 

n  z  n  z 


Our  first  attempt  to  solve  this  equation  was  done  by  dropping  the 
k^-dependence .  Then  Eq.  (13)  becomes  a  system  of  algebraic  equations. 
Physically,  this  approximation  contains  an  overestlmatlon  of  transitions  due 
to  a  change  in  the  spatial  variable  X.  and  an  underestimation  of  transitions 


i 


due  to  a  change  in  the  quantum  number  k  .  This  approach  is  shown  in  Ref.  /IS/ 

z 

for  3  Landau  levels,  including  polar  optical  and  acoustic  deformation 
potential  phonon  scattering,  which  is  adequate  for  electrons  in  GaAs  or  InSb. 
For  the  first  time,  the  possibility  of  a  population  inversion  between  Landau 
levels  is  demonstrated  theoretically  in  a  quantum  mechanical  frame. 

Recently,  we  have  extended  this  approach  by  taking  into  account  the 

k  -dependence  and  solving  the  coupled  integral  equations  (13)  for  3  Landau 
z 

levels.  This  work  is  contained  in  Ref.  /16/,  for  optical  and  acoustic 
deformation  potential  scattering  (adequate  e.g.  for  the  light  holes  in 
Germanium) . 

A  population  inversion  is  predicted  also  for  this  case.  In  addition,  the 
one-dimensional  distribution  functions  in  the  different  Landau  levels  are 
discussed,  and  special  attention  is  paid  to  a  comparison  of  the  quantum 
mechanical  picture  with  the  classical  picture  of  streaming  motion. 

The  most  important  feature  for  a  further  development  of  our  model  is  the 
inclusion  of  the  influence  of  the  heavy  hole  band  to  the  population  of  the 
light  hole  band  in  p-Ge.  As  the  redistribution  of  carriers  from  the  heavy  hole 
bond  can  enhance  the  inversion  in  the  light  hole  band. 
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Part  II 


A)  FWFRCT  KEIAXATTGH 


The  mechanism  of  the  energy  relaxation  and  the  form  of  the  distribution 
function  are  of  fundamental  interest  for  the  understanding  of  electric  field 
effects  in  two-dimensional  electron  systems.  At  low  temperatures  where  lattice 
scattering  is  weak,  hot  electron  phenomena  can  be  produced  by  fields  as  low  as 
V/cm  in  Si  as  well  as  GaAs  inversion  layers. 


In  the  present  report,  two  different  methods  are  used  to  determine  the  elctron 
temperature  and  the  energy  relaxation  time: 

a)  The  electron  gas  is  heated  up  by  electric  field  pulses  and  the  emitted 

broadband  far  infrared  radiation  is  analyzed  in  terms  of  a  black  body 

source  with  a  temperature  T  .  The  derived  electron  temperatures  are 

© 

Independent  of  sample  parameters  when  plotted  against  the  input  power  per 

electron  (P@  =  epE^).  A  nearly  linear  dependence  of  Te  over  Pfi  is  found  for 

electron  temperatures  between  10  K  and  40  K.  At  higher  T  -values  a  weaker 

© 

increase  is  found  which  tends  to  saturate  at  T@  ~  100  K.  From  an  energy 
balance  consideration  energy  relaxation  times  are  obtained  which  decrease 
linearly  with  increasing  electron  density  /1/./2 /. 


As  a  dominant  energy  relaxation  mechanism  acoustic  deformation  potential 
scattering  is  found  for  electron  temperatures  up  to  40  K.  At  higher 
electron  temperatures  optical  phonon  emission  dominates  the  energy 
relaxation.  However,  the  obtained  relaxation  times  are  longer  than  the 
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theoretically  expected  times  for  a  bulk  like  phonon  process.  This  is  an 
indication  for  a  possible  influence  of  hot  phonon  effects. 

b)  A  direct  way  to  obtain  the  energy  relaxation  rate  is  the  measurement  of  the 
incoherent  saturation  of  the  cyclotron  resonance  transition.  This  is  done 
under  the  assumption  that  a  magnetic  field  does  not  change  the  relaxation 
rates  significantly.  The  relaxation  times  derived  are  comparable  to  the 
values  found  in  the  case  without  magnetic  field  /3/.  As  the  CR-saturation 
is  a  direct  technique,  these  experiments  can  be  used  as  a  prove  for  the 
validity  of  the  black-body  emission  analysis. 

Both  techniques  give  energy  relaxation  times  between  1  ns  and  0.2  ns  depending 
on  the  2-dimensional  carrier  density.  These  times  agree  with  values  obtained 
in  the  bulk  if  compared  with  the  equivalent  3-dlmensional  density. 

B)  ANALYSIS  OF  THE  CARRIER  DISTRIBUTION  FUNCTION  THROUGH  SMITH-PUHmj.  FFFBCT 
IN  GaAs/GaAIAs  HETEROSTKUCTUKES 


The  two-dimensional  electrons  in  selectively  doped  GaAs/GaAIAs  heterojunctions 
exhibit  extremely  high  mobilities  at  low  temperatures  /4-6/.  The  electrons  are 
easily  heated  up  in  an  external  electric  field  leading  to  a  rapid  rise  in 
electron  temperature  /1,7,8/  and  a  reduction  of  electron  mobility  /9/.  No 
detailed  analysis  of  the  form  of  the  distribution  function  has  been  performed 
for  extremely  high  mobilities.  In  this  case  normal  Fermi  distributions  cannot 
be  used  since  drift  energies  comparable  with  the  Fermi  energy  can  be  achieved 


AJUXId'U'lLlltj 


already  at  low  electric  fields. 


In  the  present  paper  a  direct  analysis  of  the  electron  distribution  function 
for  extremely  high  mobility  samples  as  a  function  of  the  external  electric 
field  is  presented.  The  Smith-Purcell  effect  /10/  is  measured  in  a 
semiconductor  for  the  first  time:  A  periodic  potential  of  period  L  Induces  an 
energy  loss  of  a  drifting  carrier  distribution  via  photon  emission.  The  angle 
between  the  drift  and  the  grating  momentum  is  varied  by  using  different  sample 
geometries.  The  geometry-dependent  emission  spectra  are  momitored  through  a 
magnet ic-f ield- tuned  InSb-detector  /ll/. 


Through  the  selective  momentum  transfer  (±  q  =  2r/L)  the  spectral  analysis  of 
the  energy  loss  yields  detailed  information  on  the  drifted  electron 
distribution  function.  For  a  momentum  transfer  much  smaller  than  the  electron 
momentum,  the  frequency  u  of  the  emitted  light  is  directly  proportional  to  the 
particle  wavevector  k: 


u  as  hkq/m*  =  v^q 


(q  «  k). 


m*»  is  the  electron  mass.  Thus  the  emitted  spectra  reflect  directly  the 
distribution  of  the  carrier  velocities  v^.  In  the  following  analysis  the 
electron  distribution  function  will  be  described  as  a  shifted  Fermi  function: 


f (E)  .  1/(1  +  exp(h2(k  -  k D)2/2m»  -  p)/Te), 


introducing  the  drift  momentum  k^  and  the  drift-dependent  electron  temperature 
T  as  phenomenological  parameters.  Our  model  assumes  that  the  distribution 
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function  is  stationary,  depending  only  on  the  electric  field  and  the  sample 
mobility.  The  FIR  coupling  represents  only  a  neglibible  cooling  process. 

Theory  of  the  emission  process 

FIR  emission  from  a  free  electron  system  requires  a  second  scattering  process 
for  momentum  conservation,  which  can  be  either  impurity  or  phonon  scattering 
/12/.  An  electromagnetic  calculation  of  the  FIR  emission  from  heated  carriers 
through  a  grating  was  recently  performed  in  ref.  1.  For  the  electric  field 
driven  electrons,  the  lateral  potential  modulation  UQCos(qx)  represents  a 
coherent  scattering  source  (the  lattice  vector  has  been  taken  in  the  x 
direction).  Formally,  it  can  be  described  as  one  "extended"  impurity, 
scattering  the  electrons  by  ±q.  The  emission  spectrum  of  a  drifted  electron 
distribution  depends  on  the  angle  between  the  lateral  structures  and  the 
electric  field  and  thus  the  current  direction. 


fig.  1i 

(a)  Scattaring  diagrams  used  for 
the  theoretical  calculation. 

(b)  Calculated  emission  Intensity 
versus  frequency  from  drifted  elec¬ 
tron  distributions  with  a  tempera¬ 
ture  of  35 K,  comparing  jxq  (full 
curve)  with  j//q  (broken  curve). 
The  spectra  are  weighted  by  u  to 
account  roughly  for  the  detector 
sensitivity. 
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We  calculate  the  emission  rate  in  the  lowest  order,  including  only  the 
scattering  processes  shown  in  Fig.  (la).  This  diagrammatic  approximation  can 
be  controlled  analytically  also  for  "badly"  shaped  electron  distributions. 
Static  screening  is  included  in  the  impurity  potential.  Evaluating  the 


diagrams  Fig.  (la)  and  performing  the  sum  over  the  photons  the  energy  loss 
rate  is  obtained  as: 


dW  _  e1  Ujj.q* 
dt  4ir,m*c,fi* 


1 


/dkyf  (k*,ky)  (1  -  f  <k*  ±q.  ky)  ) 


ki 


*<? 


iiim* 

fiq 


(3) 


The  sum  over  (1)  describes  the  two  possible  momentum  transfers,  and  the 

^-integration  has  already  been  performed  by  exploiting  the  delta  function  for 

energy  conservation.  The  influence  of  the  electron  distribution  f(k  ,k  )  is 

x  y 

contained  in  the  ky  sum,  which  is  evaluated  numerically. 

Fig.  (lb)  shows  the  calculated  energy-dependent  emission  intensity  for  the 
parameters  of  the  sample  described  below,  for  an  electron  temperature  of  35  K. 
The  curves  are  weighted  by  u  in  order  to  account  roughly  for  the  detector 
sensitivity.  The  full  curve  describes  a  carrier  distribution  which  is  heated 
by  an  electric  field  normal  to  the  grating  wavevector  q.  In  this  direction  the 
result  is  independent  of  the  drift  velocity  if  a  shifted  Fermi  distribution  is 
used.  In  contrast,  the  broken  line  shows  the  result  for  the  electric  field 
parallel  to  q,  with  v^  =  0.7  Vp  :  A  substantial  change  of  the  emission 
characteristics  is  obtained.  The  spectrum  splits  into  two  broad  structures. 
The  shifts  to  higher  and  lower  energies  demonstrate  the  transition  from  normal 
free  carrier  emission  to  the  Smith-Purcell  regime  (the  low  energy  side 
disappears  for  still  higher  drift  velocities).  While  the  lower  part  of  the 
spectrum  is  out  of  the  range  of  our  detector,  the  change  on  the  high  energy 
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side  should  be  observable  when  changing  the  orientation  of  the  current  with 
respect  to  the  grating. 

Experiment 

We  have  plasma  etched  a  periodic  grating  with  lattice  constant  L  into  the  top 

layer  of  a  GaAs/GaAlAs  heterostructure  as  shown  in  Fig.  (2a).  The  height  of 

the  modulation  was  of  the  order  of  200  A  and  the  electron  density  is  2.8  x 
11  2 

10  /cm  .  The  modulation  Induces  a  periodic  potential  acting  on  the  20 

electrons,  which  is  evident  from  Shubnikov-de  Haas  (SdH)  oscillations.  At  4.2 

6  2 

K  the  samples  had  a  zero  field  mobility  of  2.0  x  10  cm  /Vs  before  plasma 

etching,  and  very  narrow  minima  in  a  .  After  etching  the  mobility  was  only 

6  2 

reduced  to  1.5  x  10  cm  /Vs,  however  the  minima  became  much  broader.  We 

interpret  this  broadening  as  the  appearance  of  a  sinusoidal  potential,  leading 

to  a  local  variation  of  the  electron  density.  From  the  width  increase  of  the 

SdH  oscillation  we  estimate  the  size  of  the  periodic  potential  to  be  of  the 

* 

order  of  1  meV. 
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Fig.  3 
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Fig.  2; 

Schema t leal  drawing  of  the  GaAs-GaAlAs  structures  with  periodic  grating  used  in 
the  experiments  (a).  Sample  1  with  the  contacts  parallel  to  the  grating  (b)  and 
sample  2  with  contacts  both  parallel  and  normal  to  the  grating  wavevector  (c) . 

EAii-li 

Mobility  and  drift  velocity  of  the  electrons  as  a  function  of  the  electric  field. 
Open  symbols:  Before  plasma  etching.  Full  symbols:  After  plasma  etching.  Due  to 
the  high  sample  mobility  the  drift  velocity  saturates  already  at  low  electric 
fields. 
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Tvo  different  sample  geometries  were  realized:  Sample  1  had  the  grating 
grooves  parallel  to  the  contacts  (Pig.  2b).  and  as  a  result  q//j.  In  sample  2 
both  geometries  (q//J  and  qij)  are  realized.  The  drift  velocity  and  mobility 
versus  the  external  electric  field  as  derived  from  direct  conductivity 
measurements  of  the  small  sample  structure  (Fig.  2b)  at  4.2  K  before  and  after 
the  plasma  etching  are  shown  in  Fig.  (3).  It  is  evident  that  the  saturation 
drift  velocity  is  already  obtained  for  electric  fields  in  the  order  of  100 
V/cm.  The  observed  dependences  of  and  p  are  in  good  agreement  with  the 
theoretical  calculation  of  Lei  et  al.  /13/  based  on  a  non-Boltzmann  balance 
equation  approach  and  using  a  drifted  Fermi/Dirac  distribution  with 
temperature  T  . 

The  analysis  of  the  emitted  radiation  in  the  presence  of  a  periodic  potential 
allows  us  to  answer  the  question  whether  the  drift  velocities  are  realistic 
and  what  is  the  average  temperature  of  .  the  carrier  distribution.  Emission 
spectra  for  a  plasma  etched  sample  1  (Fig.  2b)  for  several  electric  fields  are 
shown  in  Fig.  (4).  For  the  lowest  electric  field  the  detector  response  is  only 
due  to  broad  band  thermal  emission  from  the  heated  electron  distribution  /l/. 
Above  an  electric  field  of  30  V/cm  the  emission  spectra  start  to  show  a 
structure  around  a  magnetic  field  of  0.3  T.  which  develops  to  a  clear  but 
broad  peak  at  100  and  150  V/cm.  The  spectra  show  a  clear  shift  between  50  and 
100  V/cm.  indicated  by  the  arrow  which  marks  the  decay  of  the  signal. 


For  samples  without  a  grating  the  spectra  show  a  similar  fora  as  the  curve  for 
10  V/cm  for  all  electric  fields,  but  the  broad  peak  around  0.3  T  is  missing. 
Consequently  the  additional  structure  can  be  directly  attributed  to 


Smith-Puree 11  type  emission. 


To  deaonstrate  that  the  observed  emission  is  due  to  the  grating  Induced 
potential  we  prepared  a  sample  where  the  grating  wavevector  q  is  oriented  both 
parallel  and  normal  to  the  current  direction.  In  Fig.  (5)  a  comparison  of  the 
emission  for  q//J  (dashed  curve)  and  qlj  (full  curve)  is  shown  for  an  electric 
field  of  100  S/cm.  In  both  orientations  radiation  due  to  the  Smith-Purcell 
effect  is  observed.  It  is  evident  that  the  spectrum  for  q//J  is  shifted  to 
higher  energies.  In  the  qlj  case  the  scattering  depends  essentially  on  the 
carrier  temperature,  independent  of  the  drift  velocity.  The  emission  comes 
mainly  from  carriers  moving  with  the  Fermi  velocity  perpendicular  to  the 
current  direction. 

riq.  5 

Fig.  4 
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Emission  spectra  from  drifted  electron  distributions  with  j//q,  for  different 
electric  field  strengths.  The  broad  structure  appearing  at  high  fields  is  a  con¬ 
sequence  of  the  grating-induced  emission. 


Fig.  5i 

Emission  spectra  from  drifted  electron  distributions,  comparing  geometries  with 
j//q  and  J±q.  *  magnetic  field  normal  to  the  interface  destroys  the  selective 
emission . 


From  this  experiment  it  is  directly  evident  that  the  shift  of  the  spectra  is 
due  to  the  drift  of  the  carriers.  The  spectra  thus  show  the  behaviour 
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predicted  by  the  theory  (Fig.  lb).  For  the  given  situation  the  shift  can  be 
converted  into  the  drift  velocities,  although  we  have  to  keep  in  mind  that  the 
detector  has  a  rather  limited  resolution  of  1  meV  and  that  the  emitted 
spectrum  is  quite  broad. 

As  an  additional  proof  that  the  emission  is  due  to  the  high  carrier  velocity 
we  reduced  the  carrier  velocity  with  an  applied  magnetic  field  perpendicular 
to  the  2D-plane.  The  selective  emission  and  the  difference  in  the  spectra 
between  J//q  and  jiq  disappeared  as  shown  in  Fig.  (5)  for  a  magnetic  field  of 
1.5  T. 

7 

As  a  rough  estimate  we  obtain  a  drift  velocity  of  1.5  x  10  cm/s  for  100  V/cm. 
This  drift  velocity  is  somewhat  smaller  than  that  derived  from  current 
'measurements.  In -addition  we  can  estimate  the  carrier  temperature:  It  has  to 
be  in  the  order  of  »  50  K  from  a  fit  of  the  decaying  part  of  the  spectra  which 
are  marked  by  arrows.  <Xir  calculations  show  that  for  considerably  higher 
temperatures  (T  >  80  K)  the  spectra  become  extremely  broadened  and  it  would 
not  be  possible  to  observe  well-defined  structures  with  the  given  resolution 
of  the  InSb  detector. 

Cl  SUMMARY. 

A  significant  difference  in  electron  temperatures  is  found  from  the  black-body 

emission  analysis  and  the  Smlth-Purcell  effect.  In  both  cases  samples  with 

6  2 

mobilities  higher  1  x  10  cm  /Vs  at  4.2  K  were  used.  In  the  Smlth-Purcell  case 
a  rather  significant  drift  velocity  is  found  which  results  in  a  shifted 
Fermi-Dlrac  distribution.  In  the  hot  electron  emission  experiment  a 


t hemodynamic  equilibrium  distribution  was  assumed,  which  is  not  a  good 
description  of  the  real  situation.  In  the  case  of  the  emission  experiments  we 
ran  therefore  conclude  that  the  electron  temperature  model  gives  problematic 
electron  tesiperatures  (T  »  100K)  if  a  significant  drift  energy  is  present. 

The  obtained  drifted  Fermi-Dirac  distribution  and  the  rather  low  electron 
temperatures  derived  from  the  Smith-Purcell  effect  are  well  explained  with  the 
theory  of  Lei  et  al.  /13/. 
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We  observe  far  infrared  broadband  emission  from  the  hot  two-dimensional  electron  plasma  in 
GaAs/AIGaAs  single  heterostructures  grown  by  molecular  beam  epitaxy.  The  radiation  is 
analyzed  in  two  different  frequency  regimes  (around  35  and  100  cm  _ ').  From  the  relative 
dependence  of  the  intensities  on  the  applied  longitudinal  electric  held  the  hot-electron 
temperatures  are  determined.  From  the  absolute  emission  intensities  values  of  the  free-carrier 
absorption  coefficient  (down  to  10 ~ 5  in  ultrahigh  mobility  samples)  are  measured. 


The  topic  of  carrier  heating  in  the  lateral  transport  of 
two-dimensional  (2D)  systems  along  semiconductor  inter¬ 
faces  and  heterojunctions  has  attracted  much  interest  in  the 
last  years,  especially  due  to  the  applications  of  these  systems 
for  fast  field-effect  transistors'  and,  recently,  for  hot  carrier 
and  real  space  transfer  devices. 2  There  are  several  theoretical 
papers3  on  the  problem  of  carrier  heating  in  GaAs/AIGaAs 
heterojunctions.  Two  main  experimental  approaches,  name¬ 
ly,  hot-carrier  photoluminescence4  and  magnetotransport, 3 
have  been  made  to  measure  the  carrier  distribution  function 
or  the  carrier  temperature,  respectively. 

In  this  letter  we  report  on  far  infrared  (FIR)  emission 
experiments  that  were  performed  by  passing  a  lateral  current 
through  the  2D  electron  system  of  GaAs/AIGaAs  single 
heterojunctions.  This  heats  up  the  carriers  to  a  carrier  tem¬ 
perature  T.  It  has  been  shown  in  the  experiments  of  Shah  et 
al.*  and  in  femtosecond  studies  of  the  intraband  relaxation 
by  Erskine  et  a/.b  that  the  carriers  thermalize  very  fast  to  a 
thermal  distribution  characterized  by  the  carrier  tempera¬ 
ture  T,  via  carrier-carrier  interaction  within  less  than  100  fs.6 
Since  this  is  by  orders  of  magnitude  less  than  the  energy 
relaxation  time  (10~  lo-10~<>  s)  the  carrier  system  itself  rep¬ 
resents  a  system  in  thermal  equilibrium. 

The  FIR  emission  therefore  can  be  described  by  ther¬ 
modynamical  considerations.  In  equilibrium  the  spectral 
emission  intensity  /  |<u,r)of  electromagnetic  radiation  from  a 
system  with  temperature  T  is  given  by 

/(fti.T)  =*  Im(co,T)A  (<u),  (1) 

where  I BB  is  the  spectral  emission  intensity  of  a  black  body 
and  A  is  the  spectral  absorptivity  of  the  system,  in  a  2D 
system  for  normal  incidence 


A  M  = 


4  Re  F 

1^+1  +F  |2’ 


(2) 


with  F  =  o\b>)/€<p  and  e  being  the  dielectric  constant  of  the 
substrate.  A  Drude  form  of  the  dynamical  conductivity  aita) 
in  n-Si  inversion  layers  has  been  confirmed  by  microwave 
and  FIR  transmission  experiments  by  Allen  et  al .7  up  to 
frequencies  of  40  cm  ~ 1 .  This  allowed  the  determination  of  T 
from  the  absolute  FIR  emission  intensity  at  one  frequency  as 


shown  in  Ref.  8.  In  GaAs/AIGaAs  heterostructures,  how¬ 
ever,  the  free-carrier  absorption  has  not  been  measured  yet. 

In  order  to  eliminate  the  unknown  free-carrier  absorption 
and,  furthermore,  to  avoid  the  problem  of  measuring  abso¬ 
lute  FIR  intensities,  a  relative  method  is  applied  in  this  work. 

The  signals  of  two  FIR  detectors  (high-purity  n-type  GaAs 
at  35  cm"1,  Ga-doped  Ge  around  100  cm- ')  are  measured 
as  a  function  of  the  applied  electric  field.  The  signals  of  the 
two  detectors  are  evaluated  in  the  following  way:  at  a  field  Et 
(corresponding  to  a  carrier  temperature  T, )  the  detector  sig¬ 
nal  from  the  nth  detector  is  given  by 

U*.t  =  J"  (»)In[*>,T,)d*>.  (3) 

The  detector  responsivity  R„  (to)  can  be  written  as/?„M 
=  R  ®  rjat)  as  well  as  A  (<u)aS/4  (cv\  —  A  °  <j„M.  This  means 
that  the  responsivity  J7„M  and  the  absorptivity  of  the  elec¬ 
tron  system  A  M  are  mathematically  expressed  by  a  product 
of  the  relative  frequency  dependent  functions  r.(&>),  a, (to) 
and  absolute  factors  A  °  and  R  “  at  the  detector  peak  fre¬ 
quencies.  The  following  assumptions  can  be  made,  (a)  A 
Drude-type  relative  frequency  behavior  of  a.M  within  the 
frequency  interval  of  each  detector  n  is  assumed,  which  is  a 
good  approximation  even  if  the  dynamical  conductivity  de¬ 
viates  from  Drude  behavior  over  a  wide  range  of  frequency. 

(b)  The  field  dependence  of  A  (<u)  is  assumed  to  follow  the  field 
dependence  of  the  momentum  relaxation  time  r  via  the  dyn¬ 
amical  conductivity.  This  change  of  the  dynamical  conduc¬ 
tivity  plays  a  role  at  high  electric  fields  where  the  mobility 
changes  with  electric  field.  Then  Eq.  (3)  contains  two  un¬ 
known  parameters,  namely,  the  carrier  temperature  T,  and 
the  product  A°  R°,  since  the  relative  detector  response 
r„  (<u)  is  well  known.9  Measuring  the  detector  signals  UK  t  of 
the  two  detectors  with  two  different  applied  electric  fields, 
gives  a  numerically  soluble  set  of  four  equations  for  the  four  i 

unknown  variables.  Measurements  at  more  values  of  the 
electric  field  yield  more  data  than  unknown  variables  and  I 

thus  allow  us  to  prove  the  consistency  of  the  assumptions.  \ 

Thus  a  determination  of  the  electron  temperature  can  be  i 

made  from  the  relative  emission  signals,  without  knowing  ’ 

the  absolute  detector  response  and  the  optical  properties  of 
the  system.  I 
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FIO.  I .  Carrier  mobility  as  a  function  of  lateral  electric  field  for  the  two 
samples:  the  high  mobility  sample  (sample  2)  is  measured  before  (dark)  and 
after  (light)  illumination  with  visible  light. 

The  samples  used  in  the  experiments  are  modulation- 
doped  GaAs/Alo  ,5  Ga„  „  As  single  heterostructures  grown 
by  molecular  beam  epitaxy.  The  electron  mobilities  are 
1.0X10*  cm2/Vs  at  electron  concentration  n,  =  2.4 X 10" 
cm  “  2  at  4.2  K  (sample  2).  By  illumination  the  electron  con¬ 
centration  and  mobility  can  be  changed  persistently  to  n, 
=  3.8 x  10"  cm-2  and  p  »  1.3 x  10*  cmVVs  as  a  conse¬ 
quence  of  the  persistent  photoeffect  found  by  Stormer  el  al. 10 
Ohmic  contacts,  which  must  be  of  extremely  low  resistance 
in  order  to  avoid  FIR  emission  from  the  contacts,  were  made 
by  alloying  evaporated  Au/Ge  films  of  1000-A  thickness  at 
400  ‘C.  For  comparison  also  a  low  mobility  sample  ( 1 )  with 
n,  =  8.7x  10"  cm2/Vs and /i  =  1.6x  104cm2/Vs  was  used. 
All  experiments  were  performed  at  liquid  helium  tempera¬ 
tures.  The  electric  field  was  pulsed  in  order  to  avoid  sample 
heating  and  to  allow  the  use  of  correlation  techniques.  The 
current  through  the  device  and  the  detector  signals  are  mea¬ 
sured  simultaneously.  The  geometry  of  the  waveguide  and 
the  detectors  favors  the  detection  of  the  radiation  emitted  in 
the  direction  normal  to  the  surface  of  the  sample.  The  radiat¬ 
ing  area  is  2  mm  X  3  mm. 

Figure  1  displays  the  carrier  mobility  of  both  samples  as 
a  function  of  the  longitudinal  electric  field.  With  increasing 
electric  field  the  mobility  of  sample  2  decreases  due  to  elec¬ 
tron  heating  and  increased  phonon  scattering.  The  low  mo¬ 
bility  sample  ( I )  does  not  show  a  significant  field  dependence 
of  the  mobility  since  it  is  dominated  by  impurity  scattering. 
In  Fig.  2  the  obtained  values  of  the  detector  signals  as  a 
function  of  the  applied  electric  field  are  plotted.  Here  the 
characteristic  feature  of  this  experiment  can  be  seen  very 
clearly:  the  signals  of  the  two  detectors  have  remarkably  dif¬ 
ferent  slopes,  the  increase  of  the  Ge(Ga)  detector  signal  (full 
symbols)  with  electric  field  is  much  stronger  than  that  of  the 
GaAs  detector  signals.  This  is  a  consequence  of  the  exponen¬ 
tial  dependence  of  the  Bose  quantum  statistics  on  the  quanti¬ 
ty  ho/kg  T.  Therefore,  at  different  frequencies  the  tempera¬ 
ture  dependence  of  the  emission  intensity  is  different. 

From  Eqs.  ( 1 H^)  quantitative  values  of  the  electron 


FIG.  2.  Detector  signals  as  a  function  of  applied  electric  Held  for  two  sam¬ 
ples  (sample  I  and  sample  2.  dark)  and  two  different  detectors:  GaAs  detects 
the  radiation  of  longer  wavelength  (35  cm  ')  than  Ge  (Ga|  laround  100 
cm  '  \  full  symbols).  v 

temperature  are  determined  using  numerical  fitting  proce¬ 
dures.  The  results  are  shown  in  Fig.  3,  where  the  electron 
heating  AT  =  T  —  4.2  K  is  plotted  versus  the  input  power 
per  electron:  the  electron  heating  increases  with  the  input 
power  and  reaches  100  K.  when  the  input  power  exceeds  107 
eV.  s~ 1  per  electron.  The  values  of  heating  for  the  three  dif¬ 
ferent  values  of  carrier  concentration  and  mobility  are  sur¬ 
prisingly  close  together,  although  the  electron  concentration 
varies  by  a  factor  of  3  and  the  mobility  by  a  factor  of  100. 
Also  the  slopes  of  the  curves  are  similar  except  sample  2 
(light)  which  shows  a  somewhat  stronger  increase  of  heating 
with  input  power.  The  change  of  slope,  however,  is  too  weak 
and  almost  within  the  error  bars,  so  there  can  be  no  conclu- 


FIG.  3.  Electron  heating  T  -  4  2  K  as  obtained  from  the  evaluation 
of  the  emission  intensities  at  the  two  frequencies  of  the  detectors  The  ab¬ 
scissa  is  the  input  power  per  electron.  The  lines  (dashed  for  sample  2  light) 
are  drawn  only  to  guide  the  eye 
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FIG.  4.  Measured  values  of  the  free -earner  absorptivity  vs  far  infrared  wave 
number.  The  full  curve  is  calculated  assuming  Drude  behavior  with  mo¬ 
mentum  relaaation  limes  taken  from  the  dc  mobility. 

sions  drawn  from  this  behavior.  These  results  show  that  the 
heating  in  GaAs/AlGaAs  heterostructures  in  this  range  of 
lattice  temperature  and  carrier  heating  is  a  function  of  the 
input  power  per  carrier  rather  than  an  explicit  function  of 
carrier  density  and  mobility.  This  observation  is  similar  to 
the  results  obtained  in  Si  MOSFETV  and  is  also  consistent 
with  magnetoconductivity  data  of  Sakaki  et  a/.s  showing 
also  only  a  small  dependence  on  carrier  concentration  and 
mobility.  This  behavior  is  expected  when  acoustic  phonon 
scattering  determines  the  energy  loss  of  the  carriers  to  the 
lattice.5 1 1  The  absolute  values  obtained  by  this  optical  emis¬ 
sion  method  are  in  good  quantitative  agreement  with  the 
results  of  photoluminescence  experiments  by  Shah  et  al *  at 
high  electron  temperatures  and  the  values  obtained  by  Sa¬ 
kaki  et  al .*  from  magnetoconductance  measurements  at  low 
electric  fields. 

A  comparison  of  the  detector  signals  with  those  of  n- 
InSb  cyclotron  emission  at  saturation  electron  tempera¬ 
tures12  allows  an  absolute  determination  of  the  emission  in¬ 
tensity  as  described  in  Ref.  8.  Detector  responsivities  of 
~  1  X  107  V/W  (GaAs)  and  ~  2  X 107  V/W  (Ge)  with  an  ac¬ 
curacy  of  30%  are  obtained  for  the  used  detectors.  The  emit¬ 
ted  power  from  the  GaAs/AIGaAs  heterostructures  there¬ 
fore  ranges  from  10” 14  W,  which  represents  the  noise  level  of 
the  detector  arrangement,  up  to  10" 10  W  within  the  detector 
intervals.  Knowing  the  absolute  emission  intensity  and  the 
electron  temperature,  we  can  separate  the  product  Rj 
and  give  values  of  the  absorptivity  of  the  2D  carrier  system  at 
the  detector  frequencies.  The  results  are  shown  in  Fig.  4.  The 
absorptivity  of  the  carrier  plasma  is  plotted  for  the  different 
samples  at  the  two  detector  frequencies.  The  free-carrier  ab¬ 
sorption  (evaluated  for  low  electric  fields)  strongly  depends 
on  the  carrier  mobility;  the  absorption  varies  by  two  orders 
of  magnitude  from  the  low  mobility  sample  (high  absorptiv¬ 
ity  =  emissivity)  to  the  high  mobility  sample  which  reveals 


a  very  low  absorptivity  of  below  10_:>at  100cm-1.  Note  that 
these  low  values  of  absorption  are  measurable  without  gate 
modulation  techniques  as  usually  applied  for  absorption 
measurements  in  2D  systems.7  The  full  curve  in  Fig.  4  dis¬ 
plays  the  theoretical  frequency  and  mobility  dependence  as¬ 
suming  Drude  conductivity  with  the  values  of  the  relaxation 
time  t  taken  from  the  dc  mobility.  The  agreement  (within  the 
experimental  error  range  coming  from  the  uncertainty  in 
measuring  the  absolute  detector  responsivity)  obviously 
shows  that  the  FIR  absorption  in  the  2D  plasma  is  deter¬ 
mined  by  the  same  scattering  processes  as  the  dc  transport. 
This  means  that  Coulomb  scattering  by  residual  impurities, 
which  is  the  dominant  scattering  process  in  these  structures, 
also  provides  the  numienltini  transfer  for  the  free-carrier  op¬ 
tical  transitions. 

In  conclusion,  we  have  observed  broadband  far  infrared 
radiation  from  the  hot-carrier  plasma  in  GaAs/AIGaAs  he¬ 
terojunctions  with  ultrahigh  as  well  as  low  electron  mobility. 

A  relative  method  of  determining  the  carrier  heating  in  the 
electric  field  is  introduced,  which  allows  the  measurement  of 
T  over  a  wide  range  of  carrier  heating  ( 10  K  <  T<  1 50  K).  The 
method  can  be  also  applied  to  investigate  carrier  heating  in 
small  devices  and  integrated  circuits.  At  room  temperature, 
however,  one  has  to  work  in  the  near  infrared  range  in  order 
to  obtain  a  similarly  strong  dependence  of  the  emission  in¬ 
tensity  on  carrier  heating.  Furthermore,  the  FIR  absorption 
of  the  2D  plasma  could  be  estimated  for  the  first  time,  show¬ 
ing  a  strong  dependence  on  the  amount  of  impurity  scatter¬ 
ing  which  permits  the  momentum  transfer  for  the  optical  * 
transition. 
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The  far  infrared  (FIR)  broadband  emission  from  the  two-dimensional  (2D) 
electron  plasma  in  GaAs/AlGaAs  heterostructures  is  investigated  for  the 
first  time  by  means  of  narrowband  detectors.  The  carrier  heating  in  the 
electric  field  and  the  free  carrier  optical  properties  are  measured.  In 
low  mobility  samples  the  free  carrier  absorption  can  be  described  by  a 
Drude  behaviour  calculated  from  the  DC  electrical  conductivity,  in  high 
nobi.l  samples  the  free  carrier  absorption  extremely  decreases  as  ex¬ 
pect'd,  leviations  from  the  Drude  behaviour,  however,  are  observed. 

Free  carrier  optical  transitions  in  an  external  periodical  potential 
(periodic  gate  structure)  are  calculated,  experimental  attempts  on  FIR 
emission  from  such  structures  are  discussed. 


I.  FIR  Broadband  Emission  from  a  "warm"  2D  Plasma 


In  the  thermodynamical  equilibrium  the  spectral  emission  Intensity  X (co) 
of  electromagnetic  radiation  of  a  system  with  temperature  T  is  given  by 

I  (u)  ■  IjjB  (w)  •  A  (u)  ( 1 ) 

where  IBB  is  the  spectral  emission  intensity  of  a  black  body  and  A(u>)  is 
the  spectral  absorptivity  of  the  system,  in  a  2D  system  given  by 

AM - 1^252-5  ,  '  (21 

with  F  »  o  (u>)  / cQc  and  e  being  the  dielectric  constant  of  the  substrate. 

A  Drude  form  of  the  dynamical  conductivity  c(w)  in  n-Si  inversion  layers 
has  been  confirmed  by  microwave  and  FIR  transmission  experiments  per¬ 
formed  by  Allen  et  al.  /I/  up  to  frequencies  of  40  cm-1.  In  GaAs/AlGaAs 
heterostructures,  however,  the  free  carrier  absorption  has  not  been  mea¬ 
sured  yet.  In  order  to  evaluate  FIR  emission  intensities  from  the  free 
carrier  plasma  -  as  it  has  been  performed  in  Si-MOSFET  structures  in  ref. 
/2/  -  a  relative  method  is  applied  in  this  work.  The  signals  of  two 
narrowband  detectors  (high  purity  n-type  GaAs  at  35  cm-1,  Ga-doped  Ge  at 
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100  cm-1)  are  measured  as  a  function  of  the  electric  field.  From  the  re¬ 
lative  signals  the  electron  temperature  is  determined  using  a  fit  to 
equation  (1),  from  the  absolute  intensities  the  values  of  the  free  car¬ 
rier  absorption  coefficient  can  be  obtained. 

The  samples  used  in  our  experiments  are  MBE  grown  GaAs-AlGaAs  single  he- 

6  2 

terostructures  with  mobilities  of  1x10  cm  /Vs  at  electron  concentration 
11  -2 

n  *  2.4x10  cm  at  4.2  K.  By  LED  illumination  the  electron  concentre- 

S  11  —2 
tion  and  mobility  changes  persistently  to  n.  =  3.8x10  cm  and  u  * 

6  2  s 

®  1.3x10  cm  /Vs.  Ohmic  contacts,  which  are  important  to  be  of  extremely 
low  resistance  in  order  to  avoid  FIR  emission  from  the  contacts,  were 
made  by  alloying  Au/Ge  evaporated  films  of  1000  A  at  400°C.  All  experi¬ 
ments  were  performed  at  liquid  helium  temperatures,  experimental  details 
are  contained  in  ref.  /2/,  the  detectors'  spectral  responsivities  are 
described  in  ref.  /3/. 

The  experimental  results  are  shown  in  fig.  1  and  fig.  2.  With  increasing 
electric  field  the  mobility  decreases  due  to  electron  heating  and  in¬ 
creased  phonon  scattering.  The  evaluation  of  mobility  data  for  electron 
temperatures  -  as  shown  in  ref.  /4/  by  Shah  at  al.  in  comparison  to 
pnotolumine scenes  experiments  -  is  somewhat  problematic  due  to  sample  de¬ 
pendent  offects.  Therefore  this  decrease  of  mobility  with  carrier  heating 
is  not  quantitatively  evaluated. 

The  values  for  electron  heating  as  a  function  of  the  input  power  eu [I?  !2 
per  electron,  as  obtained  from  the  relative  FIR  emission  intensities  at 
two  frequencies,  is  plotted  in  figure  2  for  two  values  of  the  carrier 
concentration.  It  is  emphasized  that  no  absolute  intensities  have  to  be 
known  in  order  to  evaluate  the  carrier  heating  from  the  relative  signals 
of  the  GaAs-  and  Ge-detector,  respectively.  The  carrier  heating  above 
liquid  helium  temperature  is  in  the  same  range  as  measured  by  photolumi¬ 
nescence  by  Shah  et  al.  /4/,  where  the  carrier  distribution  in  multi¬ 
quantum  wells  is  directly  measured.  At  low  electric  fields  the  carrier 
heating  is  more  effective  for 'the  lower  electron  concentration  -  which 
also  can  be  seen  in  the  mobility  behaviour  qualitatively  (fig.  1).  At 
higher  fields  the  behaviour  changes  and  the  two  lines  cross. 

—13  — 1 2 

From  the  absolute  emission  intensities  (10  -10  W)  ,  which  were  esti¬ 

mated  by  comparison  with  InSb  cyclotron  emission  at  saturation  electron 
temperatures  -  see  ref.  /2/  for  details  -  first  informations  about  the 

free  carrier  optical  properties  can  be  obtained:  The  experiments  show 

4  5  2 

that  the  free  carrier  absorption  in  low  mobility  samples  (u^lO  -  10  cm /Vs, 
—  1  —3 

A  (at  v*35cm  )  ^10  )  can  be  quantitatively  described  by  a  Drude  be¬ 

haviour  using  the  T-values  from  DC  electrical  transport.  The  extremely 
high  mobility  samples  (fig.  1  and  2)  reveal  extremely  low  absorptivity 


i 
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4.  FIR-Emission  from  Free  Carrier  Plasma  in  GaAs/AlGaAs  Heierostructures 
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Fiq.  Is  Electron  mobility  determined  Fig.  2:  Electron  heating  AT  »  Te  - 
from  DC  transport  as  a  func-  -4.2K  determined  from  re- 

tion  of  the  electric  field.  lative  FIR  emission  as  a 

function  of  the  input  pcwer 
eu|Ei’  per  electron. 

values  below  10  4 (at  ~ « 35  cm  1 )  and  below  10-^  (at  v«100  cm-1).  Note, 
that  these  low  values  of  absorption  are  measurable  in  emission  experi¬ 
ments  without  using  gate  modulation  techniques!  The  dependence  of  the 
emission  intensity  on  the  electric  field,  however,'  cannot  be  described 
by  a  Drude  form  of  a(w)  with  the  relaxation  time  x  taken  from  the  field- 
dependent  mobility.  A  quantitative  fit  is  obtained  only,  when  a  constant, 
not  field-dependent  value  of  x  (in  this  case  the  low-field  momentum  re¬ 
laxation  time  xffl)  is  used.  The  decrease  of  xm  at  higher  electric  fields, 
as  seen  in  mobility  data,  is  not  relevant  for  the  free  carrier  emissivity. 
From  this  we  conclude  that  free  carrier  optical  transitions  are  made  pos¬ 
sible  mainly  by  impurity  scattering,  which  enables  the  momentum  transfer. 

II.  Free  carrier  optical  transitions  in  an  external  periodic  potential 

Two-dimensional  electron  systems  represent  ideal  systems  for  observing 
optical  transitions  due  to  periodically  accelerated  motion  (*Smith- 
Purcell-effect"  /5/) :  (1)  The  carriers  are  confined  near  the  surface  and 
can  be  brought  close  to  a  periodic  gate  potential;  (2)  the  carrier  system 
is  extremely  degenerate  and  high  Fermi  velocities  are  present  at  simul¬ 
taneously  high  mobilities  giving  extremely  long  mean  free  paths  up  to 
more  than  1  urn  at  the  Fermi  surface. 

In  the  following  the  emission  intensity  of  a  GaAs/AlGaAs  heterostructure 
due  to  the  free-carrier  radiation  by  the  motion  of  the  electrons  with  ve¬ 
locity  v  near  the  Fermi  surface  along  a  periodic  potential  with  amplitude 
AU  and  period  (a)  is  estimated  in  a  semi-classical  model  using  the  dipole 
radiation  theory  and  Fermi- statistics  for  the  transition  probability. 
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The  emission  frequency  is  given  by  w  »  (2ir/a)  .vcosj,  for  the  estimation 
of  the  spectral  emission  intensity  the  following  assumptions  are  made: 

(a)  The  drift  of  the  carriers  along  the  grating  is  neglected,  only  the 
thermal  velocity  is  considered,  (bj  only  the  part  of  the  carriers  within 
the  polar  angle  interval  {  — tt/4  ,  +it/4)  with  nearly  angle-independent  emis¬ 
sion  frequency  (cos$  ^  1 )  is  taken  into  account.  Thus  the  spectral  emis¬ 
sion  intensity  per  unit  area  is  given  by 


I  (u)  da  -  e4(MJ)2-a,-,f(E)  [l-f  (t_-  h,)]  ^  _  (2) 

6e0  c3h2 

In  fig.  3  the  spectral  absolute  intensity  of  FIR  emission  from  a  real 
GaAs/AlGaAs  heterostructure  is  plotted  (ng  «  4x1011cm-2,  a  «  2200  A). 

The  spectral  intensity  per  frequency  interval  is  for  a  realistic  periodic 
potential  of  40  »  1  mV  in  the  order  of  some  10  W/cm  cm  and  should 
therefore  be  observable.  The  intensity  is  proportional  to  the  square  of 
the  periodic  potential,  and  thus  may  reach  values  of  up  to  uW/cn2cm-1. 

By  assuming  only  carriers  with  Jc-orientation  near  to  the  grating  orienta¬ 
tion,  we  neglect  the  low-frequency  emission  of  carriers  moving  in  small 
angles  along  the  grating.  Thus  the  real  emission  intensity  should  be 
higher  than  estimated  here.  A  classical  treatment  of  this  contribution 
would  give  higher  values  than  the  black-body  emission  in  the  low-frequency 
regime  and  thus  break  the  thermodynamic  equilibrium  model. 


Experimentally  the  system  is  realised  by  evaporating  a  semitransparent 
(NiCr)  gate  electrode  upon  a  sinusoidal  photoresist  structure  and  apply¬ 
ing  small  voltages  and  a  "compensation"  voltage  along  the  gate  in  order 
to  avoid  pinch-off  effects  during  the  electric  field  heating  of  the  2D 
system  as  described  in  ref.  /2/.  Experimental  results  on  this  new  source 


will  be  published  in  a  forthcoming  paper. 


Spectral  emission  intensity 
of  "Smith-Purcell"-radia- 
tlon  in  GaAs/AlGaAs  assum¬ 
ing  a  periodic  potential 
&U  ■  1  mV  and  taking  only 
into  account  electrons  in 
the  interval  (  —tt/4  ,  ♦n/4) 
in  respect  to  the  grating 
orientation.  Frequency 
interval  1  cm”  . 
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We  report  the  first  measurement  of  inter-Landau  level  lifetimes  in  GaAs/GaAl As-heterostructures . 
Saturation  cyclotron  resonance  has  been  measured  using  a  high  intensity  optically  pumped  far 
infrared  laser.  The  results  are  analyzed  on  the  basis  of  a  three  level  model.  The  obtained  life¬ 
times  vary  from  0.18  to  1.2  ns  for  samples  with  different  carrier  concentration. 


1.  INTRODUCTION 

In  the  recent  past  high  power  far  Infrared 
spectroscopy  has  been  used  to  determine  the  life¬ 
times  of  impurity  and  Landau  states  in  semi¬ 
conductors1'^.  Up  to  now  these  investigations 
have  been  restricted  to  bulk  material.  Another 
method  to  determine  energy  relaxation  times  was 

4 

employed  by  Bauer  and  Kahlert  for  degenerate 
bulk  material  and  recently  by  Sakaki  et  al5  for 
GaAs/GaAl As-heterostructures:  The  energy  relaxa¬ 
tion  time  in  a  magnetic  field  was  deduced  from 
the  strength  of  the  Shubnikov-de  Haas  (SdH)  os¬ 
cillations  as  a  function  of  the  electric  field. 

It  was  found  by  Gornik  et  al.1  and  Allan  et 

3 

al.  that  in  bulk  material  electron-electron 
scattering  plays  a  dominant  role  for  the  inter- 
Landau  level  (LL)  lifetime.  A  connection  between 
electron  density  in  the  excited  level  and  life¬ 
time  was  established. 

The  aim  of  the  present  paper  is  to  find  out 
whether  a  similar  mechanism  controls  the  relaxa¬ 
tion  in  20-layers  as  present  in  heterol ayers . 

We  have  performed  cyclotron  resonance  (CR)  ab¬ 
sorption  saturation  measurements  on  GaAs/GaAlAs- 


heterostructures.  This  technique  is  a  direct 
method  for  probing  the  inter-LL  lifetime.  The 
absorption  process  is  described  on  the  basis  of 
a  three  level  model.  By  comparison  with  the  ex¬ 
perimental  data  the  lifetime  is  deduced  and 
found  to  vary  from  0.18  ns  to  1.2  ns  for  samples 
with  different  carrier  concentrations  indicating 
the  Influence  of  an  electron-electron  scattering 
process  also  in  the  2D  situation.  A  direct  com¬ 
parison  with  theoretical  calculations  is  diffi¬ 
cult,  since  most  of  the  calculations  concerning 
/ 

energy  relaxation  in  2D-systems  do  not  include 
a  magnetic  field6*8. 

2.  EXPERIMENTAL 

The  samples  were  grown  by  molecular  beam 

epitaxy,  their  mobilities  vary  between  105  and 
£  2 

10  cm  /Vs,  their  carrier  concentrations  between 
1 1  -2 

1.2  and  3-10  cm  .  A  high  power  quasi  cw  opti¬ 
cally  pumped  FIR  laser  was  used  in  the  pulsed 
mode.  The  pulse  duration  (0.3  ms)  was  adjusted 
to  be  much  longer  than  the  expected  relaxation 
time,  ensuring  steady-state  conditions  during 
optical  excitation.  In  the  present  work  only  the 
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118.8  pm  methanol  line  was  used,  because  it 
yielded  the  highest  intensities.  Moreover,  at 
the  corresponding  energy  of  10.43  meV  the  non- 
parabol icity  and  polaron  effect  are  large 
enough  to  enable  saturation.  The  FIR  intensity 
in  the  sample  was  determined  by  a  calibrated 
pyroelectric  detector  taking  into  account  the 
losses  of  the  light  pipe  and  the  reflection  at 
the  substrate  surface.  Due  to  the  light  pipe 
optics  the  sample  was  subjected  to  unpolarized 
radiation.  Therefore  the  effective  intensity  in 
the  sample  was  one  half  of  the  total  intensity 
in  the  sample  as  the  geometry  was  rather  a  Fara¬ 
day  geometry.  An  n-type  InSb  photoconductive 
detector  was  used,  situated  below  the  sample  in 
the  liquid  He-bath.  The  electron  concentration 
was  determined  from  SdH-measurements. 

3.  THEORETICAL  MODEL 

Saturation  of  cyclotron  resonance  transitions 

is  not  possible  in  a  equidistant  Landau  ladder. 

However,  if  a  sufficient  amount  of  nonparaboli- 

city  and  polaron  contribution  to  the  effective 

mass  is  present,  the  level  shifts  can  be  larger 

than  the  linewidth  of  the  individual  transitions. 

In  this  case  saturation  can  be  achieved.  Recent 
g 

experiments  have  shown  that  the  Influence  of 
the  LO-phonon  (hu  u  36  meV)  on  the  energy  levels 
In  GaAs/GaAlAs-heterostructures  is  negligible 
below  25  meV,  but  becomes  strong  above  this 
value.  In  the  present  situation  the  photon  ener¬ 
gy  is  hu>  «  10.43  meV,  which  means  that  the  LL 
n  >  0,  1  and  2  undergo  small  polaron  shifts, 
whereas  the  n  ■  3  LL  is  shifted  considerably. 
Thus  for  a  not  too  broad  absorption  line  the 
analysis  has  to  include  up  to  three  levels.  If 
the  absorption  linewidth  is  smaller  than  the 
nonparabol icity  (T  <  (hwc)  /eg),  only  two  levels 
need  to  be  considered. 

The  rate  equations  can,  therefore,  be  written 


df.  1 


co^1“r(vAi)*  rv’-v-vi-v:- (2) 

at  n»J  T 


df,  I 

°o  — #  *  7“  A1  '  f2(1  '  f1}  ~ 

dt  hu  t 


with  I  the  effective  laser  intensity  in  the 
sample,  oQ  *  eB/rh  the  electron  concentration 
per  LL,  and  x  the  relaxation  time,  x  is  chosen 
to  be  the  same  for  both  the  n  *  1  and  the  n  «  2 
LL  since  the  experimental  information  is  not 
enough  to  distinguish  a  separate  x1  and  ~2' 

An  is  the  absorptivity  of  a  2D  electron  gas 
given  by 

,  .  u. 

""  (^♦1  +  Reo(n,/c0c)2  assumed) 


with  Reo^nVe0c  » 

c-K(VW(n  +  ,) 


i  n 


the  real  part  of  the  dynamical  conductivity. 

C  is  a  dimensionless  constant  determined  from 
the  value  of  the  low-intensity  absorption.  The 
/i  in  the  numerator  of  Eq.  (4)  is  due  to  the 
experimental  set-up  since  the  laser  beam  passes 
the  GaAs  substrate  prior  to  the  2D  system.  The 
usual  expansion  of  Ar  for  Reo^/EgC  «  1  is 
not  employed,  because  the  changes  in  trans¬ 
mission  reach  65 %  in  the  CR-active  polariza¬ 
tion.  In  steady  state,  Eq.  ( 1 )  -  (3 )  are  solved 
numerically  together  with  the  condition 


f0  ♦  f  r  f2  •  v  ( 

where  v  is  the  number  of  (spin  degenerate) 
filled  LL,  v  *tinsh/eB.  For  a  comparison  with 
the  experiment  the  normalized  transmission 
change  1-T/T(o  =  0)  is  calculated,  with 
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4  .c 

T  1  .  r^’W  <6> 


with  the  best  fit  from  Eq.  (1)  - (5)  for  sample 
1360. 


4.  RESULTS  AND  DISCUSSION 

Typical  experimental  traces  are  shown  in 
Fig.  1  for  sample  1433  (ns  *  1.2  •  lO^cnf2)  and 
sample  1360  (ns  .  3  •  lO^cm"2). 


FIGURE  1 

Cyclotron  resonance  absorption  for  two  different 
GaAs/GaAlAs-heterostructures  at  118.8  pm  for 
several  laser  Intensities. 

(a)  sample  1433  (b)  sample  1360 

The  transmission  zero  corresponds  to  CR-active 
polarization. 

Two  features  can  be  clearly  seen:  In  both  samples 
saturation  has  been  achieved,  but  more  easily  in 
the  low  concentration  sample.  Furthermore,  a 
shift  of  several  hundred  Gauss  is  observed  only 
in  Fig.  lb,  indicating  a  transfer  of  electrons 
to  the  n  »  2  LL,  which  also  makes  the  saturation 
more  difficult.  This  clearly  shows  the  require¬ 
ment  of  including  this  level  into  the  analysis 
as  described  before.  In  Fig.  la  the  linewidth  is 
so  small  that  there  is  no  overlap  between  the 
two  transitions  and  accordingly  no  shift. 

In  Fig.  2  the  intensity  dependence  of  the 
relative  transmission  change  is  plotted  together 


FIGURE  2 

Intensity  dependence  of  the  relative  trans¬ 
mission  change  for  a  laser  line  of  118.8  pm 
for  sample  1360.  The  solid  line  represents  the 
best  fit  from  Eq.  (1)  -  (5). 

An  onset  of  saturation  is  observed  with  a  re¬ 
duction  of  1  -  T/T  up  to  30  J.  A  further  in- 
crease  in  intensity  was  limited  by  the  experi¬ 
mental  system. 

Table  1  gives  the  resulting  lifetimes  for  all 
samples  Investigated: 


TABLE  1 


sample 

ns(10ncm'2) 

x(ns) 

1360 

3.0 

0.18 

1355 

3.0 

0.25 

1234 

2.0 

0.65 

1395 

1.5 

0.55 

1433 

1.2 

1.20 

The  lifetimes  were  all  obtained  under  conditions, 
where  the  n  *  0  LL  is  just  filled  or  partially 
filled. 

Within  the  experimental  error  the  lifetime 
increases  with  decreasing  carrier  concentration. 
In  an  earlier  paper  Sakaki  et  al.^  reported 
somewhat  larger  t- values,  but  the  same  depen- 
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dence  on  carrier  concentration.  This  was  ex¬ 
plained  by  a  concentration  dependent  acoustic 
phonon  rate.  His  experiment,  however,  was  per¬ 
formed  in  the  weak  heating  regime  (T  <  30  K)  and 
with  small  magnetic  fields  (B  <  3T).  In  the 
present  work  a  considerable  part  of  electrons 
is  excited  to  higher  LL,  corresponding  to  inter¬ 
subband  temperatures  well  above  100  K.  In  addi¬ 
tion,  for  a  LI  separation  of  10  meV  the  acoustic 
phonon  emission  is  most  likely  not  the  dominant 
relaxation  process.  The  influence  of  an  elec¬ 
tron-electron  scattering  process  similar  to  that 
found  in  bulk  material1,3  seems  to  be  evident 
for  the  observed  density  dependence.  To  examine 
this  process  in  detail  the  excited  level  carrier 
concentration  has  to  be  determined.  Up  to  now 
there  are  not  enough  data  available  reaching  out 
in  the  strong  saturation  regime  to  make  a  mean¬ 
ingful  evaluation  about  the  density  dependence. 

We  found  that  the  obtained  lifetimes  depend 
Strongly  on  the  analysis  of  the  data.  The  factor 
(1  -  f^ )  In  Eq.  (1 )  -  (3)  gives  rise  to  a  much 
slower  decrease  in  the  saturation  curve  (Fig.  2) 
and  to  an  earlier  threshold  for  saturation,  es¬ 
pecially  in  the  case  of  a  filled  LL.  Neglecting 
this  factor  would  result  in  a  longer  lifetime 
up  to  a  factor  of  five.  From  this  we  conclude 
that  a  reliable  evaluation  needs  more  data  and 
data  to  higher  power  levels,  where  this  factor 
plays  a  less  important  role. 

A  large  uncertainty  is  found  in  the  eval uation 
for  the  low  density  samples  1395  and  1433  caused 
by  experimental  problems:  The  magnitude  of  the 
low-power  absorption  was  not  well  reproducible 
according  to  the  unstable  experimental  condi¬ 
tions.  This  seems  to  come  from  the  flat  inver¬ 
sion  channel  in  these  samples.  Similar  effects 

in 

have  been  observed  with  DC  fields  ”.  It  was, 
therefore,  difficult  to  determine  the  constant 
C  in  Eq.  (4)  accurately. 

In  summary,  we  have  used  FIR  saturation  spec¬ 
troscopy  to  determine  inter-Landau  level  life¬ 


times  in  GaAs/GaAlAs-heterostructures.  We  have 
found  clear  evidence  for  an  electron  concentra¬ 
tion  dependent  lifetime  ranging  from  0.18  ns  to 
1.2  ns.  This  suggests  that  electron-electron 
scattering  is  playing  a  dominant  role  in  the 
energy  loss  mechanism.  A  detailed  description 
of  the  process  needs  more  experimental  and 
theoretical  work. 
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Landau-level  population  inversion  in  crossed  electric  and  quantizing  magnetic  fields 
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The  problem  of  the  electron  distribution  in  crossed  electric  and  quantizing  magnetic  fields  is 
treated  for  free  electrons  in  a  parabolic  band.  The  master  equation  for  the  diagonal  elements  of 
the  density  matrix  is  solved,  neglecting  the  k,  dependence  of  the  electron  distribution.  It  is  shown 
that  population  inversion  between  Landau  levels  can  be  achieved  if  scattering  by  optical  phonons 
and  a  quasielastic  scattering  process  are  taken  into  account. 


It  is  well  known  that  in  crossed  electric  and  magnetic 
fields  the  electron  distribution  function  can  be  inverted,1*9 
when  the  interaction  between  electrons  and  optical  pho¬ 
nons  is  strong  enough.  The  electron  motion  is  usually 
described  in  the  picture  of  “streaming  motion":4*9  The 
electrons  are  assumed  to  be  accelerated  until  they  have 
enough  energy  to  emit  an  optical  phonon;  after  that  they 
return  into  the  so-called  passive  region  (inside  the  circle 
o<v^-y/2ha>oj  m* ).  If  they  are  scattered  into  the  part 
of  the  passive  region  where  the  trajectories  are  closed  (ac¬ 
cumulation  area),  they  will  stay  there  until  they  are  re¬ 
moved  by  acoustic-phonon  or  impurity  scattering.  The 
motion  is  characterized  by  the  parameter  (—u0 pB/E.  For 
{<  l  all  electrons  are  streaming,  and  for  (>  2  all  elec¬ 
trons  are  accumulated  (more  precisely,  the  main  trajectory 
through  the  point  |  v  |  —0  is  closed),  while  for  1  <  { <  2 
the  electrons  are  partly  streaming  and  partly  accumulated. 

It  was  shown  that  for  £S 2  a  torus-type  distribution 
function  can  arise  in  the  accumulation  area,  which  corre¬ 
sponds  to  a  population  inversion  between  Landau  levels 
(LL’s).2,3,4,7  Theoretical  studies  based  on  Monte  Carlo 
calculations7  as  well  as  analytical  approaches  based  on  the 
Boltzmann  equation,4  were  carried  out.  The  accumulation 
of  carriers  in  this  area  was  proved  experimentally  in  dif¬ 
ferent  materials,4  but  it  was  not  clear  if  a  population  inver¬ 
sion  between  LL’s  could  be  achieved.  Very  recently,  Che- 
botarev  and  Murzin  reported  the  observation  of  a  strong 
emission  signal  in  n- type  Ge,  which  was  ascribed  to  a  LL 
inversion. 

Up  to  now,  to  our  knowledge,  all  theoretical  studies 
were  performed  in  the  limit  of  nonquantizing  magnetic 
fields;  that  means  the  use  of  a  classical  distribution  func¬ 
tion  smeared  out  over  the  LL’s  and  the  use  of  the  zero- 
magnetic-field  scattering  rates.  On  the  other  hand,  if  one 
wants  to  study  optical  transitions  between  (inverted)  LL’s, 
the  discrete  structure  of  the  energy  levels  is  of  fundamen¬ 
tal  importance.  In  this  Rapid  Communication  we  treat  the 
problem  of  electrons  in  crossed  fields  in  a  fully  quantum 
mechanical  frame,  which  means  that  we  start  from  the 
Liouville  equation  for  the  density  operator.11,12 

We  consider  the  Hamiltonian  (with  obvious  notations) 

,  (I) 

where  the  electric  field  is  incorporated  in  H *  The  elec¬ 
tronic  eigenvalues  and  eigenfunctions  ( B  in  z  direction,  E 


in  x  direction)  are  given  by12 
A2*2 

Etmm(n  +  -^)haie  + - f-+e£.Tr+-|-m• 

2 m* 


£ 
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.  (2) 


where 


hky  m*£ 
eB  eB 2 


4\(r )  -const  x  elk'reU,I<t>m 


x-Xr 
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(3) 


where  l  m‘yJh/eB  and  d>„[(jr  —Xv)/l)  is  the  nth  harmonic 
oscillator  wave  function.  The  index  v  summarizes  the 
quantum  numbers  n,  ky,  k,.  The  stationary  equation  for 
the  density  operator  is 


lH,p]  -0  . 


(4) 


The  electron  distribution  function  corresponds  to  the  diag¬ 
onal  matrix  elements  p,— Tr^vlpl  v),  where  the  trace  is 
taken  over  the  phonon  states.  In  the  case  of  weak 
electron-phonon  coupling  it  is  possible  to  express  the  non¬ 
diagonal  elements  by  diagonal  elements,  which  then  obey 
the  master  equation12 

2>/»Vr-pr»V-0  .  (5) 

¥ 

The  W „<  are  the  transition  probabilities  from  state  v  to 
state  v'  induced  by  the  electron-phonon  interaction.  Equa¬ 
tion  (5)  is  a  system  of  n  coupled  integral  equations  (when 
n  LL’s  are  taken  into  account)  for  the  density-matrix  ele¬ 
ments.  As  we  look  for  a  homogeneous  distribution  func¬ 
tion,  we  assume  that  the  pr  are  independent  of  X,  and, 
therefore,  of  kr ,M2  Then  the  ky  integration  can  be  per¬ 
formed  directly,  and  we  are  left  with  k,  as  the  only  in¬ 
tegration  variable. 

At  this  point  we  make  a  crude  but  not  unrealistic  ap¬ 
proximation  to  enable  a  simple  solution.  We  drop  the 
quantum  number  ktt  which  means  that  all  electrons  are 
concentrated  at  k,  —0  (according  to  the  peak  in  the  densi¬ 
ty  of  states).  This  assumption  is  justified  by  an  experi¬ 
ment  of  Eaves  el  al. 12  who  showed  that  in  high  electric 
fields,  transitions  with  a  change  in  k,  become  less  impor¬ 
tant  than  transitions  associated  with  a  change  of  the  spa¬ 
tial  coordinate  x  (see  below  and  Fig.  1).  However,  we 
want  to  point  out  that  the  aim  of  this  paper  is  not  a  de- 
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FIG.  I.  Two  LL's  are  plotted  vs  the  direction  of  the  electric 
field.  The  mechanism  of  acoustic-  and  optical-phonon  transi¬ 
tions  is  illustrated. 


tailed  quantitative  study,  but  rather  the  first  attempt  to  get 
information  about  the  qualitative  behavior  of  the  distribu¬ 
tion  function  when  the  magnetic  field  is  quantizing.  As  a 
result,  Eq.  (5)  is  changed  into  a  system  of  purely  algebraic 
equations. 

In  the  following  we  want  to  deal  with  three  LL's,  which 
is  the  minimum  number  necessary  to  study  streaming 
motion  and  population  inversion:  two  “unbroken"  LL’s 
(n  “0,1 )  inside  the  circle  o  “Uop  and  the  n  “2  LL,  which 
is  broken  by  the  emission  of  an  optical  phonon.  The  sys¬ 
tem  of  three  linear  equations  is  solved  and  yields  po,  p\,  pi 
as  a  function  of  the  ^-integrated  transition  ratds  Wm-. 
We  take  into  account  polar  optical-phonon  scattering  and 
acoustic  deformation  potential  scattering;  thus  Wm- 
“  W'JJ- +  The  transition  probabilities  are  given  by 
Fermi's  golden  rule. 11,12  It  should  be  noted  that  the  elec¬ 
tric  Held  only  appears  in  the  energy  S  function  and  not  in 
the  squared  matrix  elements. 

In  order  to  enable  an  analytical  evaluation  of  the 
scattering  rates  we  consider  only  spontaneous  emission  of 
optical  phonons  (/V^saO)  and  apply  the  equipartition  ap¬ 
proximation  for  acoustic  phonons  (I  ). 

Further,  we  regard  the  scattering  by  acoustic  phonons  as 
elastic  These  approximations  can  be  assumed  to  be  valid 
at  intermediate  temperatures  (T®  20-80  K).  The  expli¬ 
cit  transition  rates  are  then  given  by 

W'R-*'ff-**el8+8Be2]e'*?/I  , 

H'JS  -WSL  -^“16+4(2®, )2+2<B,)4J«_a*',Vj  ,  (6) 

W'!?-^f5-W“l7+9®2-5SJe4+®‘]e'*',/2  , 

W'JU  -Aa'  te~<au,*§'m  , 

^£«“/<*l2+2(Bu>*2Bc)2]e_(^:’:2*')V2  ,  (7) 

W'jT.ia  --4*111  -bfBto*  Bf)2 

♦  (Su,*®,)4]*"^*^  , 

where 


Bi n “ p/eEl  , 

A~-  D2k,TB  'Jlxe1a>\£>Bl  i _ j_ 

li'Jlxdpc} ^El  ’  \6dhEAxso  <•  i« 


(D  . . .  deformation  potential,  p . . .  mass  density,  c, . . . 
sound  velocity,  o>lo  -  •  LO-phonon  frequency,  eu,£- . . . 
static  and  high  frequency  dielectric  constant,  respectively). 
d  is  the  thickness  of  the  two-dimensional  slab  (it  cancels 
out  in  the  final  expressions  for  pi). 

Obviously,  the  acoustic  rates  depend  critically  on  the  ra¬ 
tio  hctJeEl.  This  quantity  characterizes  the  possibility  of 
phonon  transitions.  As  long  as  eEL<Kho)e,  the  spatial 
overlap  between  wave  functions  of  different  LL’s  and  same 
energy  is  too  small  to  enable  transitions  >0).  In¬ 

stead  of  this,  transitions  associated  with  a  change  in  k„ 
which  are  not  included  in  our  model,  would  become 
relevant.  However,  no  population  inversion  is  expected  in 
this  regime,  but  we  must  be  aware  of  the  fact  that  for 
eEl  <  ho>c  out  model  is  no  longer  valid.  When  eEl 
>  ha>c,  the  exponential  in  (6)  tends  to  unity,  and  direct 
transitions  between  the  bottom  of  different  LL’s  become 
possible.  In  Fig.  I  this  mechanism  of  phonon-induced 
transitions  between  LL’s  is  shown.  Transitions  caused  by 
optical  phonons  also  yield  a  vertical  contribution  in  this 
scheme,  corresponding  to  ho>LO.  The  required  overlap  de¬ 
pends  on  the  values  of  o>io  and  <oc  and  is  much  larger  for 
transitions  to  a  higher  LL.  Hence,  H'j’P  >  IV*  (n'  >  n )  as 
expected  for  thermodynamic  reasons.  In  a  realistic  situa¬ 
tion  a>e<&a>LO'  so  the  optical-phonon-induced  transitions 
set  in  at  higher  electric  fields  than  the  acoustic  transitions. 

At  this  point  it  is  meaningful  to  analyze  the  physical  sit¬ 
uation  based  on  a  three-level  scheme.  Population  inversion 
should  build  up  between  the  two  lower  levels,  so  the  system 
has  to  be  pumped  from  n  “0  to  n  “2.  The  main  require¬ 
ments  for  the  transition  rates  are  W1{  >  W'20,  and  Wm  as 
large  as  possible  in  order  (o  ensure  a  fast  depletion  of  the 
lowest  level  (more  precisely,  the  expression  W^WqxW^J 
W'ljH'ioH'jotnust  be  large). 

A  central  result  of  our  calculation  is  that  no  population 
inversion  can  be  achieved  when  only  optical  phonons  are 
taken  into  account,  and  no  additional  quasielastic  scatter¬ 
ing  mechanism  comes  into  play.  This  stands  in  contrast  to 
previous  calculations  performed  in  the  limit  of  nonquantiz¬ 
ing  magnetic  fields4,7  and  is  probably  a  consequence  of  the 
discrete  structure  of  the  energy  spectrum.  The  usual  ex¬ 
planation  for  the  inversion  between  the  n  “0  and  n  “  1  LL 
is  that  the  transition  probability  IVj f  is  larger  than  IP'S 
(in  other  words,  after  emission  of  an  optical  phonon  the 
streaming  electron  rather  returns  into  a  state  near  v  “0,* 
which  corresponds  to  the  n  “  1  LL  in  the  present  case). 
This  argument  turns  out  to  be  correct  in  terms  of  the  expli¬ 
cit  phonon  rates.  An  inversion,  however,  is  prohibited,  be¬ 
cause  the  rates  H'of.lfo?  are  too  small  compared  to 

The  situation  changes  dramatically  when  an  elastic  pro¬ 
cess  is  also  considered  (IV*-  “  W*K).  In  Fig.  2  optical  and 
acoustic  rates  are  plotted  versus  the  electric  field  for 
two  different  magnetic  field  strengths.  The  parameters 
are  chosen  suitable  for  GaAs  (httiLO*36  meV,  m* 
“0.066mo,  D“I6  eV,  c,  "5000  m/s,  p“ 5.3  g/cm2,  c« 
“12.5,  «•“  10.9  (Ref.  1 4)1.  In  order  that  only  three  LL’s 
are  relevant,  the  magnetic  field  has  to  be  relatively  large, 
which  in  turn  requires  quite  high  electric  fields.  Figure  3 
shows  the  population  difference  p\—po  corresponding  to 
the  rates  of  Fig.  2.  It  is  clearly  seen  that  inversion  can  be 
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FIG.  2.  Transition  rates  W'jf  .W'jj!  and  W'JJ  as  a  function  of 
the  electric  field  for  fl  “2  T  (solid  line)  and  T  (dashed 
line).  AM  parameters  are  chosen  suitable  for  GaAs.  d  is  taken 
100  A,  T -30  K. 


achieved  when  and  fV$  come  into  the  same  order  of 
magnitude  as  Wg&;  more  exactly  when  IFjf  >  tV^  > 
Hence  the  mechanism  of  inversion  can  be  described  as  fol¬ 
lows:  The  system  is  pumped  by  the  transition  rate  H'oj. 
Since  IV  u  >  W2 o,  an  overpopulation  is  built  up  in  the 
n  —  I  level.  To  maintain  the  inversion  in  a  steady  state  the 
n  —0  level  has  to  be  depicted  by  the  rate  W02,  which  in 
turn  is  identical  to  the  pump.  The  inversion  is,  therefore, 
induced  by  the  selective  interplay  of  different  scattering 
processes. 

The  relative  strengths  of  the  acoustic-  and  optical- 
phonon  coupling  influence  the  regime  of  inversion:  When 
the  acoustic-phonon  interaction  increases,  the  inversion  re¬ 
gion  is  shifted  towards  higher  electric  fields  and  vice  versa. 
Consequently,  inversion  is  not  restricted  to  —  1 .5  <  £  <  2 
as  stated  in  the  literature,4  but  can  also  occur  at  somewhat 
smaller  or  larger  values  of  £.  It  can  disappear  completely 
when  the  relative  coupling  strengths  differ  too  much.  In 
addition,  the  absolute  value  of  the  magnetic  field  has  an 
influence  on  the  regime  of  inversion:  A  higher  magnetic 
field  results  in  a  shift  towards  higher  £  (see  Fig.  3). 

An  increase  of  the  electric  field  beyond  the  maximum 
value  of  pi  —po  destroys  the  inversion,  because  the  n  “1 
LL  starts  streaming  (p\  —po - l ).  At  still  higher  elec¬ 

tric  fields,  when  all  LL’s  are  in  a  streaming  state,  the  elec¬ 
trons  are  distributed  equally  between  the  LL’s  (pf  “-$•). 
The  opposite  case  of  low  electric  fields  also  leads  to  an 
equidistribution,  but  this  stems  from  the  assumption  of  to¬ 
tally  elastic  acoustic-phonon  scattering  and  has  no  physi¬ 
cal  significance.  The  region  of  very  low  and  very  high 
fields  must  be  regarded  with  caution  anyway,  because 
more  than  three  LL’s  become  relevant. 

At  this  point  we  have  to  discuss  the  validity  of  our  ap¬ 
proximations  at  different  magnetic  fields.  Our  calculation 
is  restricted  to  a  regime  where  only  three  LL’s  participate 


FIG.  3.  Population  difference  p\  —po  vs  electric  field  for  B  “2 
T  (solid  line)  and  i")T  (dashed  line).  For  B  “2  T  the  max¬ 
imum  occurs  at  4100  V/cm  (£"2.15,  e£/“2.1hw„  pi  “0.48, 
po“0.37)  and  for  B  -3  T  at  4900  V/cm  (£“2.7,  e£/  “l.4ho>*. 
pi  -0.48,  p»  “0.38). 


in  the  pumping  process.  We  have  neglected  transitions  in¬ 
volving  a  change  of  k,\  however,  these  transitions  become 
more  important  for  eEl  <  ho>e,  which  is  a  range  suffi¬ 
ciently  below  the  inversion  regime  (see  Fig.  3).  These 
transitions  tend  to  flatten  the  steep  rise  of  the  optical- 
phonon  rates,  leading  to  a  small  shift  of  the  inversion  max¬ 
imum  to  lower  electric  fields.  But  as  long  as  fVoj  >  fVffi 
the  inversion  will  probably  still  hold. 

Experiments  have  been  performed  at  rather  low  mag¬ 
netic  fields  to  fulfill  the  condition  £;S2  for  reasonable  elec¬ 
tric  fields.  In  this  case  a  large  number  of  LL’s  would  have 
to  be  considered  in  our  analysis.  This  makes  an  explicit 
evaluation  very  difficult.  On  the  other  hand,  eEl  becomes 
much  larger  than  ha>e  (if  E/B  is  kept  constant)  in  this  sit¬ 
uation.  Consequently,  the  scattering  process  considered  in 
this  work  becomes  clearly  dominant  due  to  the  increasing 
overlap  of  the  wave  functions  with  decreasing  magnetic 
field.  This  is  further  evidence  that  this  process  gives  an 
essential  contribution  to  the  electron  distribution  function 
between  LL’s. 

In  summary,  we  have  treated  the  problem  of  the  elec¬ 
tron  distribution  in  crossed  electric  and  quantizing  mag¬ 
netic  fields  with  approximations  allowing  an  explicit  solu¬ 
tion.  The  possibility  of  LL  inversion  has  been  demonstrat¬ 
ed  when  polar  optical-phonon  scattering  and  another 
quasielastic  scattering  process  are  taken  into  account. 

A  treatment  which  is  directly  comparable  to  an  experi¬ 
mental  situation,  however,  should  include  up  to  5  or  6 
LL’s.  Furthermore,  the  finite  energy  of  acoustic  phonons 
should  be  considered  and  the  solution  of  the  coupled  in¬ 
tegral  equations  ( k,  dependence)  is  necessary.  But  also, 
the  present  simple  approach  reveals  the  qualitative 
behavior  of  the  system. 

We  gratefully  acknowledge  helpful  discussions  with 
Professor  W.  Zawadzki  and  Dr.  R.  Lassnig.  This  work 
was  supported  by  the  European  Research  Office  of  the 
U.S.  Army. 
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ENERGY  RELAXATION  PHENOMENA  IN  GaAs/GaAlAs  STRUCTURES 
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Institut  fur  Exper imentalphysik ,  Univ.  Innsbruck 
A-6020  Innsbruck,  Austria 


The  energy  relaxation  of  2D  electrons  in  GaAs/GaAlAs  structures  has 
been  investigated  by  analysing  the  electric  field  dependence  of  Shubnikov- 
de  Haas  oscillations,  the  far  infrared  emission  and  photoluminescence 
spectra.  A  quite  general  behavior  of  the  electron  heating  AT  ^T^  -  T^  as 
a  function  of  the  input  power  per  electron  Pg  is  found:  T  <*  »Pe.  The  corres¬ 
ponding  energy  relaxation  times  in  the  range  of  nsec  are  independent  of 
the  electron  temperature  up  to  30  K  and  inversly  proportional  to  the  elec¬ 
tron  density.  At  higher  electron  temperatures  the  energy  relaxation  is 
governed  by  optical  phonon  emission.  However,  the  onset  depends  on  electron 
concentration  and  is  different  for  heterostructures  and  quantun  wells. 

From  intensity  dependent  cyclotron  resonance  transmission  experiments 
Landau  level  lifetimes  between  0.2  ns  and  1  ns  depending  on  the  electron 
density  are  found  in  agreement  with  data  from  time- resolved  photolunines- 
cence . 


INTRODUCTION 

The  mechanism  of  energy  relaxation  is  of  fundamental 
interest  for  the  understanding  of  electric  field  effects  in 
two-dimensional  (2D)  electron  systems.  At  low  temperatures, 
where  lattice  scattering  is  weak,  hot  electron  phenomena  can 
be  produced  by  fields  as  low  as  V/cm  in  Si  as  well  as  GaAs 
inversion  layers. 

The  first  investigations  of  energy  relaxation  in  2D  systems 
were  performed  by  Hess  et  al.1.  The  energy  loss  was  determined 
from  Shubr.ikov-de  Haas  (SdH)  experiments  as  a  function  of  elec¬ 
tric  field  for  temperatures  up  to  20  K  in  p-Si  inversion  layers. 
Similar  investigations  were  performed  by  Englert  and  Landwehr2 
for  Si-n- inversion  layers.  Sakakai  et  al.  investigated  with 
the  same  technique  n-GaAs  inversion  layers  in  heterostructures. 
For  temperatures  below  30  K  a  dominance  of  acoustic  phonon  re¬ 
laxation  is  found  for  both  materials. 

Emission  techniques  including  the  analysis  of  the  broad-- 
band  hot  electron  emission'*'^,  subband  emission®  and  the  plas- 
mon  emission7'®  have  been  used  to  determine  the  energy  loss 
rate  as  a  function  of  input  jower.  It  was  found  that  this  rate 
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is  sample-independent  as  long  as  acoustic  phonon  scattering  is 
dominant.  A  theoretical  description  of  this  phenomenon  was 
given  by  Vass9'1^  and  Pries"' 


Very  fundamental  information  on  the  hot  electron  distri¬ 
bution  and  energy  loss  rate  is  obtained  from  time  resolved. 
pHotolumincscence  experiments  performed  by  rhah  et  al.  1  . 

Over  a  wide  range  of  electric  fields  a  hot  electron  temperature¬ 
like  distribution  function  was  found.  For  temperatures  above 
50  K  the  energy  loss  rate  is  dominated  by  the  emission  of 
longitudinal  optical  (LO)  phonons.  Evidence  for  screening  of 
the  electron  LO  phonon  interaction1®'1®  and  for  hot  phonon 
effects1^'1®  is  found. 

The  energy  relaxation  in  the  presence  of  a  strong  magnetic 
field  has  been  studied  only  recently.  The  saturation  behavior 
of  the  cyclotron  resonance  (CR)  absorption  as  a  function  of 
laser  power  yields  an  electron  density  dependent  energy  relaxa¬ 
tion  time19.  From  time  dependent  luminescence  experiments  in 
GaAs/GaAlAs  quantum  wells  a  decrease  in  the  relaxation  rate 
with  magnetic  field  is  found  for  fields  below  10  T20  and  an 
increase  for  fields  of  20  T  * 1 . 

In  the  present  paper  a  summary  on  experimental  techniques 
to  determine  the  energy  relaxation  in  GaAs/GaAlAs  heterostruc- 
turcs  in  the  temperature  range  up  to  100  K  is  given.  The  results 
will  be  compared  with  recent  theoretical  results.  In  the  pre¬ 
sence  of  a  magnetic  field  a  qualitative  description  of  the 
experimental  findings  will  be  presented  due  to  the  lack  of 
theoretical  work. 


ENERGY  RELAXATION  AT  ZERO  AND  LOW  MAGNETIC  FIELDS 

The  electron  heating  and  the  energy  relaxation  in  2D 
GaAs  was  first  investigated  by  analysing  Shubnikcv-dc  Haas  (Sdll) 
oscillations  at  low  magnetic  fields  (B  <  3  T)  as  a  function  of 
the  input  power.  Optical  techniques  as  photoluminescence  and 
far  infrared  (FIR)  emission  at  zero  magnetic  field  followed. 

The  classical  technique  evaluating  the  time  dependence  of  the 
hot  electron  current  following  a  small  step  variation  in  elec¬ 
tric  field22  has  not  been  applied  to  a  2D  system  so  far. 

a)  Electric  field  dependence  of  SdH  oscillations 

A  damping  of  SdH  oscillations  in  magneto-resistance  is 
observed  when  the  2D  electron  gas  (2DEG)  is  heated  in  an  elec¬ 
tric  field.  To  determine  the  electron  temperature  T  from  the 
temperature  and  electric  field  dependent  oscillations,  the 
Dingle  temperature  must  be  constant23.  If  that  is  the  case  the 
change  in  amplitude  can  be  described  by  one  parameter  Tg.  Fig.  1 
shows  a  typical  result  of  SdH  measurements  after  Sakaki  et  al.3 
at  4.2  K.  The  decay  in  SdH  amplitude  is  clearly  observed  with 
increasing  current  (corresponding  to  increasing  electric  field). 
Significant  differences  in  the  heating  behavior  are  found  as  a 
function  of  zero  field  mobility.  However,  a  very  general  be¬ 
havior  is  obtained  when  Te  is  plotted  against  the  input  power 
per  electron.  The  weakness  of  this  technique  lies  in  the  fact 
that  Pe  =  euE2  is  calculated  using  the  electric  field  dependent 
mobility  at  zero  magnetic  field  while  the  Te  values  are  deter¬ 
mined  from  SdH  oscillations  at  finite  B. 
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Fig.  1:  Shubnikov-de  Haas  oscil¬ 
lations  of  Sample  K-6  after  Ref.  3. 
The  sample  has  a  mobility  of 
2.05  x  105cm2/Vs. 
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Fig.  2:  Photoluminescence  spectra 
of  GaAs/AlGaAs  heterostructures 
for  three  in-plane  electric  fields. 
The  electron  temperatures  deduced 
from  the  high-energy  slope  are  indi¬ 
cated.  Relative  intensities  of  curves 
are  arbitrary  (after  Ref.  13). 


The  results  from  Sakaki  et  al .  will  be  shown  in  Fig.  4 
together  with  the  data  from  the  other  techniques.  Similar 
results  were  obtained  by  other  authors  studying  the  mobility 
and  electron  temperature  as  a  function  cf  electric  field34'“'' 


b)  Hot  electron  luminescence  analysis 


Luminescence  experiments  on  GaAs/GaAlAs  multiple  quantum 
well  (MQW)  structures  were  performed  by  Shah  et  al.13,T  .  A 
cw  infrared  dye  laser  (7700  k  '  X  <  7900  k,  P  =  2  mW)  was  used 
to  weakly  excite  carriers  in  the  GaAs  across  the  gap.  Lumines¬ 
cence  spectra  as  a  function  of  electric  field  were  analysed 
with  a  double  monochromator.  Typical  spectra  after  Shahetal.14 
are  shown  in  Fig.  2  for  three  different  electric  fields.  The 
high  energy  tail  is  attributed  to  electron  hole  recombination. 
This  tail  can  be  well  fitted  by  a  single  exponential  (see  Fig. 2) 
at  all  fields  up  to  the  highest  field  (150  V/cm) .  The  spectra 
are  fitted  with  a  Maxwell-Boltzmann  type  distribution  with  a 
temperature  TQ  higher  than  T^. 

The  obtained  electron  heating  AT  is  also  plotted  in  Fig. 4 
over  the  input  power  derived  from  I-V  measurements.  In 
a  steady  state  situation  the  power  input  must  be  equal  to  the 
power  loss  to  the  lattice. This  technique  gives  direct  means 
to  determine  energy  loss  rates  at  zero  magnetic  field.  It  is- 
found  that  the  electron  heating  shews  a  somewhat  weaker  depen¬ 
dence  cn  input  power  than  the  previous  technique  and  that  there 
is  also  practically  no  dependence  on  the  carrier  concentration. 
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The  disadvantage  of  this  technique  is  the  simultaneous 
excitation  of  holes  which  can  influence  the  total  energy  loss 
rate.  In  addition  we  have  to  be  aware  of  the  fact,  that  the 
experiments  were  performed  on  MQW  which  may  give  a  somewhat 
different  behavior  than  heterostructures.  A  comparison  will 
only  be  meaningful  for  rather  wide  wells  (as  in  the  present 
case:  well  thickness  250  A) . 


c)  Far  infrared  (FIR)  emission 

The  basic  idea  of  this  experiment  is  to  heat  up  the  2D 
carrier  gas  by  electric  field  pulses  and  to  measure  the  spec¬ 
trum  and  intensity  of  the  emitted  broadband  FIR  light  and  to 
correlate  it  with  an  electron  temperature. 


FIR  broadband  emission  has  first  been  used  in  Si-MOSFBTs 
to  determine  the  energy  loss  of  2D  carriers’.  The  electron 
temperature  Tg  was  determined  from  the  absolute  value  of  the 
emitted  power.  For  2D  electrons  in  GaAs  a  somewhat  different 
experimental  technique  based  on  a  relative  intensity  measure¬ 
ment  was  applied.  The  signals  of  two  narrowband  FIR  detectors 
(n-GaAs  at  35  era'1,  Ga  doped  Ge  at  100  cm-1)  are  measured  as 
a  function  of  the  applied  electric  field  and  calculated  for 
each  detector  according  to 

V  -  F  R  (w)  A  (a;)  IgB  (uj  ,Te )  "do  (1) 

where  R(oj)  =  R  r(uj)  and  A(u >)  =  A  a  (<■>)  are  the  frequency  depen¬ 
dent  and  knowndetector  response0  and  the  absorptivity  of  the 
electron  system  respectively.  In  equilibrium  the  spectral 
emission  intensity  of  a  system  with  absorptivity  A(oj)  is  given 
by 

I  (w,T)  =  Ibb(w,T)  •  A ( oj)  (2) 


with  I0B  the  black  body  emission  intensity  (Planck-function) . 
Assuming  a  quasi  equilibrium  of  the  electron  system  at  a  tempe¬ 
rature  Te,  we  can  express  the  broadband  emission  from  the  2D 
carrier  gas  by  I(<,),Te). 


The  absorptivity  can  be  expressed  for  normal  incidence 
according  to 

4ReF 

A  (  id)  =  - 5  (3) 

(/f;  ♦  i  ♦  Fr 

with  F  =  o(u))/e  *c  and  e  being  the  dielectric  constant  of  the 
substrate.  For  °o(ci))  a  Drude  type  frequency  behavior  is  assumed. 


The  unknown  parameters  in  equ.  (1)  are  determined  by 
measuring  with  two  detectors  at  two  different  electric  fields 
giving  four  equations  for  four  unknown  factors.  Fig.  3  shows 
the  observed  detector  signal  in  Volts  as  a  function  of  elec¬ 
tric  field  for  two  different  samples:  a  low  mobility  sample 
with  ns  =  8.7  x  lO^cm'2  and  u  =  1.6  x  lO^cmVVs  (Sample  1)  and 
a  high  mobility  sample  with  ns  =  2.4x1011cm-2  and  u  =  1.0  x 
10®cm2/Vs  (Sample  2,  dark).  It  is  directly  evident  that  the 
lower  mobility  sample  shows  a  considerably  higher  emission 
signal . 
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Fig.  3:  Detector  signals  as  a 
function  of  applied  electric 
field  for  two  different  detec¬ 
tors  after  Ref.  5. 

The  sample  and  detector  proper¬ 
ties  are  given  in  the  text. 


Fig .  4 :  Electron  heating  versus  input 
power  per  electron  obtained  from  the 
different  techniques: 

Ref.  3:  (A)  sample  described  in  Fig.  1, 
(O)  Ns*3.5xl0ncm-a,  n*7.0xl04cm2/Vs. 
Ref.  5:  (x)  Sample  1,  (o)  Sample  2  (see 
text);  dashed  curve:  data  from  Ref.  13. 


Quantitative  values  for  the  electron  temperature  are  de¬ 
termined  from  the  fit  of  Fig.  3  according  to  equ.  (1)  and  also 
plotted  as  a  function  of  the  input  power  in  Fig.  4.  The  basic 
behavior  of  AT  is  the  same  as  determined  with  the  other  tech¬ 
niques.  The  emission  data  have  a  somewhat  steeper  slope  than 
the  luminescence  data.  There  is  also  a  weak  dependence  on 
density  evident  in  the  emission  data. 

The  emission  technique  seems  to  be  the  most  appropriate 
method  to  determine  T„.  However,  the  analysis  requires  several 
assumptions  which  might  introduce  sample  dependences. 

A  similar  FIR  emission  technique  is  the  analysis  of  the 
emission  from  2D  plasmon  excitations  as  demonstrated  first  in 
Si2/,/  and  later  in  GaAs28.  In  a  recent  paper  Sambe  et  al.29 
have  analyzed  the  emitted  power  from  2D  plasmon  excitation  in 
GaAs  and  derived  electron  temperatures  as  a  function  of  elec¬ 
tric  field.  As  mobility  on  electric  field  data  were  not 
published  an  inclusion  of  their  data  in  Fig.  4  is  not  possible. 

The  summary  of  the  experimental  findings  is  shown  in 
Fig.  4  as  a  plot  of  AT  =  T  -  TL  over  the  input  power  per  elec¬ 
tron  Pe.  It  is  clearly  evident  that  the  emission  data  continue 
very  well  the  SdH  results.  The  data  can  be  described  by  a  re¬ 
lation  AT  «  /euE2  up  to  electron  temperatures  of  30  K.  For 
higher  temperatures  the  emission  data  show  a  weak, while  the 
luminescence  data  show  a  considerably  stronger  change  in  slope. 
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<D  Determination  of  the  energy  relaxation  time 


To  derive  an  energy  relaxation  time  from  the  data  an 
energy  balance  equation  for  the  average  electron  energy  e (T) 
using  Fermi-Dirac  statistics  with  an  electron  temperature  Tc 
is  used: 

2  2 

|c(Te)  -e(TL)|/te  =  — 6|p—  (Tg  -  T2 )  /t£=  euE2  =  Pg  (4) 

where  e„  is  the  Fermi  energy  and  u  the  electric  field  depen¬ 
dent  mobility.  It  is  directly  evident  that  the  expression 
AT  «  Tg  -  =  const.  /Pe ■ Ep • t £  describes  the  experiments 

quite  well  indicating  that  for  a  given  Pe  the  product  cp‘T 
has  to  be  sample  independent.  As  Ns  =  D(c)*cp  with  D(c)  =  const 
t£  is  inversely  proportional  to  the  2D  carrier  density. 

The  resulting  energy  relaxation  times  as  a  function  of 
density  as  obtained  from  Fig.  4  for  the  range  of  electron 
temperatures,  where  the  slope  is  constant,  is  plotted  in  Fig. 5. 
As  long  as  the  slope  of  AT  over  Pe  is  constant  we  obtain  a 
Te-value  independent  of  Te  for  a  given  sample.  It  is  clearly 
evident  that  the  energy  relaxation  time  is  a  linear  inverse 
function  of  the  density  as  a  consequence  of  the  analysis  with 
equ.  (4) .  Both  methods,  the  SdH  oscillations  and  the  emission 
technique,  give  within  experimental  accuracy  the  same  rc-values 
so  that  we  can  be  quite  confident  about  the  obtained  results. 


Fig.  5:  Calculated  energy  relaxation 
time  versus  2D  electron  density  for 
three  samples  from  Ref.  3  (a,  •  ,  A) 
and  two  samples  from  Ref.  5  (o,  x)  . 

The  symbols  represent  the  same  samples 
as  in  Fig.  4. 


Fig.  6:  Energy  relaxaticn  times 
versus  electron  temperature  for 
Sample  1  (x)  and  Sample  2  (o) 
from  Ref.  5  and  for  a  sample  with 
Ng  =  3.9 x 101 ’em'2  and  u  =  7.9 
x  104cm2/Vs  from  Ref.  3.  The  full 
curve  represents  the  loss  due  to 
optical  phonons14. 


An  evaluation  of  x£  over  Te  is  only  meaningful  for  tempe¬ 
ratures  above  30  K,  where  the  slope  starts  changing.  Eut  also 
in  the  higher  temperature  range  we  use  equ.  (4)  to  evaluate  xE 
which  means  that  we  use  degenerate  Fermi-Dirac  statistics 

with  an  electron  temperature  Te.  The  use  of  the  statistics  has 
a  strong  influence  on  the  results.  A  plot  of  xE-values  from 
the  emission  data  for  two  different  samples  and  for  the  photo- 
luminucoccncc  data  is  shown  in  Fig.  6  as  a  function  of  the 
inverse  electron  temperature  (and  electron  temperature) .  The 
energy  relaxation  times  derived  from  the  emission  data  show 
a  similar  behavior  for  both  concentrations:  For  temperatures 
above  40  K  te  decreases  significantly  reaching  values  of  a 
few  hundred  psec.  The  analysis  of  the  photoluminescence  data 
gives  a  considerably  stronger  dependence  of  tc  on  Te  which  is 
already  evident  in  the  weaker  AT  over  Pe  curve.  The  full 
curve  shows  calculations  of  the  energy  relaxation  rate  due  to 
optical  phonons  for  a  3D  case'^.  This  curve  has  a  slope  pro¬ 
portional  to  the  optical  phonon  energy  (36  moV)  in  GuAu.  It 
is  evident  from  Fig.  6  that  for  both  techniques  there  is  a 
tendency  toward  the  optical  phonon  line.  However,  there  is  a 
clear  difference  in  the  behavior  of  the  results  from  the  two 
techniques.  The  photoluminescence  data  seem  to  be  more  strongly 
influenced  by  the  optical  phonon  emission.  The  reason  for  this 
is  not  clear  but  might  be  due  to  the  fact  that  in  the  photo¬ 
luminescence  experiment  electrons  are  always  excited  above 
the  optical  phonon  energy  by  the  laser  excitation. 

In  Ref.  14  it  was  argued  that  the  energy  loss  rate  for 
electrons  was  considerably  smaller  than  for  holes.  However, 
we  do  not  find  this  result  since  the  energy  loss  rates  shown 
in  Fig.  6  are  comparable  with  loss  rates  of  holes.  The  differ¬ 
ence  in  the  data  comes  from  the  analysis.  Shah  et  al.14  used 
non-degenerate  statistics.  The  xE-values  obtained  this  way  are 
longer  by  a  factor  Cp/ AT.  For  samples  with  Fermi  energies  of 
the  order  of  20  meV  the  use  of  non-degenerate  statistics  seems 
not  appropriate. 

A  critical  analysis  of  the  energy  relaxation  was  performed 
by  Tsubaki  et  al.  who  analysed  transport  studies  as  a  function 
of  electric  field  from  several  authors.  The  influence  of  the 
used  statistics  in  the  evaluation  was  examined.  It  is  shown 
that  the  application  of  degenerate  statistics  is  meaningful 
down  to  ep/kBT-values  of  3.  However,  the  obtained  ie-values  are 
independent  of  Ng  and  are  nearly  an  order  of  magnitude  shorter 
than  the  data  shown  in  Fig.  6.  At  low  temperatures  the  values 
are  comparable,  however,  the  slope  of  the  xe  versus  Te  plot  is 
considerably  steeper  (in  Ref.  30)  than  the  optical  phonon  curve 
(see  full  curve  in  Fig.  6). 

A  theoretical  description  of  the  energy  relaxation  rate 
for  2D  electrons  in  GaAs  in  the  case  of  acoustic  phonon 
scattering  has  been  published  by  Vass9.  Vass  derived  an  ex¬ 
pression  for  the  energy  loss  rate:  AT  =  const.  /Pg  where  the 
calculated  value  of  the  constant  agrees  well  with  the  experi¬ 
mental  findings  in  GaAs  and  Si4.  The  power  loss  i&  dominated 
by  acoustic  deformation  potential  scattering  up  to  tempera¬ 
tures  of  40  K.  In  a  recent  paper1 0  values  for  the  energy  re¬ 
laxation  time  as  a  function  of  electric  field  and  electron 
concentration  were  reported  including  optical  phonon  emission 
and  electron-electron  scattering.  The  calculated  values  show 
the  same  electric  field  (or  Te)  dependence  as  the  experimental 
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values  but  are  an  order  of  magnitude  to  high.  The  predicted 
carrier  density  dependence  is  also  too  weak. 

Price  has  investigated  the  mobility  and  energy  loss  rate 
for  piezo-elcctric  and  deformation  potential  acoustic  phonon 
scattering^.  For  temperatures  up  to  10  K  he  predicts  piezo¬ 
electric  scattering  to  be  dominant  giving  a  considerably 
steeper  dependence  of  AT  on  Pe  than  from  equ.  (4).  For  higher 
temperatures  he  predicts  a  linear  dependence  of  AT  on  Pe  which 
is  also  dependent  on  the  sample  concentration.  These  results 
are  not  in  agreement  with  the  experimental  observations. 


ENERGY  RELAXATION  IN  HIGH  MAGNETIC  FIELDS 

A  direct  way  to  obtain  the  energy  relaxation  rate  in  the 
presence  of  a  strong  magnetic  field  is  the  measurement  of  the 
incoherent  saturation  of  the  CR  transmission.  This  technique 
has  been  applied  to  bulk  n-type  InSb  and  GaAs2^*22.  it  was 
found  that  an  interplay  between  electron-electron  scattering 
and  optical  phonon  emission  governs  the  energy  relaxation. 

The  same  technique  was  applied  to  GaAs/GaAlAs  heterostructures 
by  Helm  et  al.1*.  A  high  power  cw  optically  pumped  FIR  laser 
was  used  to  perform  the  experiments.  Typical  transmission 
spectra  for  two  different  samples  are  shown  in  Fig.  7.  It  is 
evident  that  powers  below  1  W/cm2  reduce  the  transmission  al¬ 
ready  significantly  and  that  the  lower  density  sample  saturates 
stronger.  In  the  analysis  rate  equations  for  three  Landau 
levels  and  a  constant  relaxation  time  Td  was  used.  Fig.  8 
shows  the  obtained  relaxation  time  tb  plotted  against  the 
total  2D  carrier  concentration.  A  systematic  decrease  of  xB 
with  Ns  is  observed  even  though  there  is  considerable  scatter 
in  the  data.  The  dependence  of  tB  on  Ns  indicates  that  the 
mechanism  may  be  similar  to  that  found  in  bulk  semiconductors^ 1 . 
It  is  also  interesting  to  note  that  in  the  zero  magnetic  field 
case  a  similar  dependence  is  found. 

Data  on  relaxation  rates  between  Landau  levels  were  also 
reported  by  Ryan  et  al.  0  and  Hollering  et  al.21  using  psec 
photoluminescence  techniques.  Ryan  et  al .  find  relaxation 
times  in  the  order  of  0.3  ns  between  the  first  excited  and 
lowest  Landau  level  for  MQW  samples  with  Ns  =  5x 10''cm"2. 

The  relaxation  times  in  magnetic  fields  of  7  and  8  T  are  longer 
than  without  magnetic  field.  On  the  other  hand  Hollering  et  al. 
find  a  considerably  shorter  time  for  fields  of  20  T  as  compared 
to  the  zero  magnetic  field  case. 

The  whole  behavior  seems  to  be  consistent:  In  fields 
below  10  T  the  energy  relaxation  time  derived  is  increased 
as  compared  to  the  field  free  case.  The  data  from  Fig.  8  have 
to  be  compared  with  tG-values  in  Fig.  6  at  the  highest  tempera¬ 
tures.  This  comparison  shows  that  they  are  somewhat  longer  in 
agreement  with  Ryan's  observation  on  the  same  sample.  The  con¬ 
siderably  shorter  values  for  20  T  are  somewhat  surprising.  In 
Ref.  21  they  are  briefly  explained  by  a  reduced  screening  of 
the  electron  LO  phonon  interaction. 

Further  experimental  and  theoretical  studies  in  the  magne¬ 
tic  field  are  necessary  to  get  clear  evidence  for  the  relaxation 
mechanism.  In  a  very  recent  paper  Rodriguez  et  al.^3  report  an 
intensity  dependent  CR  transmission  experiment  at  rather  low 
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Fig.  7:  cylotron  resonance  absorption 
for  two  different  GaAs/GaAlAs  hetero¬ 
layers  at  118  urn  for  different  laser 
intensities:  left  sample  with  Ns  * 

1.2  x  lO^cm-^,  right  sample  with 
Ns  *  3.0  x  10* 
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Fig.  B:  Calculated  lifetimes 
from  saturation  behavior  as  a 
function  of  Ng:  (o)  cyclotron 
r e  sonan  ce ,  ( A )  pho  to 1 umi ne scence . 
The  dashed  curve  represent  a 
linear  dependence. 


frequencies.  Saturation  is  observed  at  considerably  higher 
intensities  than  in  Ref.  19  for  magnetic  fields  below  4  T. 

This  indicates  a  non-resonant  heating  process  similar  to  direct 
electric  field  heating.  The  observed  relaxation  times  in  the 
order  of  10  psec  probably  directly  reveal  the  optical  phonon 
relaxation  time.  The  main  reason  for  this  behavior  is  the 
nearly  equidistant  Landau  level  spacing  at  low  magnetic  field. 
In  this  situation  a  resonant  saturation  of  a  transition  bet¬ 
ween  two  levels  is  not  possible.  Only  if  a  sufficient  amount 
of  nonparabol icity  is  present  the  level  shifts  are  larger 
than  the  individual  Landau  level  linewidths  and  a  saturation 
of  an  individual  transition  becomes  possible  as  demonstrated 
in  Ref.  19. 
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ABSTRACT 

The  spec  true  of  the  stimulated  far  Infrared  amission  from  p-Ceraanlua  In  crossed  electric 
and  magnetic  fields  Is  studied  by  means  of  a  tunable  narrowband  GeAs-detector.  A  multimode 
spectrum  Is  observed  from  alrrorless  samples  and  quantitatively  explained  In  terms  of 
waveguide- like  modes. 
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nmeooucriGM 

During  the  last  three  years  It  has  bacomo  possible  to  realize  high  power  far  Infrared  (FIR) 
semiconductor  lasers  /1-7/.  The  stimulated  amission  is  caused  by  a  population  Inversion 
between  the  light  and  heavy  hole  bands  In  lightly  p -doped  Geranium  crystals  subjected  to 
crossed  electric  and  magnetic  fields. This  population  Inversion  builds  up.  when  the  heavy 
holes  interact  strongly  with  the  optical  phonons  ("streaming  notion"),  whereas  the  light 
holes  accumulate  In  momentum  space  /8A  Up  to  now.  lasing  has  been  observed  for  electric  and 

aagnetlc  fields  satisfying  the  relation  1  S  S  1.9.  where  Cj,»  v®iT^7ag.(B/E)  »  vjjpB/E. 

with  the  optical  phonon  energy  hw _ »  37  aeV  and  the  effective  ansa  of  the  heavy  holes 

op 

eg  a  0.35  ap.  Emission  wavelengths  between  75  im  and  250  pm  have  been  reported  with  output 

powers  up  to  10  Matts.  The  radiation  has  been  observed  both  from  alrrorless  saaples  which 
form  a  waveguide- like  cavity  and  from  samples  with  external  mirrors  /7/. 

A  point  not  yet  clarified  Is  the  node  structure  of  the  emitted  light.  Different  authors 
reported  a  variety  of  different  spectra  /4-7/,  and  no  unambiguous  explanation  could  be 
given,  partly  due  to  Insufficient  resolution  of  the  spectrometers  used  In  these  experiments. 

We  have  studied  the  spectrum  of  the  stimulated  emission  from  p-Ce  by  means  of  an  extresmly 
narrowband,  aagnetlc  field  tunable  n-CaAs  detector.  By  analyzing  samples  with  different 
dimensions,  we  are  able  to  give  a  quantitative  explanation  of  the  node  structure  In  the 
configuration  without  external  cavity.  Recent  results  with  external  alrrors  are  also  briefly 
discussed. 


EXPERIMENT 

Three  different  saaples  were  used  for  the  experiments,  all  cut  from  a  Ga-doped  Ce-crystal 
with  pWI^-Njj-1.2'1014  cm*3.  The  saaple  dimensions  are  given  in  table  1: 


TAH1F  1  3—iM  dimensions 


Sample 

d 

A 

1 

1 

3.7 

2.4 

46.9 

2 

4.2 

3.15 

43.1 

3 

5.2 

3.75 

4S.4 

The  distance  between  the  electrical  contacts  It  d.  1  is  the  sample  length  In  the  direction 
of  the  aagnetlc  field  (In  our  case  the  (lll)-dlrectlon).  and  a  Is  the  width.  Electrical 


contacts  were  weile  by  evaporating  Indium  (with  Maw  cold)  and  al loy  ing  at  400“C.  All  faces 
of  the  saaples  wore  polished  sad  are  parallel  within  X".  which  turned  out  to  be  crucial 
for  the  achlevaaent  of  laser  oscillation.  Voltage  pulses  of  up  to  1000  V  with  lps  duration 
were  applied  to  the  saaples  with  a  low  l^edance  high  power  pulse  generator.  The  typical 

current  through  the  saaples  was  200  A.  in  accordance  with  the  relation  J^Jpevj^.  where  J  is 

the  current  density,  and  Kvh  Is  the  approximate  drift  velocity  of  the  heavy  holes  under 

op 

"streaming  motion"  conditions  /8/.  The  pulse  repetition  rate  was  chosen  to  be  around  SHz. 
For  rates  higher  than  10Hz  lasing  braaks  down  due  to  the  saeg)le  heating.  The  experiment! 
were  performed  at  4.2K.  so  the  saaple  was  loeersed  In  liquid  He  at  the  center  of  a 
superconducting  solenoid.  To  check  the  range  of  electric  and  angnetlc  fields  where  lasing 
can  be  achieved,  the  radiation  was  detected  with  a  broad  band  Ce:Ca-detector  with  peak 
sensitivity  at  about  110  pa.  A  magnetic  field  tunable  narrowband  GaAs-detector  was  used  for 
the  detailed  analysis  of  the  emission  spec trua .  It  was  placed  In  the  center  of  a  second 
Mgnet  In  the  cryostat  below  the  saaple.  The  Zeeman-spilt  ls-2p+  shallow  donor  line  Is  used 
for  detection  /9/.  This  line  is  tunable  between  45  cm" ‘  and  IX  cm"1  for  angnetlc  fields  of 
1  T  to  €.5  T.  Due  to  the  fact  that  the  CaAs  contains  only  one  chemical  species  of  shallow 
donors,  the  absorption  line  Is  extremely  narrow  (0.25  ca~().  and  the  CaAs  can  be  used  as  a 
high  resolution  spectrometer. 


RESULTS  AND  DISCUSSION 

Figure  1  shows  typical  traces  of  stimulated  amission  from  saaple  2.  recorded  with  a 
Ge:Ca-detector.  The  detector  signal  due  to  stimulated  emission  Is  a  factor  of  10a-103  above 
the  level  of  spontaneous  amission.  The  region  of  lasing  extends  from  somewhat  less  than 
Bal  T  for  the  lowest  electric  field  of  E  a  1000  V/ca  up  to  more  than  2  T  for  E  =t  2000  V/ca. 
This  corresponds  to  1.3  (  [  (  2.4.  Ve  observe  lasing  at  scwaewhat  higher  angnetlc  fields, 
that  is  for  larger  values  of  f.  and  also  In  a  wider  range  of  C  than  has  been  reported 
previously  /1-5/.  A  possible  explanation  of  this  discrepancy  is  the  effective  lnhoaogenelty 
of  the  (agnatic  field  In  the  saaple.  However,  this  Is  only  about  10X  over  the  whole  saaple 
length.  Another  reaaon  could  be  the  crystallographic  direction  of  the  electric  field,  which 
Is  not  known  In  our  experiment.  Comparing  the  magnitude  of  the  signal  with  the  signal 
obtained  from  a  FIR  gaa  laser  and  a  InSb-Landau  source,  the  emitted  power  was  estimated  to 
reach  100  a I. 


B 

Fig.  1:  Ge:Ca-detector  signal  from  sample  2  as  a  function  of  nmgnetic 
field  for  different  electric  field  values. 


To  analyze  the  spectrum  of  the  stimulated  emission  the  electric  and  magnetic  fields  applied 
to  the  sasple  were  set  to  a  fixed  value,  while  the  aagnetlc  field  of  the  GaAs-detector  was 
tuned.  Using  the  ls-2pr  shallow  donor  transition,  the  frequency  of  the  detected  radiation  Is 
linearly  related  to  the  magnetic  field  over  a  wide  range.  Figure  2  shows  three  spectra  from 
saaple  2.  for  three  different  values  of  the  electric  and  magnetic  field.  The  spectra  consist 
of  10  -  X  Unas  with  a  regular  node  spacing  of  about  1.4  cm'1.  Their  positions  remain  fixed 
with  changes  in  the  applied  fields.  In  the  hlgh-galn  region  some  additional  lines  appear, 
disturbing  the  regular  structure.  The  measured  widths  are  about  0.4  cm*'.  Ve  conclude  that 
the  real  line  width  of  the  emission  is  of  the  same  order  of  aagnitude  as  the  detector 
Unewldth  (0.25  cm"’).  By  varying  the  electric  and  magnetic  fields  it  Is  possible  to  obtain 
lasing  between  75  pm  and  110  pa.  There  Is  a  general  trend  towards  emission  at  shorter 
wavelengths  for  both  higher  E-  and  B-fields  (see  Fig. 2).  In  contrast  to  Andronov’s  /6/  and 
Koaiyaaa’s  /5/  work,  no  long-wave length- lasing  (150  -  250  pa)  could  be  found,  not  even  at 
lower  magnetic  fields. 


VX(errr) 


Fig.  2:  Spectrum  of  stimulated  Mission  froa  saaple  2  for  dlfforont 
electric  and  magnetic  f is Ids. 


Andronov  ot  al.  /I/  and  Koalynaa  at  al.  /2/  Suggested  that  the  lasing  modes  are  due  to  total 
reflections  at  the  saapla  surface,  since  without  mirrors  the  quality  factor  of  the  cavity  is 
far  to  low  to  allow  longitudinal  or  transverse  nodes  to  oscillate  (reflectivity  of  the 
surface  at  no  real  incidence:  36X) .  In  the  following  we  want  to  clear  up  the  origin  of  these 
nodes  and  their  spacing.  In  order  to  got  more  information,  wa  have  analyzed  two  other 
saaples  with  different  dlaansions  (see  Table  1).  The  characterization  was  performed  in  the 
sane  nay  as  described  above.  As  expected,  a  wider  node  spacing  is  observed  for  the  thinner 
saaple  and  a  smaller  spacing  for  the  thicker  saapla;  As  «  1.7cm'1 .  1.4cm'*.  and  l.lca'1 
for  saaples  1.2  and  3.  respectively.  Coopering  these  values  with  the  dimensions  a  and  d  of 
the  saaples.  it  appears  to  be  very  likely  that  the  width  a  is  the  crucial  parameter 
determining  the  node  spacing.  This  is  also  consistent  with  the  physical  reason  that  the 
surfaces  coated  with  the  electrical  contacts  are  not  serving  as  alrrors  as  well  as  the  other 
surfaces. 
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Fig.  3:  Sketch  of  the  ray  path  in  the  Cermanl un-crystal . 

For  the  calculation  of  the  node  spacing  we  consider  an  electromagnetic  wave  in  the  germanium 
crystal  propagating  with  components  of  the  wavevector  in  x-  and  z -directions,  according  to 

Fig. 3:  k  m  kx?x+  kjSj.  where  «x  and  e^  are  the  respective  unit  vectors.  The  transverse 
resonance  condition  is  given  by  /10/ 

k^n  a  k  n  sin  9  a  Hw/a  .  (1) 

where  n  is  the  refractive  index.  M  is  an  integer,  and  9  is  indicated  in  Fig.  3.  In  addition, 
we  have  the  geometrical  condition 

i  tan  9  a  Na/1  (2) 

(N  integer),  which  ensures  that  the  path  of  the  ray  is  closed  within  the  saaple.  This  leads 
to  a  node  spacing  (in  wavenumbers) 

An  e  (2a  n  sin  9)'*.  (3) 

9  must  satisfy  the  condition  H.5°  <  9  <  75.5°.  because  14.5°  is  the  angle  of  total 
reflection  for  Germanium.  One  can  assume  that  modes  with  the  smallest  possible  9  have  the 
lowest  losses,  because  they  undergo  the  fewest  reflections.  Thus,  for  each  sample  the  node 


spacing  la  determined  first  by  choosing  tha  angle.  0.  with  tha  smallest  possible  N 
(Eq.  (2)).  and  then  this  angle  is  used  in  Eq.  (3).  This  procedure  gives  a  aode  spacing.  4o. 
of  1.77ca*‘  (Nm6).  1.41ca"r  (N«5).  and  1.13  ea*1  (N»4).  for  saaple  1.2  and  3.  respectively, 

in  good  agreeaent  with  the  aaasured  values.  Tha  additional  nodes  appearing  at  the  high-gain 
regions  can  be  understood  as  nodes  with  a  larger  value  of  6.  According  to  Eq.  (3).  their 
spacing  is  saaller  than  the  fundaaental  spacing,  which  wakes  the  node  structure  in  the 
high-gain  region  rather  coaplex.  The  longitudinal  resonance  condition  analogous  to  Eq.  (1) 
is  not  considered,  because  it  would  yield  a  very  narrow  node  structure  which  is  not  observed 
in  the  experiment . 

Froa  the  above  considerations  it  appears  reasonable  to  produce  a  very  thin  saaple  in  order 
to  obtain  a  single  aode  oscillation  /ft/ .  We  prepared  a  saaple  with  awl  ns.  expecting 
oscillation  on  two  lines.  However,  it  was  not  possible  to  achieve  lasing  froa  this  sample. 
Probably  the  losses  are  too  high  because  of  the  high  order  of  the  nodes  necessary  for  this 
geometry  (N>12). 

Ve  have  also  investigated  saaples  with  an  external  sealcanfocal  resonator.  Both  airrors 
consist  of  pure  Geraanlua  plates.  Tha  curved  airror  la  coated  with  gold,  the  plane  mirror 
with  a  wash  or  a  Cr-flla  to  yield  90®  -  99X  reflection.  When  tha  reflectivity  Is  about  90X. 
no  lasing  can  be  observed.  This  is  a  further  proof  for  the  waveguide- like  modes  discussed 
above:  They  are  destroyed  by  attaching  the  airrors  on  both  ends.  On  Increasing  the 
reflectivity  to  99X  lasing  occurs  again,  aoat  likely  on  longitudinal  modes  in  this  case.  The 
analysis  of  the  aalssion  spectrum  will  be  discussed  elsewhere  /ll/. 

In  suanary,  we  have  investigated  tha  aalssion  spectrum  of  p-Ce  lasers  with  different 
dimensions.  The  aode  structure  of  saaples  without  mirrors  could  be  quantitatively  explained 
in  terms  of  waveguide-like  modes. 
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The  mode  structure  of p-Ge  far-infrared  lasers  with  and  without  mirrors  is  investigated 
Without  external  mirrors  a  multimode  spectrum  is  observed  and  quantitatively  explained  in 
terms  of  waveguidelike  modes.  An  external  resonator  drastically  reduces  the  number  of 
oscillating  modes.  For  the  first  time  single  line  operation  is  demonstrated  in  this  configuration. 


Recently,  it  has  become  possible  to  realize  high-power 
far-infrared  (FIR)  semiconductor  lasers  from p-Ge. '**  The 
stimulated  emission  is  due  to  transitions  between  the  light 
and  the  heavy  hole  bands,  which  can  form  an  inverted  sys¬ 
tem  when  the  crystal  is  subject  to  crossed  electric  and  mag¬ 
netic  fields.  Emission  wavelengths  between  75  and  250  pm 
have  been  reported  with  output  powers  up  to  10  W. Sever¬ 
al  attempts  have  been  made  to  analyze  the  spectrum  of  the 
emitted  radiation  and  to  understand  the  mode  structure,  but 
no  unambiguous  interpretation  of  the  various  spectra  ob¬ 
served  could  be  given  so  far.4'1 

We  have  studied  the  spectrum  of  the  stimulated  emis¬ 
sion  from  p-Ge  by  means  of  an  extremely  narrowband,  mag¬ 
netic  field  tunable  n-GaAs  detector.  By  investigating  sam¬ 
ples  with  different  dimensions,  we  are  able  to  give  a 
quantitative  explanation  of  the  mode  structure  in  the  config¬ 
uration  without  external  cavity.  After  attaching  mirrors  we 
observe,  for  the  first  time  in  this  configuration,  a  single  line 
oscillation. 

For  the  experiments  three  different  samples  were  used, 
all  cut  from  a  Ga-doped  germanium  crystal  with 
p  =  ,V«  —  Nd  =  1.2X  10l4cm-\  The  dimensions  are  given 
in  Table  I.  Electrical  contacts  were  made  by  evaporating 
indium  (with  5%  gold)  and  alloying  at  400  *C.  All  faces  of 
the  sample  were  polished  and  are  parallel  within  30  arcse- 
conds.  To  form  a  semiconfocal  resonator,7  two  external  mir¬ 
rors  made  of  pure  Ge  are  used.  The  plane  mirror  consists  of  a 
0.9-mm-thick  parallel  Ge  plate  coated  with  a  1500- A  chro¬ 
mium  layer.  The  reflectivity  was  determined  with  a  FIR  gas 
laser  to  be  99%.  The  curved  mirror  is  produced  by  deposit¬ 
ing  2000  A  gold  on  the  curved  face  of  a  plane-convex  Ge  lens 
(thickness  1.1  mm)  to  serve  as  a  100%  reflecting  mirror. 
The  radius  of  curvature  is  100  mm,  which  provides  a  stable 
cavity  configuration.  Both  mirrors  were  attached  to  the  sam¬ 
ple  with  silicone  oil  and  fixed  with  nail  enamel. 

Voltage  pulses  of  up  to  1200  V  with  I  ps  duration  were 
applied  with  a  high-power,  low-impedance  pulse  generator. 


The  pulse  repetition  rate  was  chosen  to  be  around  5  Hz.  For 
rates  higher  than  10  Hz,  lasing  breaks  down  due  to  sample 
heati.ig.  The  experiments  were  performed  at  4.2  K,  so  the 
sample  was  immersed  in  liquid  helium  at  the  center  of  a 
superconducting  solenoid.  The  FIR  radiation  was  guided 
through  a  brass  light  pipe  ( 12  mm  bore)  to  the  detectors, 
which  were  placed  in  the  same  cryostat  below  the  sample  in 
the  center  of  a  second  superconducting  magnet.  To  check 
the  range  of  electric  and  magnetic  fields  where  lasing  can  be 
achieved,  the  radiation  was  detected  with  a  broadband 
Ge.Ga  detector  with  peak  sensitivity  at  about  1 10  pm.  A 
magnetic  field  tunable  GaAs  detector  was  used  for  the  de¬ 
tailed  analysis  of  the  emission  spectra.  The  Zeeman-split  ls- 
2 p*  shallow  donor  line  was  used  for  detection.4  This  line  is 
linearly  tunable  between  45  and  130  cm'1  with  magnetic 
fields  between  I  and  6.5  T.  The  GaAs  detector  having  an 
intrinsic  linewidth  of  0.25  cm  “ 1  is  thus  used  as  a  high-resolu¬ 
tion  spectrometer. 

For  samples  without  mirrors,  lasing  was  observed  for 
magnetic  fields  between  I  and  2  T  and  electric  fields  between 
1000  and  2500  V/cm.  The  output  power  was  estimated  to 
reach  200  mW  by  comparing  the  magnitude  of  the  signal 
with  the  signal  obtained  from  a  FIR  gas  laser  and  an  InSb- 
Landau  source.  Figure  1  shows  spectra  ( intensity  versus  fre¬ 
quency)  from  two  samples  for  different  values  of  the  electric 
and  magnetic  fields.  The  decreasing  background  level  is  due 
to  the  variation  of  the  detector  characteristic  with  magnetic 
field.  The  spectra  consist  of  10-20  lines  with  a  regular  mode 
spacing.  Their  positions  remain  fixed  with  changes  in  ap¬ 
plied  fields.  The  measured  linewidth  is  about  0.4  cm  1  We 

TABLE  I  Sample  dimensions  and  observed  mode  ^pacing*  /  is  ihe  sjmple 
length  (til,)  direction  | .  d is  the  distance  between  the  electrical  coniacts.  a 
is  Ihe  width,  and  Av  is  the  measured  mode  spacing. 


j  At 

Sample  /  d  (mm)  (cm  ') 
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FIG.  I.  Spectrum  of  stimulated  emission  fer  different  eleclricand  magnetic 
fields.  Sample  I:  £=  1 .8  kV/cm.  B  =  1.4  T  Sample  3:  E  =  1.9  kV/crn. 
B  =  1  7  T  The  insert  schematically  shows  the  ray  path  in  the  Ge  crystal. 

conclude  that  the  real  width  of  the  emission  is  of  the  same 
order  of  magnitude  as  the  detector  linewidth  ( 0.25  cm  ~ 1 ) . 
By  varying  the  electric  and  magnetic  fields  it  is  possible  to 
obtain  lasing  between  75  and  1 10  fi m.  There  is  a  general 
trend  towards  emission  at  shorter  wavelengths  for  both 
higher  electric  and  higher  magnetic  fields.  In  contrast  to  An¬ 
dronov  et  al.*'  and  Komiyama,’  no  long  wavelength  lasing 
( 1 50-250 /im )  could  be  found  in  our  samples,  which  is  prob¬ 
ably  because  of  the  somewhat  higher  carrier  concentration. 
Andronov  et  al.'  and  Komiyama  et  air  suggested  that  the 
lasing  modes  are  due  to  total  reflections  at  the  sample  sur¬ 
faces,  since  without  mirrors  the  quality  factor  of  the  cavity  is 
far  too  low  to  allow  longitudinal  or  transverse  modes  to  os¬ 
cillate  (the  Ge  surface  has  36%  reflectivity  at  normal  inci¬ 
dence).  According  to  Fig.  1  the  observed  mode  spacings  are 
1.74  and  1.20  cm-1  for  samples  1  and  3,  respectively,  and 
1 .40  cm  ~ 1  for  sample  2  ( not  shown ) .  Comparing  the  spac¬ 
ings  with  the  dimensions  of  the  samples,  it  appears  to  be  very 
likely  that  the  width,  a,  is  the  crucial  parameter  determining 
the  mode  spacing  (see  Table  I). 

The  mode  spacings  can  be  explained  by  considering  an 
electromagnetic  wave  in  the  germanium  crystal  ( see  insert  in 
Fig.  1 )  propagating  with  components  of  the  wave  vector  in 
the  direction  of  /  ( z  direction)  and  a  ( x  direction):  k 
—  k,e,  +  k,  e, ,  where  e,  and  e,  are  the  respective  unit  vec¬ 
tors.  The  transverse  resonance  condition  is  given  by10 
k,n  =  kn  sin  8  =  Mir /a,  where  n  is  the  refractive  index,  M 
is  an  integer,  and  8  is  the  angle  between  k  and  the  longitudi¬ 
nal  axis  of  the  crystal.  The  analogous  longitudinal  resonance 
condition  is  not  considered,  because  it  would  yield  a  very 
narrow  mode  structure  which  is  not  observed  in  the  experi¬ 
ment.  In  addition,  we  have  the  geometrical  condition 
tan  6  =  Na/l  (A  integer),  which  ensures  that  the  path  of  the 
ray  is  closed  within  the  sample.  This  leads  to  a  mode  spacing 
(in  wave  numbers)  Av  =  ( Ian  sin  8)~'.  9  must  satisfy  the 
condition  14.5*  <  8  <  75.5*,  because  14.5’  is  the  angle  of  total 
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reflection  for  germanium.  One  can  assume  that  modes  with 
the  smallest  possible  8  have  the  lowest  losses,  because  they 
undergo  the  fewest  reflections.  Thus,  for  each  sample  the 
mode  spacing  is  determined  first  by  choosing  the  angle  8 
with  the  smallest  possible  A\  and  then  this  angle  is  used  in  the 
expression  for  Av.  This  procedure  gives  a  mode  spacing  Av 
of  1.77  cm- 1  {N  =  6),  1.41  cm-1  (N  =  5),  and  1.14  cm-1 
(N  —  4),  for  samples  1,2,  and  3,  respectively,  in  good  agree¬ 
ment  with  the  measured  values.  This  interpretation  is  con¬ 
firmed  by  the  work  of  Andronov  et  a/.*  who  reported  emis¬ 
sion  of  a  single  line  by  using  a  very  thin  (a  =  0.1  mm) 
sample,  where  the  mode  spacing  is  as  wide  as  the  gain  region. 

After  attaching  external  mirrors  to  sample  3  (see  insert 
in  Fig.  2),  the  region  of  electric  and  magnetic  fields  where 
lasing  occurs  slays  roughly  the  same.  The  detector  signal  is 
also  of  the  same  order  of  magnitude  with  and  without  mi; 
rors,  in  contrast  to  Komiyama  and  Kuroda,7  because  our 
light  pipe  is  only  24  cm  long  and  therefore  the  attenuation  of 
the  radiation  is  negligible  also  for  the  strongly  deflected  light 
from  the  mirrorless  samples.  The  observed  spectra,  however, 
exhibit  a  dramatic  change.  Figure  2  shows  spectra  of  sample 
3  for  different  electric  and  magnetic  fields,  recorded  with  the 
GaAs  detector.  A  considerably  reduced  number  of  lines  ap¬ 
pears  with  wider  mode  spacing  and  higher  amplitude,  ac¬ 
companied  by  some  other  less  dominant  structure.  Varying 
the  electric  and  magnetic  fields  allows  us  to  switch  on  and  off 
lines  on  the  long  or  short  wavelength  side  of  the  spectrum 
due  to  a  change  in  the  gain  spectrum.  In  this  way  it  was 
possible  to  obtain  operation  of  a  single  line  (lower  part  of 
Fig.  2).  This  is  the  first  report  on  forcing  a  multimode  p-Ge 
laser  to  oscillate  on  a  single  line  by  attaching  external  mir¬ 
rors.  The  emission  wavelength  is  87.3  pm  ( 1 14.5  cm  “ 1 ),  the 
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FIG.  2.  Emission  spectra  from  sample  3  after  attaching  external  mirrors  for 
different  electric  and  magnetic  fields.  The  insert  shows  the  cavity  configura¬ 
tion. 
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adjacent  line  at  83.3  /rm  ( 120.0  cm" 1 )  is  just  at  threshold. 
As  is  seen  from  Fig.  2  the  single  line  contains  approximately 
the  same  total  intensity  as  the  three  lines  together.  However, 
one  cannot  exclude  the  coexistence  of  several  cavity  modes 
within  this  line,  because  their  spacing  would  be  below  the 
resolution  of  the  detector.  The  spectra  observed  by  An¬ 
dronov  et  at.*  from  samples  with  external  mirrors  showed  a 
very  irregular  and  complicated  structure.  Komiyama  and 
Kuroda7  reported  some  changes  in  the  spectrum  and  were 
able  to  reduce  the  number  of  oscillating  lines  by  sandwiching 
a  thin  layer  of  poly-4-methylpent-l-ene  (manufactured  un¬ 
der  the  trade  name  TPX  by  Imperial  Chemical  Industries) 
between  the  sample  and  one  mirror. 

The  mode  spacing  in  Fig.  2  is  equal  to  Av  =  5.5  cm  **  \ 
which  is  too  large  for  fundamental  longitudinal  or  transverse 
cavity  modes.  A  possible  explanation  is  the  formation  of  an 
intracavity  interference  filter  by  the  two  Ge  plates:  Since 
they  are  not  subject  to  the  electric  field  like  the  sample,  they 
are  expected  to  have  a  slightly  different  index  of  refraction 
and  also  the  thin  layer  between  sample  and  mirror  induces  a 
change  of  the  refractive  index,  which  provides  the  possibility 
of  Fabry-Perot  like  reflections.  Thicknesses  of  d ,  =  0.9  mm 
and  d:  =  1.1  mm  (the  slight  curvature  of  the  convex  plate  is 
neglected)  would  give  a  spacing  of  only  1.4  and  1.1  cm-1, 
respectively.  If,  however,  the  combined  resonance  condition 
for  both  plates  is  considered,  the  spacing  is  calculated  to  5.5 
cm*1.  This  result  is  based  upon  a  ratio  of  d2/d,  equal  to  5/4. 
Considering  an  uncertainty  of  the  measured  thicknesses  of 
0.05  mm  and  a  low  finesse  of  the  Fabry-Perots,  one  gets 
spacings  between  4  and  6  cm  - '.  A  definitive  explanation  of 
the  present  spectra  would  require  knowledge  about  the  re¬ 
fractive  index  of  the  lasing  sample  as  well  as  the  influence  of 
the  thin  layer  between  sample  and  mirrors.  The  investigation 
of  different  mirror  configurations,  especially  with  different 
thicknesses,  appears  necessary. 

We  also  checked  the  emission  when  the  plane  mirror 


had  a  reflectivity  of  only  90%.  In  this  case  no  lasing  was 
achieved,  which  confirms  the  interpretation  as  waveguide¬ 
like  modes  in  the  mirrorless  configuration:  These  modes  are 
destroyed  by  attaching  mirrors,  which  prevents  total  reflec¬ 
tion  on  both  end  faces.  Considering  only  the  mirror  output 
and  no  other  losses  the  gain  is  calculated  to  be  about 
5x  10* '  cm~ '.  which  indicates  that  a  high-quality  factor  of 
the  cavity  is  necessary. 

In  summary,  we  have  quantitatively  explained  the  mode 
structure  of p-Ge  lasers  without  external  mirrors  in  terms  of 
waveguidelike  modes.  With  external  mirrors  laser  operation 
on  a  single  iine  is  demonstrated. 
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HOT  CARRIER  QUANTUM  DISTRIBUTION  FUNCTION  IN  CROSSED  ELECTRIC  AND 

MAGNETIC  FIH  .PS 
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Institut  fiir  Experimentalphysik.  Universitat  Innsbruck, 
A-6020  Innsbruck.  Austria 


ABSTRACT 

The  problem  of  the  hot  carrier  distribution  function  in  crossed 
electric  and  magnetic  field  is  treated  for  free  carriers  in  a 
parabolic  band.  To  account  for  the  discrete  structure  of  the  energy 
spectrum  (Landau  levels)  a  fully  quantum  mechanical  treatment  is 
employed.  The  master  equation  for  the  diagonal  elements  of  the  density 
operator  is  solved  for  three  Landau  levels  considering  optical  and 
acoustic  deformation  potential  scattering.  Numerical  results  are 
presented  for  the  case  of  the  light  holes  in  p-Germanium.  The  total 
populations  of  the  different  Landau  levels  are  discussed  as  well  as 
the  one-dimensional  distribution  functions  in  the  Landau  levels,  and 
their  dependences  on  electric  and  magnetic  field  strengths  are 
studied.  Special  attention  is  paid  to  the  occurrence  of  a  population 
inversion  between  Landau  levels  and  to  a  comparison  of  the  present 
approach  to  the  classical  picture  of  streaming  motion. 

^Present  address:  Bell  Comnunl cat ions  Research,  Inc.,  Red  Bank. 

New  Jersey  07701 

PACS  72.20  My,  72.20  Dp 


1 . INTRODUCTION 


The  transport  problem  of  carriers  in  semiconductors  under  presence  of 
crossed  electric  and  magnetic  fields  has  been  discussed  in  the 
literature  for  many  years  (for  review  articles  see  Refs.  1  and  2).  In 
most  of  the  previous  work  the  electric  field  was  assumed  to  act  as  a 
weak  perturbation,  i.e.  its  main  role  was  to  probe  the  system. 
According  to  this,  the  central  aim  was  to  calculate  the  diagonal  part 
of  the  conductivity  .  a '  ,  using  the  equilibrium  carrier  distribution 
function.  In  this  way  it  was.  for  example,  possible  to  explain 

3 

magnetophonon  oscillations  successfully  .  Another  part  of  the  work 

published  in  the  past  was  devoted  to  hot  electron  transport  in  crossed 
4-9 

fields  .  which  means  the  presence  of  a  rather  high  electric  field 

and  requires  the  calculation  of  the  nonequilibrium  distribution 

function.  Most  of  the  authors  assumed  a  very  small  magnetic  field  to 

justify  the  use  of  a  classical  distribution  function,  calculated  by 

7  8  9 

means  of  the  Boltzmann  equation  or  Monte  Carlo  techniques  ’  .  The 

other  limiting  case,  the  extreme  magnetic  quantum  limit,  has  been 

studied  in  the  frame  of  quantum  transport  equations*®.  In  this 

situation  only  the  lowest  Landau  level  (LL)  is  assumed  to  be  occupied. 

Only  very  little  attention  has  been  paid  to  the  intermediate  range, 

where  several  LL’ s  are  occupied,  but  the  discrete  structure  of  the 

11  12 

energy  spectrum  is  nonetheless  Important 

Very  recently,  this  case  turned  out  to  be  of  great  practical 

interest  due  to  the  observation  of  population  inversion  between  LL's 

13 

and  stimulated  cyclotron  resonance  emission  in  p-Germanium  .  This 


phenomenon  (also  discussed  as  negative  dynamic  conductivity)  has  been 
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predicted  already  fifteen  years  ago  .  The  build-up  of  an  inverted 

distribution  function  has  been  described  in  terms  of  streaming  motion 

18 

and  accumulation  of  carriers  in  momentum  space  .  This  concept  is 

based  on  the  classical  motion  of  carriers  in  crossed  fields 

considering  a  sufficiently  strong  interaction  with  optical  phonons. 

The  classical  trajectories  of  the  carriers  are  circles  in  velocity 

space  centered  at  v=-E/B.  The  energy  of  the  optical  phonons,  is 

represented  by  another  circle,  v  =  vQp=  ^2hfc>op/m*  (m*. . .effective 

mass),  centered  at  the  origin  v=0.  For  high  electric  fields.  E/B»vQp, 

all  carrier  trajectories  cross  the  circle  v=v  ,  so  the  carriers 

op 

continuously  get  accelerated  and  scattered  back  by  the  emission  of  an 

optical  phonon  (streaming  motion).  For  higher  magnetic  fields 

(E/B£vop)  a  part  of  the  trajectories  is  closed  without  crossing  the 

circle  v=v  .  The  corresponding  carriers  are  not  affected  by  optical 
•  * 

phonon  scattering  and  accumulate  in  a  limited  area  of  v-  or  momentum 
space,  because  they  can  be  removed  from  there  only  by  acoustic  phonon 
or  impurity  scattering,  which  are  both  much  slower  than  the  optical 
mode  scattering.  It  has  been  shown  that  a  torus-type  distribution 
function  can  arise  in  this  accumulation  area1  ,  which  corresponds  to  a 
population  inversion  between  LL.  Flg.l  shows  electron  trajectories  in 
velocity  space.  The  accumulation  area  is  bounded  near  the  second 
circle.  The  circles  can  also  be  regarded  as  Landau  levels,  then  the 
n=0  and  n=l  LL’s  lie  in  the  accumulation  area,  whereas  the  carriers  in 
all  higher  LL's  interact  strongly  with  optical  phonons  (streaming; 
indicated  by  the  dashed  part  of  the  cycles).  Up  to  now,  most  of  the 
calculations  analyzing  this  situation  dealt  with  the  limit  of 
nonquantizing  magnetic  fields;  that  means  the  use  of  a  classical 
distribution  function  and  the  use  of  the  zero  magnetic  field 
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scattering  rates  .  Although  much  physical  insight  can  be  gained  by 

this  approach,  it  would  be  interesting  to  illuminate  the  problem  from 

the  opposite  side  by  calculating  the  occupation  of  LL's  in  the  frame 

of  quantum  mechanical  transport  theory.  A  first  attempt  has  been  made 
12 

by  Helm  and  Gornik  .  who  solved  the  quantum  mechanical  transport 
equation  for  three  LL’s,  though  neglecting  the  k^-  dependence 
(z. .direction  of  magnetic  field)  of  the  distribution  function. 

In  this  paper  we  want  to  extend  this  approach  by  solving  the 

master  equation  for  the  density  operator  Including  the  k  -  dependence 

z 

of  the  distribution  function.  We  calculate,  for  the  first  time,  the 
hot  carrier  quantum  distribution  function  in  crossed  electric  and 
magnetic  fields  for  three  LL  without  any  weak-electric  field 
approximation.  The  basic  equations  and  the  way  of  their  solution  will 
be  presented  in  Sec. 2.  Analytical  results  of  the  scattering  rates  for 
acoustic  and  optical  deformation  potential  interaction  are  given.  In 
Sec. 3  the  numerical  results  are  shown  for  the  case  of  the  light  holes 
of  p-Cermanium.  The  total  populations  of  the  different  LL’s  are 
discussed  as  well  as  the  one-dimensional  distribution  functions  in  the 
LL’s.  and  their  dependences  on  electric  and  magnetic  fields  and  phonon 
coupling  strengths  are  studied.  Special  attention  is  paid  to  a 
comparison  of  the  quantum  mechanical  picture  with  the  classical 
picture  of  streaming  motion. 

2.  THEORY 


(a)  Basic  expressions 
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fe  consider  the  Hamiltonian  for  a  carrier  with  charge  e  and  effective 

mass  m*  in  crossed  electric  and  magnetic  fields  interacting  with 

.  6,10. 

phonons 


» *  sr  ♦  Hph  *  H.i-ph 


Here  H  ,  is  the  free  phonon  term: 
ph 


H  .  =  2  bu  (b  b  +  tf) 
Ph  J  q'q  q 


(1) 


(2) 


where  the  b^  and  b^  are  the  phonon  creation  and  destruction  operators, 
respectively,  and  H@j_  ^  is  the  electron-  (hole-)  phonon  coupling: 


Hel-ph  *  1  |  C<«>  % 


(3) 


C(q)  depends  on  the  specific  nature  of  the  electron-phonon  coupling. 
Eigenvalues  and  eigenfunctions  of  the  electronic  part  of  the 
Hamiltonian  (1?  in  z-direction,  if  in  x-direction)  are  given  by 

Ey  =  (n+$4)  huc  +  ft2k2/2m*  -  eEX  +  H  m*(E/B)2  (4) 


where  the  index  v  summarizes  the  quantum  numbers  n,  k  ,  and  k  , 
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u  as  je|B/m*.  and  X  =  (hk  /eB  +  m*E/eB2). 
c  y 


*  (r) 

uv  ' 


.‘V  * 

nv  1  ' 


(5) 


where  1  =v¥i/  |e  |B  and  ♦ 


is  the  n-th  harmonic  oscillator  wave 


function.  The  basic  quantum  transport  equation  is  the  equation  of 
motion  for  the  density  operator,  the  Liouvllle  -  von  Neumann  equation 


if»  §£  =  [«.  p]  (6) 

Here,  p  is  the  density  operator  of  the  total  electron  -  phonon  system. 
The  electron  distribution  function  is  given  by  the  diagonal  elements 
of  the  reduced  density  operator,  f  .  which  is  obtained  by  taking  the 
trace  over  the  phonon  states: 


f 


v 


~  Trph<w lp lu> 


(7) 


Assuming  no  correlation  between  electron  and  phonons  in  the  diagonal 

elements  of  p  and  a  weak  electron  -  phonon  coupling  the  Pauli  master 

equation  can  be  derived  for  the  diagonal  elements  of  the  reduced 
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density  operator  (neglecting  the  Pauli  exclusion  principle) 


D  U 


-  f  W  . 

U  DU 


(8) 


The  W  ,  are  the  scattering  probabilities  from  state  u  to  state  v' 
according  to  Fermi’s  golden  rule  : 


l<»|C(q)  «iqV|»'>r  [V(VE.'*t',WV1>!<VE.’-|“,>' 


q'  q 


Nq  is  the  Planck  distribution  for  the  phonons,  the  matrix  elements  are 


given  by  (see  e.g.  Ref.  21) 


( 

I 


Mounuryjr^  vihjiw  ”J rw  ■>  mn  vsrrjnsr  v  vnr',  v*  .t.1^  ■'■«  ■■«  r*  t*>  gw  it  '.«.■  gw  '.t.twt'.t.v.v  vv  v 


|<»|.l3?|.'>|*  -  |Jnn.(Vqy)|a 

•ith  |Jnn.(£)|2-jiTE"'‘ne-f 


S-V',  6ky'ky"S’ 

[L”,_n(f)]2  (n£n* ) . 


(9) 


and  f  =  l2(q^  +  q^)/2.  The  n  are  associated  Laguerre 
22 

polynomials  .  It  is  worth  noting  that  the  electric  field  only  appears 
in  the  energy  conserving  5-function  and  not  in  the  squared  matrix 
elements. 


(b)  Steady  state  solution 


We  now  proceed  to  get  a  steady  state  solution  of  the  master  equation 
(8)  which  then  reads 


0  =  2.  f  . 
u  1) 


W 


V  V 


-  f  W  , 
V  vv 


(11) 


In  the  following  we  want  to  consider  the  situation  where  the  current 

through  the  sample  is  constant.  Therefore  we  assume  a  spatially 

homogeneous  distribution  function.  Thus,  the  fv  have  to  be  independent 
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of  X,  and  therefore  of  k  ’  .  This  enables  us  to  perform  the 
summation  over  ky  directly  and  we  are  left  with  a  system  of  N  coupled 
integral  equations  (for  N  LL)  in  the  variable  k^: 


Xdk*  f  ,(k’)  W  .  (k\k  )  -  f  (k  )  W  ,(k  ,k*) 
z  n  v  z ’  n  nv  z  zJ  nv  z '  nn  v  z  z* 


(12) 


The  rates  V  ,(k  ,k')  are  now  also  k  -integrated.  The 
nn  z  z  y 

functions  for  the  different  LL’s  can  be  expressed  as 


distribution 


fc»  l*  ■»<  U»  fc’  i|»  ||*  ***.14*  _•»•.***  j.* 


W  -  r»  *.<k  T  i-Wki  fn-(ki>  Vn<kikz>  • 


n  nn  '  z 


where  W  .(k  )  =  fdk'  W  ,(k  .k*). 

nn  '  zJ  z  nn  v  z  zy 


The  proper  normalization  condition  is  I/dk  f  (k  )  =  1. 

n  z  nv  z' 


(c)  Transition  rates 


Regarding  the  explicit  calculation,  we  restrict  ourselves  to  three 

Landau  levels.  A  larger  number  would  be  desirable,  but  since  the 

integrals  are  not  symmetric  concerning  ^  and  due  to  the  electric 

field,  no  simple  general  expression  for  the  transition  ’  rates 

W  ,(k  ,k*)  could  be  found.  However,  three  LL’s  should  be  enough  to 
nn  z  z 

show  some  Interesting  features  of  the  carrier  distribution. 


For  the  sake  of  simplicity,  we  consider  acoustic  and  optical 
deformation  potential  interaction  in  the  calculation,  which  is 
appropriate  i.e.  for  the  light  holes  of  p-Germanium.  The  polar-optical 
phonon  interaction  would  also  be  of  Interest  in  order  to  treat  the 
electrons  of  GaAs  or  InSb,  but  the  different  momentum  dependence  needs 
a  numerical  evaluation  of  the  scattering  rates  in  this  case.  The 
coupling  constants  are  given  by  (see  e.g.  Ref. 23) 


,  n2  hJ(DK) 

Irf'eo Is  -  |rr!  •  lcop(i)la  =  5vT^ 


2V  p  hw 


where  is  the  optical  phonon  energy,  D  is  the  acoustic  deformation 


8 


Lr|V| 


m 


potential,  D(K  characterizes  the  optical  deformation  potential,  p  is 
the  mass  density,  s  is  the  sound  velocity,  and  V  is  the  crystal 
volume.  We  adopt  three  more  approximations  for  the  carrier-phonon 
interaction-* 

(1)  We  consider  only  the  spontaneous  emission  of  optical  phonons 
(Nqp=  0)  and 

(2)  apply  the  equipartition  approximation  for  acoustic  phonons 

(1«  N  =  kT/hu  ). 
v  ac  q' 

(3)  Furthermore  we  regard  the  scattering  by  acoustic  phonons  as 
elastic  (fc>q=  0  in  the  energy  conserving  6-function) . 

These  approximations  can  assumed  to  be  valid  at  intermediate 

temperatures  (T  =;  10  -  80  K) .  As  a  consequence,  the  optical  and 

acoustic  scattering  rates  exhibit  the  same  q-dependences .  The  explicit 

transition  rates  between  the  three  lowest  LL’s  are  then  given  by  (the 

argument  (k  ,  k’)  is  omitted) 
z  z 


w, 


oo 


=  A  e 


-Ha4 


-ifes 


wio(oi)=  A  <1/2>  <1+fea>  e 

Wn  =s  A  (1/4)  (3-2ta+e4)  e 

W20(02)=  A  {1/8)  C3*2*2**4)  •' 

W21(i2)=  A  (1/16)  (7+9ea-5e4+t6)  e 


-ifea 

-ifea 


-Vie2 
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=  A  (1/64)  (41-52ea+50e4-12e6+fe8)  e 


-ife2 


(15) 


Here  e  is  the  respective  difference  of  energies  in  state  v  and  v’ 
divided  by  eEl  for  each  transition  rate,  which  is 


h2(k’ a-ka)/2m*  +  (hu)  for  W00,  Wlt.  W22. 
z  z  op 

ha(k'a-ka)/2mw  -  to  +  (f*>)  for  Wto.  W21, 
z  z  c  op 
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I 

i5i 


h2(k,2-k2)/2m*  +  ftu  +  (ftu  ) 
v  Z  z'  C  '  Op7 

ft2(k*2-k2)/2m«  -  2fi«  +  (ftu  ) 
y  z  z'  c  v  opJ 

h2(k'2-k2)/2m*  +  2hw  +  (hw  ) 
v  z  z'  c  y  op' 


for  Wqi*  Wj2* 
for  W2o>  suid 
for  Wq2 • 


The  hu  in  parentheses  applies  only  for  the  optical  phonon  rates.  The 
op 

constant  A  is  different  for  acoustic  and  optical  scattering,  namely 


^ac _ D2kT  BvSr 

4ir2p  s2h2  E  1 


and  AH= 


op_  Bv^ 


8i r2p  hu  E  1 

op 


The  crucial  quantity  in  the  transition  rates  (15)  is  e,  the  ratio  of 

the  energy  difference  between  initial  and  final  state  and  the 

electrostatic  energy  eEl.  For  weak  electric  fields  the  energy 

difference  has  to  be  practically  zero  to  yield  a  reasonable  magnitude 

of  the  rates.  If  this  is  not  the  case,  the  rates  go  to  zero.  For 

higher  electric  fields  the  energy  conservation  is  seemingly  not  that 

strict;  as  long  as  the  energy  difference  is  of  the  order  of  eEl,  there 

is  a  considerable  probability  for  the  respective  transition. 

Physically,  the  missing  amount  of  energy  is  supplied  by  spatial 

transitions,  which  are  possible  due  to  the  overlap  of  the 

12 

wavefunctlons  :  A  change  of  the  spatial  coordinate  X  implies  a  change 

of  the  total  energy  due  to  the  term  eEX  in  Eq.(4).  Hence  at  low 

electric  fields  the  transitions  between  Landau  levels  are  dominated  by 

transitions  associated  with  a  change  of  the  quantum  number  k  ;  at  high 

electric  fields  the  contribution  due  to  spatial  transitions  (change  of 

quantum  number  k^)  becomes  more  important.  Typical  transitions 

W  ,(k  ,  k')  are  sketched  in  Fig. 2  for  acoustic  phonons, 
nn  z  z 


The  k’-integrated  rates,  W  ,(k  ),  which  are  needed  in  Eq.(13), 
z  nn  z 

are  calculated  numerically.  The  main  procedure  now  is  the  solution  of 
the  three  coupled  integral  equations  (13)  with  the  transition  rates 
(15).  This  is  done  by  simultaneous  iteration  of  the  three  equations. 
The  convergence  is  rather  fast,  the  fifth  iteration  is  accurate  up  to 
0.1%. 


At  the  end  of  this  section  we  want  to  make  some  comments  on 
previous  work  performed  on  the  basis  of  the  master  equation  (11). 
Inoue  and  Yamashita^  expanded  the  energy  conserving  6-function  in 
Eq.(9)  up  to  second  order  in  the  electric  field.  They  omitted  the 
^-integration  and  evaluated  instead  the  integral  at  the  k^-value 
which  gives  the  major  contribution.  This  leds  to  a  system  of  N  second 
.  order  differential  equations.  Within  their  approximations  they  were 
able  to  give  general  expressions  for  the  distribution  function  in  an 
arbitrary  number  of  LL’s.  Unfortunately  no  numerical  study  was 
performed. 

Calecki  et  al  transformed  the  master  equation  (11)  into  a 

Fokker-Planck  equation  and  studied  the  extreme  quantum  limit  (one  LL) 
analytically.  They  got  a  Maxwellian  distribution  function  in  certain 
limits  and  derived  an  expression  for  the  electron  temperature. 
Moreover,  they  pointed  out  that  the  (stationary)  distribution  function 

depends  only  on  ratio  of  the  coupling  strengths  of  different 

,  20 
scattering  processes 


In  comparison,  we  can  not  include  an  arbitrary  number  of  LL’s, 
but  do  not  make  any  approximation  concerning  the  electric  field  or  the 


igwin  uv’w’v’jvjy^rv.  r  v»  v»HtyT*v  TJ  vjy.v.v v  t»  ■tv  ^7^ 
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integration,  like  Inoue  and  Yamashita.  Owing  to  numerical  methods  we 
are  able  to  calculate  the  distribution  functions  of  a  finite  number  of 
LL's  directly  from  the  master  equation. 


3.  RESULTS  AND  DISCUSSION 


$ 

£ 


The  central  quantities  obtained  from  the  above  calculation  are  the 

one-dimensional  distribution  functions  f  (k  )  for  the  different  LL*s. 

nv  zJ 

Another  interesting  quantity,  especially  regarding  the  appearance  of  a 
population  inversion,  is  the  total  population  in  every  LL 


f  =  Jdk  f  (k  ). 
nz  nv  z7 

Fig. 3  shows  the  dependence  of  the  populations  f0.  fi.  f2  on  the 

electric  field  for  a  constant  magnetic  field  (B  =  0.35  T  and  1  T).  The 

24 

parameters  are  chosen  suitable  for  the  light  holes  of  Germanium 
m*  =  0.043  mo.  p  =  5300  kg/m3 .  s  =  5400  m/s,  37  meV,  D  =  4.6  eV. 

D^K  =  9*10loeV/m.  According  to  Fig. 3a,  one  can  describe  the  behavior 
as  follows:  At  low  electric  fields  all  populations  are  near  0.33,  but 
fQ>  ft>  f2.  This  is  the  case,  because  we  assumed  acoustic  phonon 
scattering,  which  is  dominant  in  this  region,  to  be  totally  elastic.  A 
finite  acoustic  phonon  energy  (in  the  energy  6-function)  would 
increase  the  differences  in  the  populations.  When  the  electric  field 
is  Increased,  the  carriers  in  the  n=2  LL  start  to  interact  with 
optical  phonons,  which  causes  a  decrease  of  f2.  This  can  be  regarded 
as  an  onset  of  streaming  motion  of  the  carriers  in  the  n=2  LL,  or,  in 
other  words,  of  the  carriers  on  the  classical  trajectory  with  the 
largest  radius  in  k-  or  v-  space.  As  all  LL’s  are  connected  by  the 
respective  scattering  rates,  also  fQ  and  ft  are  influenced  and  a 


population  inversion  between  the  n=l  and  n=0  LL  occurs  in  a  narrow 
range  of  E.  On  further  increasing  the  electric  field  also  the  carriers 
in  the  n=l  LL  start  to  interact  strongly  with  optical  phonons 
(streaming),  and  f4  decreases.  For  very  strong  electric  fields  the 
carriers  in  all  LL’s  can  emit  optical  phonons  and  the  three 
populations  again  will  approach  the  value  0.33.  However,  there  remains 
the  condition  fQ>  ft>  f2.  because  phonon  emission  preferably  occurs 
into  lower  LL’s.  The  difference  f0-  ft  is  also  shown  in  Fig.  3, 
inversion  corresponds  to  f0~  fj  <  0.  The  situation  at  population 
inversion  turns  out  to  be  closely  analogous  to  the  classical  picture: 
The  n=2  LL  corresponds  to  the  streaming  carriers,  the  n=0  and  n=l  LL’s 
represent  the  accumulation  area  (see  Fig.l).  Due  to  the  specific 
nature  of  the  scattering  rates  an  inversion  occurs  in  the  accumulation 
area,  that  is  between  the  n=l  and  n==0  LL. 

Fig. 3b  shows  the  analogous  curves  for  B  =  1  T.  The  general 
behavior  is  similar,  but  no  population  inversion  occurs  any  more  at 
this  magnetic  field.  All  the  features  are  less  pronounced  and  shifted 
to  higher  electric  fields. 

In  the  calculation  the  lattice  temperature  has  been  taken  10  K, 
which  is  near  the  lower  limit  for  the  validity  of  the  high  temperature 
approximation  for  acoustic  phonons;  however,  at  this  low  temperature 
(equivalent  to  weaker  acoustic  phonon  coupling  -  see  Eq.(16))  all  the 
features,  especially  population  inversion,  are  more  pronounced.  It 
turns  out  that  the  relative  strengths  of  optical  and  acoustic  phonon 
coupling  Influence  critically  the  possibility  of  inversion:  At  a  too 
high  or  low  ratio  of  coupling  strengths  no  inversion  builds  up  any 
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more.  By  chance,  the  realistic  coupling  constants  turn  out  to  be 
nearly  optimal. 

Let  us  now  turn  to  the  k  -dependence  of  the  distribution 

z 

functions.  Ve  have  calculated  f  (k  )  for  magnetic  fields  of  B=0.35T 

n  z 

and  B=1T,  with  values  of  the  electric  field  which  appear  interesting 

in  view  of  Fig. 3.  For  B=0.35T  (Fig. 4).  these  are  E=lkV/cm  (acoustic 

phonon  scattering  dominant) ,  E= 1 . 9kV/cm  (inversion),  and  E=2 . 8kV/cm 

(all  LL's  streaming).  The  scale  on  the  abscissa  is  transformed  from  k 

z 

to  h3k3/2m»*.  because  in  this  way  a  Boltzmann  distribution  (with 
z 

carrier  temperature  T  )  is  represented  by  a  straight  line  (the  line  in 

c 

Figs. 4  and  5  corresponds  to  T  s  50  K).  This  temperature,  however,  is 

c 

not  comparable  to  the  temperature  derived  by  Calecki  et  al*®.  because 
their  expressions  are  valid  only  in  the  limit  of  high  lattice 
temperature  (kT»ftuop)  or  extremely  high  electric  fields  (eEl»ht»op). 

For  EatlkV/cm  (Fig. 4a)  most  of  the  carriers  are  concentrated  below 
a  h3k3/2mw-value  of  about  20  meV,  where  the  emission  of  optical 
phonons  starts  to  become  effective.  The  corresponding  kink  is  observed 
at  energies  somewhat  below  the  optical  phonon  energy  (37  meV)  because 
of  the  possibility  of  spatial  transitions.  As  the  energy  is  counted 
from  each  LL-edge,  optical  phonon  emission  from  higher  LL’s  is 
possible  at  lower  k  -values  than  from  lower  LL’s  (for  example  the 
effective  energy  in  the  n=2  LL  is  given  by  2hu +  h3k3/2m»).  At  low 
energies  the  distribution  function  is  very  hot,  because  in  our 
approximation  the  acoustic  phonons  do  not  carry  energy.  At  higher 
energies,  where  optical  phonon  emission  becomes  effective,  the 
distribution  Is  cooled  down. 


Fig. 4b  shows  f  (k  )  for  1.9  kV/cm.  where  according  to  Fig. 3  the 
n  z 

inversion  is  strongest.  This  "total"  Inversion  is  reflected  by  a 

"local”  inversion  in  k  -space  near  k  =  0.  It  should  be  pointed  out. 

z  z 

however,  that  a  local  inversion  near  the  band  edge  can  appear,  even  if 

the  total  populations  are  not  inverted  (see  below).  In  connection  with 

a  sufficient  amount  of  nonparaboliclty  in  a  real  system,  this  could 

lead  to  light  amplification,  too.  (Also  for  a  total  inversion 

nonparaboliclty  is  necessary  to  get  amplification,  because  in  the 

parabolic  case  the  ratio  of  induced  absorption  to  emission  is  always 

(n+l)/n  for  carriers  in  the  n-th  LL.)  The  fact  that  the  carrier 

distributions  are  inverted  for  large  k  -values  (ftak2/2m»  >  25  meV)  is 

not  essential,  because  the  distribution  function  there  is  at  least 

three  orders  of  magnitudes  smaller  than  at  the  band  edge.  Since  the 

distributions  are  roughly  parallel  to  the  straight  line,  a  mean 

temperature  of  about  50  K  can  be  found  for  them  in  spite  of  the 

Inversion.  At  E  =  2.8  kV/cm  (Fig. 4c)  all  carriers  can  interact 

strongly  with  optical  phonons  (streaming).  The  ^-distributions  are 

hotter  and  near  to  being  thermal  (T  =  80  K). 

c 

Fig. 5  shows  the  distribution  functions  for  B  =  1  T  and  E  *  2,  3 

and  4  kV/cm.  2kV/cm  (Fig. 5a)  represents  the  low  electric  field  region. 

The  higher  magnetic  field  compared  to  Fig. 4  is  reflected  by  the  larger 

energetic  distance  of  the  optical  phonon  kink  in  different  LL’s.  At 

3  kV/cm  (Fig. 5b)  there  appears  a  local  inversion  between  the  n=l  and 

n=0  LL  near  k  =0.  Fig. 5c  shows  the  situation  for  E=4kV/cra,  when  the 
z 

distributions  are  nearly  thermal ized  by  the  optical  phonon 


Interaction. 


In  this  paper  we  have  calculated  the  hot  carrier  distribution 
function  in  crossed  electric  and  magnetic  fields  using  the  Pauli 
master  equation.  Since  we  take  into  account  three  LL's,  our  approach 
bridges  the  gap  between  the  extreme  quantum  limit  and  quasi classical 
calculations.  The  total  populations  of  the  different  LL's  have  been 
discussed  as  well  as  the  one  dimensional  distribution  functions  within 
the  LL’s. 

At  low  magnetic  fields  (B  £  0.35  T)  we  find  an  inversion  in  a 
certain  range  of  electric  fields,  where  the  motion  of  the  carriers  in 
the  two  lower  LL's  is  dominated  by  acoustic  phonon  scattering,  whereas 
the  carriers  in  the  third  (ns2)  LL  already  interact  strongly  with 
optical  phonons.  However,  the  populations  are  sensitive  to  the 
relative  coupling  strength  of  acoustic  and  optical  phonon  interaction. 
At  low  electric  fields  the  one  dimensional  distribution  funtlons  in 
the  different  LL’s  show  a  sudden  decrease  at  a  certain  energy  level 
due  to  the  onset  of  optical  phonon  emission.  At  higher  fields,  when 
all  carriers  in  one  LL  can  emit  optical  phonons,  the  distribution 
function  is  nearly  thermal.  Therefore  a  common  electron  temperature 
can  be  defined  in  this  case. 

The  predictions,  however,  are  not  claimed  to  be  quantitatively 
reliable,  because  in  the  real  p-Ge  System  also  the  heavy  holes  (which 
we  have  neglected)  will  influence  the  distribution  of  the  ligth  holes. 
For  further  improvement  of  the  theory  it  will  be  necessary  to  include 
a  larger  number  of  LL's,  especially  at  lower  magnetic  fields. 
Nevertheless  the  present  approach  is  the  first  rigorous  solution  of 
the  master  equation  in  the  crossed  field  situation  and  shows  the 
appearence  of  a  population  inversion,  which  is  due  to  the  interplay  of 


REFERENCES 


1.  R.Kubo,  S.J. Miyake,  and  N.Hashitsume,  in  "Solid  State  Physics” 
Vol.17,  Ed.  F. Seitz  and  D. Turnbull  (Academic  Press,  New  York, 

1965) ,  p.270 

2.  L.M.Roth  and  P.N.Argyres,  in  "Semiconductors  and  Semimetals" 
Vol.l,  Ed.  R.K.Willardson  and  A. C. Beer  (Academic  Press,  New  York, 

1966) ,  p. 159 

3.  R.L. Peterson  in  "Semiconductors  and  Semimetals”  Vol.10,  Ed. 
R.K.Villardson  and  A. C. Beer  (Academic  Press,  New  York,  1975), 

p.221 

4.  E.Yamada  and  T. Kurosawa.  J.Phys.Soc.  Japan  34.  603  (1973) 

5.  H. Parti,  W. MUller,  F.Kohl,  and  E.Gornik.  J.Phys.C  11,  1091  (1978) 

6.  H.F.Budd,  Phys.Rev. 175,  271  (1968) 

7.  F.Brosens  and  J.T.Devreese,  Solid  State  Conmiun.44.  597  (1982) 

8.  P.Varmenbol,  F.M.Peeters,  and  J.T.Devreese,  Phys.Rev.B  33,  1213 

(1986) 

9.  F.Abou  El-Ela  and  B.K. Ridley,  at  the  "5th  Int  Conf.  on  Hot 
Carriers  in  Semiconductors,  Boston  1987”,  to  be  published  in 
Solid  State  Electronics 

10.  D.Caleckl,  C.Lewiner,  and  P.Nozieres,  J.Phys.  (Paris)  2§,  169 
(1977) 

11.  S.Inoue  and  J.Yamashita,  Progr.Theor.Phvs.42.  158  (1969) 

12.  M.Helm  and  E.Gornik,  Phys.Rev.B  32.  7459  (1986) 

13.  Yu. B. Vasil ’ev  and  Yu. L. Ivanov,  in  "Proceedings  of  the  18th  Int 
Conf.  on  the  Physics  of  Semiconductors.  Stockholm,  1986”,  Ed. 
O.Engstrom  (World  Scientific  Publishing,  Singapore,  1987)  p.1659 

14.  H.Kaeda  and  T. Kurosawa,  in  "Proceedings  of  the  11th  Int.  Conf.  on 


the  Physics  of  Semiconductors,  Warsaw,  1972",  Ed.  Polish  Academy 
of  Sciences  (PWN-Polish  Scientific  Publishers,  Warsaw,  1972), 

p.602 

15.  Ya.I.Al'ber,  A. A. Andronov,  V. A. Valov.  V. A. Kozlov,  and 

I. P. Ryazan tseva.  Solid  State  Comnun.  19.  955  (1976) 

16.  T. Kurosawa,  Solid  State  Coimun.  24,  357  (1977) 

17.  A. A. Andronov.  V. A. Valov,  V. A. Kozlov,  and  L.S.Masov,  Solid  State 
Connun. 36.  603  (1980) 

18.  S.Komiyama,  Adv.Phys.  31.  255  (1981) 

19.  L.  Van  Hove.  Physica  21.  517  (1955);  23.  441  (1957) 

20.  D.Calecki,  in  "Physics  of  Nonlinear  Transport  in  Semiconductors", 
Ed.  D.K. Ferry.  J.R. Barker,  and  C.Jacoboni  (Plenum  Press,  New 
York,  1980),  p.289 

21.  J.R. Barker.  J.Phys.C  §.  1657,(1972) 

22.  I.S.Gradsteyn  and  I.N.Ryzhik,  "Tables  of  Integrals.  Series  and 
Products"  (Academic  Press,  New  York,  1965) 

23.  L.Reggiani  in  "Hot  Electron  Transport  in  Semiconductors”,  Ed. 
L.Reggiani,  Topics  in  Applied  Physics  §8  (Springer.  Berlin, 
1984),  p.7 

24.  "Semiconductors”,  Ed.  O.Madelung,  M. Schulz,  and  H. Weiss, 

Numerical  Data  and  Functional  Relationships  in  Science  and 
Technology,  Landolt-Bornstein  Series,  Group  III,  Vols.  17a  and 
17b  (Springer,  Berlin,  1982) 


FIGURE  CAPTIONS 


Flg.l:  Trajectories  of  electrons  in  velocity  space  in  crossed  electric 
and  magnetic  fields  (£.  in  x-direction,  if  in  z-direction) .  The 
circle  centered  at  v=0  represents  the  optical  phonon  energy. 
The  dashed  part  of  the  trajectories  indicates  the  emission  of 
an  optical  phonon. 

Fig. 2:  The  energies  of  two  LL's  are  plotted  versus  kz  and  the  carrier 
position  along  the  electric  field.  The  mechanism  of  acoustic 
phonon  transitions  is  illustrated. 

Fig. 3:  The  populations  fR  of  the  n=0  (solid  line),  n=l  (dashed)  and 
n=2  LL  (dashed-dotted)  and  the  population  difference  fD-f t 
(dotted)  as  a  function  of  the  electric  field  for  (a)  B=0.35T 

and  (b)  B=1T. 

Fig. 4:  The  one-dimensional  distribution  functions  f  (k  )  in  the  n=0 

n  z 

(solid  line),  n=l  (dashed)  and  n=2  LL  (dashed-dotted)  are 

plotted  versus  hak2/2m*»  for  B=0.35  T  and  (a)  E=lkV/cm,  (b) 

E=1.9kV/cm,  (c)  E=2 . 8kV/cm .  The  dotted  straight  line  represents 

a  Boltzmann  distribution  with  T  =50K. 

c 

Fig. 5:  Sam  as  Fig. 4.  but  for  B=1T  and  (a)  E=2kV/cra.  (b)  E=3kV/cm,  (c) 
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ABSTRACT 

Wa  have  observed  for  the  first  time  the  Smith-Pur call  affect  in  GaAs-GaAlAs- 
hatarostructuras.  A  grating  structure  induces  far  infrared  (FIR)  amission  from 
drifted  carrier  distributions,  which  is  angle  dependant  and  allows  the  determina¬ 
tion  of  intrinsic  carrier  properties. 

INTRODUCTION 

The  two-dimension-iil  electrons  in  selectively  doped  GaAa/GaAlAs  hetero junctions 
exhibit  rixtremelv  high  mobilities  at  low  temperatures  /I —3/ .  The  electrons  are 
easily  heated  up  in  an  external  electric  field  leading  to  a  rapid  rise  in  electron 
temperature  /4-6/  and  a  reduction  of  electron  mobility  /!/.  No  detailed  analysis 
of  the  form  of  the  distribution  function  has  been  performed  for  extremely  high 
mobilities.  In  this  case  normal  Fermi  distributions  cannot  be  used  since  drift 
energies  comparable  with  the  Fermi  energy  can  be  achieved  already  at  low  electric 
fields. 

In  this  paper  we  present  a  direct  analysis  of  the  electron  distribution  function 
for  extremely  high  mobility  samples  as  a  function  of  the  external  electric  field. 
The  Smith-Purcell  effect  /8/  is  measured  in  a  semiconductor  for  the  first  time: 

A  periodic  potential  of  period  L  induces  an  energy  loss  of  a  drifting  carrier  dis¬ 
tribution  via  photon  emission.  The  angle  between  the  drift  and  the  grating  momen¬ 
tum  is  varied  by  using  different  sample  geometries.  The  geometry-dependent  Mission 
spectra  are  monitored  through  a  magnetic-f laid- tuned  XnSb-detector  / 9/. 

Through  the  selective  momentum  transfer  (tq  ■  2»/L)  the  spectral  analysis  of  the 
energy  loss  yields  detailed  information  on  the  drifted  electron  distribution  func¬ 
tion.  For  a  momentum  transfer  much  smaller  than  the  electron  momentum,  the  frequen¬ 
cy  u  of  the  emitted  light  is  directly  proportional  to  the  particle  wavevector  k: 


,u  *  ftkq/m*  ■  v^q 


(q  «  k) . 


m*  is  the  electron  mss.  Thus  the  emitted  spectra  reflect  directly  the  distribu¬ 
tion  of  the  carrier  velocities  vk.  In  the  following  analysis  the  electron  distri¬ 
bution  function  will  be  describee  as  a  shifted  Fermi  function: 

f(E)  -  1/(1  ♦exp(ft*(k-kD)V2m*-»)/T#),  (2) 

introducing  the  drift  momentum  k_  and  the  drift-dependent  electron  temperature  T 
aa  phenomenological  parameters.  OOr  model  assumes  that  the  distribution  function 
is  stationary,  depending  only  on  the  electric  field  and  the  sample  mobility.  The 
FIR  coupling  represents  only  a  negligible  cooling  process. 

THBORY  OF  THE  EMISSION  PROCESS 

FIR  emission  from  a  free  electron  system  requires  a  second  scattering  process  for 
momentum  conservation,  which  can  be  either  impurity  or  phonon  scattering  /10/.  An 
electroMgnetlc  calculation  of  the  FIR  emission  from  heated  carriers  through  a 
grating  was  recently  performed  in  ref.  4.  For  the  electric  field  driven  electrons, 
the  lateral  potential  modulation  u  cos(qx)  represents  a  coherent  scattering  source 
(the  lattice  vector  has  been  taken°in  the  x  direction) .  Formally,  it  can  be  des¬ 
cribed  as  one  ’extended*  impurity,  scattering  the  electrons  by  iq.  The  emission 
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spectrum  of  a  drifted  electron  distribution  depends  on  the  an9le  between  the  la-  i 
teral  structure  and  the  electric  field  and  thus  the  current  direction. 

we  calculate  the  emission  rate  in  the  lowest  order,  including  only  the  scattering  i 
processes  shown  in  Fig.  (la).  This  diagrammatic  approximation  can  be  controlled  i 
analytically  also  for  'badly*  shaped  electron  distributions.  Static  screening  is  j 
included  in  the  impurity  potential.  Evaluating  the  diagrams  Fig.  (la)  and  per-  j 
forming  the  sum  over  the  photons  the  energy  loss  rate  is  obtained  as: 


dW 
-  ■ 

e2 

J£si_ 

so 

Z  — Ur  /dkf(k*,k)  {1  -  f  (k;  iq,  k)  }  , 

dt 

4tt  j  m*c  n 

±  |k- i  -«  y  x  y  *  y 

k*  » 

X 

♦“<! 

* 

■ 

The  sum  over  (±)  describes  the  two  possible  momentum  transfers,  and  the  k  -inte¬ 
gration  has  already  been  performed  by  exploiting  the  delta  function  for  energy 
conservation.  The  influence  of  the  electron  distribution  f  (k  ,k  )  is  contained  in 
the  ky  sum,  which  is  evaluated  numerically.  ” 

Fig.  (1b)  shows  the  calculated  energy-dependent  emission  intensity  for  the  para¬ 
meters  of  the  sample  described  below,  for  an  electron  temperature  of  35  X.  The 
curves  are  weighted  by  u  in  order  to  account  roughly  for  the  detector  sensitivity. 
The  full  curve  describes  a  carrier  distribution  which  is  heated  by  an  electric 
field  normal  to  the  grating  wavevector  q.  In  this  direction  the  result  is  indepen¬ 
dent  of  the  drift  velocity  if  a  shifted  Fermi  distribution  is  used.  In  contrast, 
the  broken  line  shows  the  result  for  the  electric  field  parallel  to  q,  with 
v_  *  0.7vri  A  substantial  change  of  the  emission  characteristics  is  obtained.  The 
spectrum  splits  into  two  broad  structures.  The  shifts  to  higher  and  lower  energies 
demonstrate  the  transition  from  normal  free  carrier  emission  to  the  Smith-Purcell 
regime  (the  low  energy  side  disappears  for  still  higher  drift  velocities) .  While 
the  lower  part  of  the  spectrum  is  out  of  the  range  of  our  detector,  the  change  on 
f  the  high  energy  side  should  be  observable  when  changing  the  orientation  of  the 

current  with  respect  to  the  grating. 


EXPERIMENT 

We  have  plasma  etched  a  periodic  grating  with  lattice  constant  L  into  the  top 
layer  of  a  GaAs/GaAlAs  heterostructure  as  shorn  in  Fig.  (2a).  The  height  of  the 
modulation  was  of  the  order  of  200  A  and  the  electron  density  is  2 .8 x 101 '/cm1 . 

The  modulation  induces  a  periodic  potential  acting  on  the  2D  electrons,  which  is 
evident  from  Shubnikov-de  Haas  (SdH)  oscillations.  At  4.2 X  the  samples  had  a  zero 
field  mobility  of  2.0 x  lO^cm’/Vs  before  plasma  etching,  and  very  narrow  minima  in 
Ojpc-  After  etching  the  mobility  was  only  reduced  to  1.5x  10*cm’/Vs,  however  the 
minima  became  much  broader.  We  interpret  this  broadening  as  the  appearance  of  a 
sinusoidal  potential,  leading  to  a  local  variation  of  the  electron  density.  From 
the  width  increase  of  the  SdB  oscillation  we  estimate  the  size  of  the  periodic 
potential  to  be  of  the  order  of  1  meV. 

Two  different  sample  geometries  were  realized!  Sample  1  had  the  grating  grooves 
parallel  to  the  contacts  (Fig.  2b),  and  as  a  result  q//j.  In  sample  2  both  geome¬ 
tries  (q//j  and  qxj)  are  realized.  The  drift  velocity  and  mobility  versus  the  ex¬ 
ternal  electric  field  as  derived  from  direct  conductivity  measurements  of  the 
small  sample  structure  (Fig.  2b)  at  4.2  X  before  and  after  the  plasma  etching  are 
shown  in  Fig.  (3) .  It  is  evident  that  the  saturation  drift  velocity  is  already  ob¬ 
tained  for  electric  fields  in  the  order  of  100  V/cm.  The  observed  dependences  of 
vQ  and  u  are  in  good  agreement  with  the  theoretical  calculation  of  Lei  et  al.  /II/, 
based  on  a  non-Boltzmann  balance  equation  approach  and  using  a  drifted  Fermi/Dirac 
distribution  with  temperature  T#. 

The  analysis  of  the  emitted  radiation  in  the  presence  of  a  periodic  potential 
allows  us  to  answer  the  question  whether  the  drift  velocities  are  realistic  and 
what  is  the  average  temperature  of  the  carrier  distribution.  Emission  spectra  for 
a  plasma  etched  sample  1  (Fig.  2b)  for  several  electric  fields  are  shown  in  Fig. (4) . 
For  the  lowest  electric  field  the  detector  response  is  only  due  to  broad  band 
thermal  emission  from  the  heated  electron  distribution  /4/.  Above  an  electric 
field  of  30  V/cm  the  emission  spectra  start  to  show  a  structure  around  a  magnetic 
field  of  0.3  T,  which  develops  to  a  clear  but  broad  peak  at  100  and  150  v/cm.  The 
spectra  show  a  clear  shift  between  50  and  100  V/cm,  indicated  by  the  arrow  which 
marks  the  decay  of  the  signal. 

For  samples  without. a  grating  the  spectra  show  a  similar  form  as  the  curve  for 
10  V/cm  for  all  electric  fields,  but  the  broad  peak  around  0.3 T  is  missing.  Con¬ 
sequently  the  additional  structure  can  be  directly  attributed  to  Smith-Purcell 
type  emission. 
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To  demonstrate  that  the  observed  emission  is  due  to  the  grating  induced  potential 
we  prepared  a  sample  where  the  grating  wavevector  q  is  oriented  both  parallel  ar.d 
normal  to  the  current  direction.  In  Fig.  (5)  a  comparison  of  the  emission  forc//j 
(dashed  curve)  and  q±j  (full  curve)  is  shown  for  an  electric  field  of  lOOV/cm.  In 
both  orientations  radiation  due  to  the  Smith-Purcell  effect  is  observed.  It  is 
evident  that  the  spectrum  for  q//j  is  shifted  to  higher  energies.  In  the  gj.j  case 
the  scattering  depends  essentially  on  the  carrier  temperature,  independent  of  the 
drift  velocity.  The  emission  comes  mainly  from  carriers  moving  with  the  Fermi 
velocity  perpendicular  to  the  current  direction. 

From  this  experiment  it  is  directly  evident  that  the  shift  of  the  spectra  is  due 
to  the  drift  of  the  carriers.  The  spectra  thus  show  the  behaviour  predicted  by  the 
theory  (Fig.  lb).  For  the  given  situation  the  shift  can  be  converted  into  the 
drift  velocities,  although  we  have  to  keep  in  mind  that  the  detector  has  a  rather 
limited  resolution  of  1  meV  and  that  the  emitted  spectrum  is  quite  broad. 

As  an  additional  proof  that  the  emission  is  due  to  the  high  carrier  velocity  we 
reduced  the  carrier  velocity  with  an  applied  magnetic  field  perpendicular  to  the 
2D-plane.  The  selective  emission  and  the  difference  in  the  spectra  between  j//q 
and  j^q  disappeared  as  shown  in  Fig.  (5)  for  a  magnetic  field  of  1.5T. 

As  a  rough  estimate  we  obtain  a  drift  velocity  of  1.5x  I07cm/s  for  lOOV/cm.  This 
drift  velocity  is  somewhat  smaller  than  that  derived  from  current  measurements. 

In  addition  we  can  estimate  the  carrier  temperature:  It  has  to  be  in  the  order  of 
40  -  50  K  from  a  fit  of  the  decaying  part  of  the  spectra  which  are  marked  by  arrows. 
Our  calculations  show  that  for  considerably  higher  temperatures  (Te  >70  X)  the 
spectra  become  extresiely  broadened  and  it  would  not  be  possible  to  observe  well- 
defined  structures  with  the  given  resolution  of  the  InSb  detector. 
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(a)  Scattering  diagrams  used  for 
the  theoretical  calculation. 

(b)  Calculated  emission  intensity 
versus  frequency  from  drifted  elec¬ 
tron  distributions  with  a  tempera¬ 
ture  of  35  K,  comparing  j*q  (full 
curve)  with  j//q  (broken  curve) . 
The  spectra  are  weighted  by  a  to 
account  roughly  for  the  detector 
sensitivity . 
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